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1.  INTRODUCTION 

The  objective  of  this  report  is  to  document  the  review  of  the  current  methods  used  to  predict  the 
system  dynamics  of  an  altered  structure  or  of  combined  structures  based  upon  a  previously  defined, 
modal  or  impedance,  model  of  the  structure(s).  Of  particular  interest  is  the  performance  of  such 
modeling  methods  with  respect  to  experimentally  based  models. 

This  report  investigates  several  system  modeling  techniques  to  determine  their  capabilities  and 
limitations  from  a  theoretical  and  practical  viewpoint.  Several  experimental  techniques,  practical 
aspects  of  the  analytical  and  approximate  techniques,  test  results,  modeling  results,  and  analysis  of 
the  results  are  presented  to  compare  and  evaluate  the  various  modeling  methods.  This  study 
presents  all  of  the  techniques  in  a  consistent  manner  from  the  same  origin,  using  consistent 
nomenclature,  to  clearly  highlight  the  similarities  and  differences  inherent  in  their  development 
which  form  the  basis  of  the  strengths  and  weaknesses  of  each  technique.  To  gain  practical 
insights,  all  of  the  techniques  presented  in  Sections  2  through  4  are  compared  with  experimental 
results. 

1.1  System  Modeling 

System  modeling  is  a  computer  based  technique  that  is  used  to  represent  the  dynamic 
characteristics  of  a  structure.  This  representation  takes  the  form  of  either  experimental  data,  modal 
data,  or  analytical  data.  Once  the  dynamic  characteristics  of  a  structure  are  used  to  form  a  model  or 
system  model,  several  uses  of  the  model  are  possible.  First,  the  effects  of  design  changes  or 
hardware  changes  to  the  original  structure  could  be  studied.  Second,  the  structure  could  be  coupled 
with  another  structure  to  determine  the  overall  resultant  dynamic  behavior.  Finally,  the  model  can 
be  used  analytically  to  apply  forces  and  determine  the  forced  response  characteristics  of  the 
structure. 

The  main  objective  of  system  modeling  is  to  use  a  mathematical  representation  of  the 
dynamic  characteristics  of  a  structure  in  a  computer  environment  to  effectively  develop  a  design 
or  trouble  shoot  a  particular  problem  of  a  design.  Several  techniques  of  various  origin  have 
evolved  with  the  advancement  of  computer  technology.  Depending  on  the  situation,  each  is  very 
effective  if  properly  utilized.  Design  development  is  generally  considered  an  extensive  long  range 
process  that  results  in  an  optimally  designed  structure  given  the  constraints  of  the  project. 
Trouble  shooting  involves  the  evaluation  of  failures  or  design  flaws  which  must  be  corrected  quickly. 

The  most  obvious  way  to  classify  system  models  is  into  analytical  and  experimental  methods. 
The  primary  approach  to  analytical  modeling  is  commonly  known  as  finite  element  Analysis.  Finite 
element  analysis  is  analytical  in  nature  because  only  knowledge  of  the  physical  properties  of  the 
structure  is  used  to  build  a  dynamic  representation.  This  is  done  by  subdividing  the  structure  into 
discrete  elements  and  assembling  the  linear  second  order  differential  equations  by  estimating  the 
mass,  stiffness,  and  damping  matrices  from  the  physical  coordinates,  material  properties,  and 
geometric  properties. 

Finite  Element  Analysis  is  extremely  useful  because  no  physical  test  object  is  necessary  to  compute 
resonant  frequencies  and  mode  shapes,  forced  response  simulation,  or  hardware  modifications. 
Therefore,  this  method  is  extremely  useful  in  the  development  cycle  where  it  can  be  used  to  correct 
major  flaws  in  the  dynamic  characteristics  of  a  structure  before  a  prototype  is  built.  Since  finite 
element  analysis  is  an  approximate  analytical  technique,  experimental  modal  analysis  data  is 
obtained  as  soon  as  possible  so  the  analytical  model  can  be  validated.  Detailed  finite  element 


-1- 


analysis  method  is  not  covered  in  this  report. 

Experimental  modeling  techniques  are  further  subdivided  into  two  groups.  They  are  modal  models 
and  impedance  models.  Modal  modeling  is  an  experimentally  based  technique  that  uses  the 
results  of  an  experimental  modal  analysis  to  create  a  dynamic  model  based  on  the  estimated  poles 
of  the  system.  The  result  is  a  computationally  efficient  model  that  is  uncoupled  due  to 
transformation  of  the  physical  coordinates  into  the  modal  coordinates  which  renders  the 
system  into  a  number  of  lumped  mass,  single  degree  of  freedom  components.  This  concept  is 
fundamental  to  modal  analysis.  This  model  is  then  used  to  investigate  hardware  changes,  couple 
structures,  or  compute  the  forced  response.  This  method  is  developed  fully  in  Section  2. 

Impedance  modeling  uses  measured  impedance  functions  or  frequency  response  functions  to 
represent  the  dynamic  characteristics  of  a  structure  in  the  frequency  domain.  This  method  uses 
the  experimentally  measured  functions  to  compute  the  effects  of  hardware  changes  or  to  couple 
together  several  components.  To  use  this  modeling  technique,  measurements  at  the  constraint  or 
connection  points,  driving  point,  and  cross  measurements  between  the  two  are  needed  to 
compute  a  modified  frequency  response  function.  Impedance  modeling  is  fully  developed  for  the 
compliance  method  and  stiffness  method  in  Section  4. 

Both  of  the  experimental  modeling  methods  are  quick  and  easy  to  use  in  their  basic 
implementations.  Therefore,  they  are  extremely  useful  in  trouble  shooting  situations  but  have 
limited  application  in  the  design  cycle. 

The  final  classification  of  system  models  is  experimental/analytical  models  or  mixture  methods.  Two 
techniques  are  considered  mixture  methods.  The  first  is  sensitivity  analysis.  Sensitivity  analysis  is 
an  approximate  technique  that  uses  the  first  term  or  first  two  terms  of  a  power  series  expansion  to 
determine  the  rate  of  change  of  eigenvalues  or  eigenvectors  with  respect  to  physical  changes  (mass, 
stiffness,  and  damping).  Therefore,  this  method  is  used  for  trend  analysis,  selection  of  hardware 
modification  location,  and  design  optimization.  This  technique  is  considered  a  mixture  method 
because  it  computes  sensitivity  values  for  modal  parameters  which  result  from  either 
experimental  modal  analysis  or  finite  element  analysis. 

The  other  experimental/analytical  method  is  component  mode  synthesis.  Component  mode 
synthesis  and  the  building  block  approach  are  techniques  that  use  experimental  or  analytical  modal 
representations  of  the  components  of  a  large  or  complex  structure  to  predict  the  resultant  dynamic 
characteristics  of  the  entire  structure.  Furthermore,  this  technique  has  evolved  to  the  point  where 
components  are  combined  in  either  physical  or  modal  coordinates.  Therefore,  this  technique  is 
truly  a  mixture  method  where  experimental  and  analytical  data  are  used  to  optimize  a  design.  This 
method  is  very  useful  in  the  design  cycle  of  industries  that  produce  large  structures,  such  as  the 
automotive  and  aerospace  industries.  A  new  component  mode  synthesis  method  (not  a  classical 
mode  synthesis  method  such  as  SYSTAN)  is  discussed  in  detail  in  Section  5  and  the  building  block 
approach  is  discussed  in  Section  4  along  with  the  development  of  the  impedance  modeling  technique. 

1.2  Practical  Application 

To  determine  the  practical  limitations  and  capabilities  of  the  various  techniques,  each  is 
compared  with  experimental  results  for  four  hardware  changes  to  an  H  frame  of  hollow  tubing, 
welded  construction,  and  rigid  end  plates.  The  four  design  changes  considered  are  a  simple 
stiffener  addition,  addition  of  two  masses,  beam  addition,  and  plate  addition.  The  modeling 
techniques  were  used  to  predict  new  dynamic  characteristics  only  when  the  changes  were  applicable 
within  the  constraints  of  their  intended  implementation. 


An  experimental  modal  analysis  was  conducted  on  the  H  frame  in  its  original  state.  The  frequency 
response  functions  and  resulting  modal  parameters  are  inputs  to  the  various  programs  to  predict  the 


effects  of  the  design  changes. 


In  summary,  system  modeling  is  a  diverse  field  that  involves  many  aspects  of  structural  dynamics. 
One  of  the  primary  goals  of  this  report  is  to  present  this  material  in  a  concise  and  consistent 
manner  to  reduce  unnecessary  confusion  and  better  relate  the  various  factions  involved. 
Furthermore,  each  technique  is  applied  to  a  structure  to  gain  further  insight  into  the  practical  aspects 
of  system  modeling. 

1.3  Boundary  Conditions 

In  the  application  of  the  system  modeling  techniques  mentioned  above,  there  are  three  test 
configurations  used  to  obtain  the  frequency  response  functions  or  the  derived  modal  data  base. 
These  three  test  configurations  involve  the  boundary  condition  and  can  be  summarized  as  free-free, 
constrained,  and  actual  operating  boundary  conditions.  In  terms  of  analytical  modeling  technique, 
such  as  finite  element  analysis,  various  boundary  conditions  can  be  easily  simulated  in  the 
mathematical  model  to  predict  the  dynamic  characteristics  of  a  system.  Therefore,  this  gives  the 
analysts  the  luxury  to  evaluate  the  model  using  desired  boundary  condition.  In  terms  of  experimental 
modeling  technique,  in  a  laboratory  environment,  usually  it  is  very  difficult  and  too  costly  to 
implement  the  test  fixture  to  simulate  the  actual  operating  condition  of  a  complete  system,  or  the 
constrained  boundary  condition  of  a  component,  'nierefore,  most  of  the  modal  tests  are  performed 
under  an  environment  simulating  free-free  boundary  conditions. 

Besides  the  boundary  conditions  mentioned  above,  modal  tests  can  be  performed  on  the  mass-loaded 
structures  to  predict  the  shifted  dynamic  characteristics  of  the  original  structure.  The  advantages  of 
adding  lumped  masses  at  the  connecting  or  attachment  points  of  a  component  under  testing  are:  (1) 
modal  coefficients  associated  with  those  connecting  degrees  of  freedom  can  be  more  accurately 
excited  and  described  under  the  mass  loading  effect,  (2)  rotational  degrees  of  freedoms  at  the 
connecting  degrees  of  freedom  can  be  computed  using  rigid  body  computer  programs  if  sufficient 
number  of  accelerometers  are  mounted  on  the  additive  masses,  (3)  analytically,  added  masses  can  be 
removed  from  the  mass-loaded  testing  configuration,  and  the  enhanced  modal  parameters  of  the 
original  structure  can  be  obtained,  (4)  if  the  dynamic  characteristics  of  the  original  structures  are 
available  through  analytical  or  experimental  method,  then  more  accurate  generalized  masses  can  be 
obtained  through  these  two  sets  of  data.  Examples  of  applying  the  mass  additive  technique  can  be 
found  in  References  [1]  and  [2], 


2.  MODAL  MODELING 


2.1  Modal  Modeling  Overview 

This  section  develops  the  system  modeling  technique  known  as  modal  modeling.  Modal  modeling  is 
also  known  as  the  Snyder  Technique,  Local  Eigenvalue  Modification,  Structural  Modifications, 
Dual  Modal  Space  Structural  Modification  Method,  and  Structural  Dynamic  Modification  t316l. 
The  common  thread  of  the  research  mentioned  is  that  all  utilize  a  model  in  generalized  or  modal 
coordinates  from  experimental  data  upon  which  to  investigate  structural  changes.  Structural 
changes  are  transformed  into  modal  coordinates  and  added  to  the  structure  and  the  result  is  resolved 
to  yield  the  modified  modal  parameters. 

The  modal  modeling  technique  was  initially  published  by  Kron  ^  in  1962  and  extended  by 
Weissenberger  Simpson  and  Taborrok  ^  ,  and  Hallquist,  Pomazal  and  Snyder  Early 
researchers  in  this  area  restricted  themselves  to  a  local  modification  eigensolution  technique. 
Several  software  packages  have  been  developed  employing  this  technique  since  the  implementation 
and  widespread  use  of  digital  Fourier  analysis.  Notably,  Structural  Measurement  Systems,  Inc.  first 
released  Structural  Dynamics  Modification  ^  using  the  local  modification  procedure.  The  local 
modification  procedure  allows  only  simple  mass,  stiffness  and  damping  changes  between  two 
general  points. 

Recent  research  has  progressed  in  several  areas.  First,  Hallquist  and  Snyder  Luk  and  Mitchell  ^ 
and  SMS  t7'  have  used  the  local  modification  technique  for  coupling  two  or  more  structures 
together.  Hallquist  l101,  O’Callahan  ^nl  and  De  Landseer  ^121  have  expanded  the  method  to  include 
complex  modes.  O’Callahan  ^13'  and  SMS  ^  have  found  methods  to  approximate  more  realistic 
modifications  using  local  modification  for  trusses  and  beams.  Mitchell  and  Elliot  t15'  and 
O’Callahan  and  Chou  have  developed  different  methods  that  use  full  six  degree  of  freedom 
representations  that  depend  on  the  approximation  or  measurement  of  the  rotational  degrees  of 
freedom  to  make  beam  or  plate  modifications. 

2.2  Modal  Model  Development 

The  equations  of  motion  of  a  mechanical  structure  are  second  order  linear  differential  equations. 
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where  M,  C  and  K  are  the  structure’s  mass,  damping,  and  stiffness  matrices  of  order  N  by  N.  The  ,-l 

displacement  vector  {x}  is  in  physical  coordinates,  {f}  is  the  applied  forcing  vector,  and  N  is  the  r 

number  of  physical  degrees  of  freedom.  j-1 


Two  techniques  are  commonly  used  to  find  the  solution  vector  {x}.  Either  a  solution  is  assumed 
of  the  form  {x}  =  {Ar}e*‘ or  Laplace  transforms  are  used  to  arrive  at  a  set  of  linear  equations. 

[[M]52  +  [C].v  +  [^]]W)}={^)}  (2) 

where  the  initial  conditions  are  equal  to  zero 
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Now  the  substitution 
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[BOO]  =  [[A/]s2  ♦  [C]5  +  [X]]  (3) 

is  made  to  Equation  (2)  yielding 

[fl(s)]{*(s)}={F(i)}  (4) 

[B(s)]  is  the  system  matrix.  The  eigenvalues  of  the  structure  are  found  by  taking  the  determinant  of 
the  system  matrix  for  the  nontrivial  case  when  {X(s)}  exists.  The  vector  {X}  that  is  associated  with 
each  eigenvalue  is  the  eigenvector  found  from  the  solution  of: 

[fl(s*)]{*(s*)}  =  {0}  (5) 

The  inverse  of  the  system  matrix  B(s)  is  called  the  transfer  function  matrix  of  the  structure. 


[H(s)  ]  =  [B(s)l 


Adj\B  ] 
Dct[B] 


where  Adj  [B]  is  the  adjoint  matrix.  The  adjoint  matrix  and  determinant  of  [B]  matrix  are 
polynomials  in  the  variables  because  the  elements  are  functions  of  s.  Therefore,  the  transfer 
function  matrix  is  conveniently  written  in  partial  fraction  form. 


IH(s>1=ii  <£>*<£>/ 


where  the  \  is  the  poles  or  zeroes  of  the  determinant  matrix  and  N  is  the  number  of  degrees  of 
freedom  included  in  the  model.  Continuous  systems  have  infinite  degrees  of  freedom;  therefore,  N 
modes  form  a  truncated  modal  model  in  practice. 

The  elements  of  the  system  matrix  are  quadratic  in  s;  therefore,  the  number  of  poles  p  is  2N. 
Furthermore,  since  M,  K,  and  C  are  real  matrices  the  poles  are  in  conjugate  pairs  and  are  purely 
imaginary.  Also,  the  residue  matrix  [A  ],  is  an  N  by  N  matrix  of  complex  values.  This  concept  is 
critical  to  experimental  modal  analysis  because  the  rows  and  columns  of  [A],  are  made  of  a  linear 
combination  of  the  homogeneous  solution  vector, 
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where  Qr  is  a  scaling  constant  and  {<fi}r  =  {.Y}  is  the  mode  shape  vector  of  mode  r. 

The  physical  significance  of  this  is  that  experimentally,  frequency  response  function  is  measured  in 
the  complex  plane  at  s  =  jw.  When  frequency  response  functions  are  measured,  either  rows  or 
columns  of  the  residue  matrix  are  defined.  Thus,  the  result  of  an  experimental  modal  analysis 
parameter  estimation  is  the  eigenvalues  (frequency  and  damping)  and  eigenvectors  (mode  shapes) 
from  the  frequency  response  functions.  Due  to  the  unique  orthogonality  relationships,  modal 
vectors  are  linearly  independent  quantities,  which  is  stated  mathematically  as: 
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(8) 


{*}.r[KM*=0  (9) 

for  i  t  k. 

Thus,  the  projection  of  one  vector  on  another  is  zero  and  is  orthogonal  with  respect  to  the  [M]  and 
[K]  matrices.  The  orthogonality  relationships  allow  a  set  of  simultaneous  equations  to  be 
transformed  into  uncoupled  equations;  this  is  perhaps  the  single  biggest  advantage  of  the  modal 
domain.  An  N  degree  of  freedom  system  is  orthogonal  in  N  space.  It  is  desired  to  transform  our 
system  in  physical  coordinates  into  modal  or  generalized  coordinates.  If  i  =  k  =  r  in  Equation  (8) 
and  (9). 


{*} rTWMr=K 

(10) 

{4>}rT[KMr=K 

(ID 

Mr  and  K,  are  scalar  quantities  related  to  the  modal  mass  and  modal  stiffness  of  mode  r. 
Magnitudes  of  mr  and  K,  are  dependent  on  the  modal  scaling  factor  Qr. 

The  coordinate  transformation  was  required  to  diagonalize  the  physical  mass  and  stiffness  matrices. 
As  already  noted  in  Equation  (8)  and  (9),  when  the  mass  or  stiffness  matrix  is  pre-multiplied  by 
the  transpose  and  post-multiplied  by  different  eigenvectors  the  result  is  zero.  If  this  same 
operation  is  performed  with  the  same  vectors,  a  constant  results.  Therefore,  if  all  the  modal 
vectors  are  formed  into  a  matrix  it  is  readily  seen  that  the  desired  transformation  matrix  exists  in 
the  modal  matrix  [$]. 

2.3  Local  Eigenvalue  Modification 

The  local  eigenvalue  procedure  is  a  simple  technique  that  transforms  modifications  into  modal 
space.  This  is  done  for  simple  structural  changes  hence,  a  very  efficient  solution  is  obtained. 
Going  back  to  the  original  equation  of  motion  and  assuming  proportional  damping: 


[Afp}*[CK*W*P}  =  {/} 


(12) 


Now,  this  set  of  equations  is  transformed  into  modal  space  by  substituting  the  coordinate 
transformation  {x}  =  (Op)  and  pre*  multiplying  by  the  transpose  of  the  modal  matrix  [#]: 


(13) 

f*]T[C][*)  =  rCj 

(14) 
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l*]T[*][*]  =  % 


(15) 


This  is  a  statement  of  orthogonality  with  respect  to  the  properties  of  the  modal  vectors. 

The  set  of  equations  in  modal  coordinates  is  written  as: 

fMj  { q }  ♦  fCj  {q }  ♦  %  {q}  =  [*f{/}  (16) 


Note  that  the  mass,  stiffness,  and  damping  matrices  are  assumed  to  be  diagonal.  For  this  case,  the 
system  is  classically  or  proportionally  damped,  and  the  set  of  equations  yields  a  real  normal  mode 
model.  Therefore,  this  development  does  not  treat  the  nonproportionally  damped  case. 


To  make  modifications  that  have  physical  significance,  the  mode  shapes  must  be  scaled  by  some 
technique  so  that  the  equations  can  be  consistently  transformed  into  modal  space  to  yield 

accurate  results.  Unity  modal  mass  scaling  is  used  by  choosing  Qr  =  -  in  Equation  (7)  to  put 

2/h. 


the  equations  in  a  convenient  form. 


[/]{<?}  ♦  r2 sAfiq)  ♦  U’j {q}  =  [*]{/} 


(17) 


where: 


tyfly  -  modal  damping  matrix 

r(\2j  -  modal  stiffness  matrix 

<7,  -  damping  coefficient  (  =  sfy) 

(, r  -  damping  ratio  for  mode  r 

(X  -  undamped  natural  frequency  of  mode  r 

Now  the  orthogonality  equations  are  written  as: 

[*f[M][*]  =  [/] 

(18) 

[*]r[C][$]  = 

(19) 

[*]Tmr*]  = 

(20) 

The  results  of  a  modal  test  are  the  eigenvalues  (frequency  and  damping): 


which  are  solutions  of  the  characteristic  equation: 


det  [[Af]s2  +  [C]s  ♦  [X]]  =  0  (22) 

and  the  eigenvectors  are: 

[*]  =  . {fis]  (23) 

which  are  solved  from  the  homogeneous  equation: 

[[Af]  A,2+[C]  V[Xj]  {fir  =  0  r  =  1, - N  (24) 

Using  the  simple  relationships  (Equation  (12)  through  (24)),  modifications  can  be  made  to  a 
structure. 

[[MMAM]]  {i‘}+|[C]+[AC]J  {ih[[X]+[AX]J  {*}={/}  (25) 

where: 

[AM]  -  mass  modification 
[AC]  -  damping  modification 
[AX]  -  stiffness  modification 

This  equation  is  transformed  into  modal  coordinates  using  the  standard  transformation  [4>]  : 

wm  ♦  [CKq}  ♦  [XKq}  =  imf)  (26) 

where: 

[M]  =  [/]  +  [$]r[AM][$] 

[C]  =  r2fta,  ♦  [*]T[AC][*] 

[X]  =  fa3j  ♦  [*F[AX][$] 
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Note  that  the  bar  over  a  variable  denotes  that  the  matrix  is  modified. 

The  eigenvalues  and  eigenvectors  of  the  modified  system  are  found  as  follows.  The  eigenvalues 
are  calculated  from  the  determinant  equation : 


H 


det  [[Mis3  +  (C]s  +  [/q]  =  0 


(27) 


The  eigenvectors  are  calculated  as  a  solution  of  the  homogeneous  equation: 


[[M]\3 +  [<:]>, +  [X]|{?}r  =  0  r  =  l . N 


This  results  in: 


[$]  =  [$][{?}l,{?}2.  . . . 


(28) 
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Therefore,  given  a  set  of  calibrated  frequency  response  functions  from  which  system  parameters 
are  extracted  and  scaled  to  unity  modal  mass,  numerous  modifications  can  be  investigated.  Mass, 
stiffness,  and  damping  modifications  may  be  made  between  any  two  measured  degrees  of 
freedom  or  between  ground  and  a  measured  degree  of  freedom. 

Frequently,  many  authors  have  presented  the  local  modification  procedure  similar  to  this 
development.  To  better  understand  the  local  modification  procedure,  the  form  of  the 
modification  will  be  considered  more  carefully.  Only  a  stiffness  modification  is  handled  to  simplify 
the  equations. 

Consider  two  modal  spaces  1  and  2;  1  is  the  original  modal  coordinates  and  2  is  the  modified  modal 
coordinates.  A  stiffness  modification  in  modal  coordinates  ^  is  written  as: 


[AXl2]  =  [*lr[AX][$]  (29) 


where  12  denotes  the  stiffness  change  and  may  be  represented  in  physical  form  as: 


[A/C12]  =[T]  rDj[rf 


(30) 
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where: 


[T  ]  -  The  matrix  containing  the  eigenvectors  of  [AX]  which  represent  the 
tying  vectors  of  the  physical  degree-of-freedoms 


a 
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rDj  -  The  eigenvalue  spectral  matrix  of  [AX) 
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Equation  (30)  can  be  written  as  a  summation  over  the  number  of  modifications  in  the  system. 

[ak13]  =  £{/}*</*{/}? 

k  =  l 

where: 

{f  }*  -  The  k*  tie  vector  of  [T] 
dk- The  k *  diagonal  element  of  j 
k  -  The  modification  tying  index 
M  -  The  number  of  modification(s) 

For  a  single  modification,  substituting  this  result  into  Equation  (29)  yields: 

[a k12] = < 

Now,  a  modal  projection  vector  which  indicates  the  modes  affected  by  the  modification  is  defined. 

This  identity  is  used  in  Equation  (31)  to  further  simplify  the  expression. 

[A I 


Equation  (33)  is  added  to  the  undamped  and  homogeneous  form  of  Equation  (17)  to  yield: 

M{9i2}*[faV[^«]]{^}  =  0 

The  vector  ^  ia  describes  the  change  from  modal  space  1  to  2. 

The  characteristic  Equation  of  (34)  is: 

(raV“*2K<7i  a}--{v}*rf*{v}ZW  (35) 


ten)  =  -  [  fa2j 


(35a) 


where: 

fi,  -  Natural  frequency  of  the  r*  mode  of  the  original  structure  in  the  modal  space  1 
u>2  -  Modified  natural  frequency 
multiply  both  sides  of  Equation  (35a)  by  {v}*  : 

Mlten}  =  -{v>n  (36) 

since  {v}*  ten}  t  0,  therefore  : 

i  =  -4t{v}nraaJ-‘‘*2]'1te}*  (37) 

which  is  manipulated  into  the  summation  form: 


Equation  (37)  is  the  modified  characteristic  polynomial  equation  in  modal  space  2  relative  to  modal 
space  1.  Although  the  use  of  tie  vectors  confuses  the  development,  it  lends  much  insight  to  the 
understanding  of  the  relationships  between  the  two  modal  spaces.  Once  the  eigenvalues  are 
determined,  the  eigenvectors  of  the  modified  system  are  found  using  standard 
eigenvalue/eigenvector  decomposition  algorithms  to  yield  the  eigenvectors  in  modal  coordinates. 
They  are  scaled  to  unity  modal  mass  in  relative  space: 

and  transformed  back  into  modal  space  by  multiplying  by  [$13]. 

[<*2]  =  [*i][*i2]  (38) 


The  local  eigenvalue  modification  procedure  allows  the  projection  of  the  physical  space  into 
modal  space  1  using  the  tie  vector  approach  originated  by  Kron  ^3L  Therefore,  a  complete 
eigenvalue  solution  is  not  necessary.  The  largest  advantage  of  this  method  is  that  it  is  very  quick 
because  a  full  rank  eigenvalue/eigenvector  extraction  is  not  needed. 
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2.4  Coupling  of  Structures  Using  Eigenvalue  Modification 


Modal  synthesis  is  used  primarily  to  find  the  resultant  dynamic  characteristics  of  a  large  structure 
by  adding  smaller  components  together.  It  is  assumed  that  the  components  have  been  measured 
and  data  exists  at  the  connection  points.  The  modal  synthesis  technique  has  been  implemented 
using  the  local  eigenvalue  modification  technique.  This  implementation  is  extremely  simple 
because  it  eliminates  the  need  for  attachment  modes,  residuals,  and  constraints  that  are  discussed 
in  Sections  4  and  5.  Only  the  mass,  stiffness,  and  damping  at  the  coupling  points  are  necessary, 
along  with  the  modal  data  of  the  two  structures. 

Recalling  the  original  differential  Equations  (12),  it  can  be  readily  seen  that  for  n  components,  the 
equation  becomes: 


mx}+[CKi}  +  mx}={h 


(39) 


where: 
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The  total  size  of  this  system  of  equations  is  the  sum  of  the  size  of  each  of  the  component 
equations. 


Now,  the  modifications  can  be  made  by  adding  the  AA/,  A K,  or  AC  terms  to  the  stacked  matrices 
(39).  This  results  in  the  following  equation: 


[M+AA/K*}  ♦  [C+ACKi}  ♦  [K+A/Cp}  =  {/}  (40) 


This  equation  is  transformed  into  modal  coordinates  as  before: 

lMX9}  +  09}}  +  l*XtfM*im}  (41) 

where: 

[M]  =  [M]  *  &F[AA/]&] 

[C]  =  [C]  +  [^flAC]^] 

£]  =  [K]  +  [^F[A^][4] 

fl^i]  [0]  [0]  . . .  [0] 

[0]  [*,]  [0]  . . .  [0] 

H-  :  :  ::::: 

[o]  [6]  to i'.'.Vi*. 

Assuming  that  the  structure  is  proportionally  damped,  the  equations  of  motion  of  the  stacked 
matrices  take  the  form  of  Equation  (40).  Therefore,  the  resultant  dynamics  depend  only  on  the 
modal  parameters  of  each  structure  plus  the  stiffness,  damping,  and  mass  of  the  connections. 

The  modal  parameters  of  the  two  structures  are  combined  by  the  coupling  system  transformed 
into  modal  space  12.  The  eigenvalues  A,  =  aT  +  ;uv  are  computed  from  the  determinant  equation: 

det  [(A/Js2  ♦  [C]s  ♦  [*]]  =  0  (42) 

and  the  eigenvectors: 

[$j)  =  [$1][{^}1,{^}2,...,{^}r]i2] 

are  found  from: 

[[A/]X3MClA;M^]]{?i2}r  =  0  r  =  l . n  (43) 


Once  again,  the  addition  of  the  connections  [A/C]  is  made  using  a  tie  vector  approach, 


Once  again,  the  addition  of  the  connections  [AKJ  is  made  using  a  tie  vector  approach. 


[A/C]  =  [r]  tojfrF  (44) 


where  the  tie  vector  {t}  designates  the  coordinates  to  be  joined.  Using  the  modal  transformation  [<£] 
the  new  eigenvalue  problem  of  the  total  modified  system  is  formulated.  The  tie  vector  provides  the 
connection  or  tie  between  the  original  modal  space  and  the  modified  modal  space.  Once  again,  the 
characteristic  equation  is  Equation  (37)  and  the  eigenvectors  are  the  vectors  associated  with  the 
nontrivial  solution  of  the  eigenvalue  problem.  Once  they  are  obtained  in  the  modified  modal 
space  12  they  are  transformed  into  modal  space  2  using  Equation  (38). 

2.5  Complex  Mode  Eigenvalue  Modification 

Thus  far,  only  real  normal  modes  have  been  considered.  This  approximation  is  very  accurate  for 
most  physically  realizable  structures.  Some  systems  have  heavier  damping  or  exhibit 
nonproportional  damping.  In  this  limited  case  complex  modes  are  needed  to  adequately  describe  the 
modal  characteristics.  Another  situation  may  arise  when  complex  modes  are  necessary.  When 
investigating  a  modification  that  tends  to  make  the  damping  nonproportional  or  increase  the 
damping  ratio,  a  real  normal  mode  method  will  not  accurately  model  the  structure.  Complex 
modes  provide  an  exact  solution  but  make  the  development  and  understanding  of  the 
technique  more  difficult. 

For  the  nonproportional  damped  case  the  undamped  normal  modes  do  not  uncouple  the 
equations  of  motion.  Frazer,  Duncan,  and  Collar  proposed  a  method  to  reduce  the 

nonproportional  damped  second  order  differential  equation  of  motion: 


{y}>»xl  ~ 


Now  the  transformation  from  physical  coordinates  y  to  modal  coordinates  q  of  rank  2m  by  2N  is: 


{y}>nxi  -  [^lanx2w{9} 


(48) 


i  The  state  vector  velocity  transformation  is  found  by  taking  the  derivative  with  respect  to  time: 

;  {y}»»xi  =[^]rAj{9}  (49) 

where: 

!Aj  -  The  diagonal  spectral  matrix  of  rank  2N  by  2N 
i  of  the  original  complex  eigenvalues  A, 

1  N  -  The  number  of  degree  of  freedom  in  modal  coordinates 

m  -  The  number  of  degree  of  freedom  in  physical  coordinates 
Applying  the  standard  coordinate  transformation  to  the  homogeneous  form  of  Equation  (46) 

[Wm*]  {?}  =  (Wmw  rAj{<7>  (so) 

Since  original  equations  are  symmetric  and  real,  the  right  and  left  hand  eigenvectors  are  equal. 

Using  the  orthogonality  conditions  the  equations  become  diagonal  which  is  the  advantage  of  the 
|  U  and  V  formulation.  U  unity  scaling  method  is  used  by  letting  Qr  =  1.0  to  make  Equation  (50) 

I  valid : 


\muim  =  rAj  (5D 

=  in  w 


The  modal 

as: 


matrix  ['I']  may  be  partitioned  with  respect  to  the  state  vector  defined  in  Equation  (48) 


[*1 


W4 

fV] 


(53) 


where  [♦]  is  the  modal  matrix  of  rank  (m  by  2N)  used  to  transform  physical  coordinates  into 
modal  coordinates  for  the  complex  case. 


(54) 


Once  again,  since  the  necessary  coordinate  transformations  exist,  a  structural  modification  or 
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coupling  of  systems  can  be  made  using  a  local  eigenvalue  technique.  Changes  are  represented  as 
AM,  AC  and  A K.  The  original  modal  space  is  modal  space  1  and  the  modified  modal  space  is  modal 
space  2.  Since  the  equations  are  reformulated  into  U  and  V  matrices  the  appropriate  A U  and  AK 
must  be  found. 

Modifications  in  the  original  system  of  Equation  (45)  appear  as: 

[[MM AM]]  {i>[[C]+[AC]]  {i}*[[Kj+[AK]]  {*}  =  {/}  (55) 


Taking  advantage  of  the  transformation  for  the  homogeneous  solution: 

where: 


The  [Af/12]  and  [AF12]  matrices  are  the  modifications  that  transform  the  system  from  modal 
['  space  1  to  modal  space  2.  Utilizing  Equations  (50)  through  (54)  on  Equation  (55)  produces: 

1 

[rAij  *  [^nAC12]f^|]{r/l}  =  [[/]  ♦  [^F[AU12][^l]  rAlJ{9l}  (57) 

» 

. 

r  This  equation  can  be  solved  to  produce  the  modified  eigenvalues  and  eigenvectors  in  space  12 

t  which  maps  space  1  into  space  2.  Therefore,  a  transformation  exists: 


3 


3 


» 


! 

i 


{<7i}  =  f’i'iaKfla} 

where:  [^12]  is  the  modal  matrix  of  rank  2N  by  2N  found  from  Equation  (56) 

The  eigenvectors  are  normalized  to  unity  [fl12]  matrix  and  transformed  into  modified  space  2  by: 

[^2]  =  fap'u] 

Thus,  [^2 1  is  the  new  modal  matrix  of  the  modified  structure.  It  was  obtained  by  projecting 
information  into  the  new  modal  space  which  eliminates  the  need  for  a  new  eigenvalue  solution. 

2.6  Realistic  Modal  Modifications 

Several  researchers  have  developed  techniques  to  make  more  realistic  modifications.  Rarely  is  a 
physical  structural  change  realistically  modeled  as  a  local  modification  between  two  points.  This  has 
led  to  the  development  of  truss  and  beam  type  modifications  in  local  modification  schemes. 
Another  major  problem  in  the  prediction  of  realistic  structural  modifications  is  the  lack  of 
rotational  degrees  of  freedom  in  a  modal  model.  Yasuda  and  Brown  Martinez  O’Callahan 
and  Lieu  ^  ,  arid  Smiley  and  Brinkman  have  proposed  various  methods  to  measure  or  predict 
the  rotational  frequency  response  at  a  point. 
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The  details  of  these  techniques  are  beyond  the  scope  of  this  investigation  so  they  are  not  fully 
developed.  Structural  Measurement  Systems  has  proposed  making  a  series  of  local  eigenvalue 
modifications  to  approximate  a  realistic  modification.  For  a  rib  element  a  three  point  technique 
is  used  to  make  three  modifications  locally  to  approximate  a  rib.  This  technique  does  not  require 
rotational  degrees  of  freedom. 

O’Callahan  and  Chou  t13]  have  proposed  two  techniques  to  make  beam  modifications.  The  first  is 
a  local  eigenvalue  modification  technique  that  requires  data  for  the  rotational  degrees  of 
freedom  through  measurement  or  calculation  at  the  modification  points.  A  three  dimensional 
symmetric  beam  stiffness  element  is  approximated  by  making  six  individual  modifications.  This 
is  accomplished  by  adding  the  axial,  shear,  and  bending  stiffness  to  each  of  the  two  connection  points. 

The  second  technique  is  a  general  technique  that  introduces  the  concept  of  group  eigenvalue 
modifications.  Rather  than  using  only  a  single  tie  vector  to  make  the  modification,  a  group  of  tie 
vectors  is  developed  for  the  stiffness  matrix  to  transform:  normalized  coordinates  into  the 
principal  axes,  principal  bending  axes  into  a  centroidal  set  of  coordinates,  centroidal  axes  to  a 
local  reference  through  the  centroid,  local  coordinates  from  the  centroid  to  a  geometric  reference 
point,  local  bending  axes  at  the  reference  point  to  a  parallel  set  through  the  attachment,  and  local 
coordinates  of  the  beam  element  into  global  coordinates. 

Also,  a  global  mass  tie  matrix  is  derived  by  using  only  the  last  three  tie  matrices  mentioned 
above.  Six  spectral  coefficients  are  needed  to  formulate  the  final  matrix.  Three  are  the 
translational  mass,  one  the  torsional  inertia,  and  the  last  two  are  rotary  inertia.  The  only  remaining 
information  needed  to  develop  the  stiffness  tie  matrix  is  the  six  spectral  coefficients  which  define 
axial  stiffness,  torsional  stiffness,  shear  about  the  two  principal  axes,  and  moments  about  the  two 
principal  axes.  This  formulation  requires  rotational  information  at  the  modification  points. 
This  formulation  results  in  the  general  mass  and  stiffness  matrices  for  an  offset  beam. 

Elliot  and  Mitchell  ^  have  proposed  an  alternate  technique  to  transfer  the  mass  and  stiffness 
matrix  of  a  beam  into  a  dynamic  stiffness  matrix.  Therefore,  the  modal  model  is  also  formed  using 
the  dynamic  stiffness. 

]  -  u?\M  ]j  =  Dynamic  Stiffness  Matrix  (58) 

The  beam  transfer  matrix  is  more  flexible  than  the  previous  method  because  a  number  of 
structural  elements  already  exist  in  this  form.  This  technique  also  transforms  the  elements  in  a 
group  sense  which  makes  it  very  efficient. 

The  need  for  more  realistic  modifications  has  led  to  research  in  the  area  of  beam  modifications 
using  modal  modeling  techniques.  The  main  drawback  to  these  techniques  is  that  rotational 
degrees  of  freedom  are  required  to  properly  formulate  the  beam  dynamic  matrix.  Currently,  several 
techniques  are  being  researched  to  provide  the  necessary  rotational  information  experimentally. 
Until  a  technique  becomes  commercially  available,  beam  modifications  will  not  be  a  useful  tool 
in  modal  modeling. 

2.7  Normalization  of  Measured  Complex  Modes 


2.7.1  Introduction  and  Background 

Experimental  modal  analysis  is  carried  out  to  extract  a  set  of  modal  parameters  from  the  measured 
time  or  frequency  domain  data  of  the  structure  under  test.  The  identified  eigenvectors  are  in  general 
complex  modes  due  to  several  possible  reasons: 

1.  The  damping  is  nonproportional,  i.e.,  [C]  matrix  is  not  proportional  to  [K]  and  [M]  matrices. 

2.  Measurement  errors  due  to  mass  loading  effects,  noise,  nonlinearities  etc. 

3.  Digital  processing  errors  due  to  finite  frequency  resolution,  leakages,  high  modal  density,  and 
frequency  response  functions  estimation  procedure  (Z/^//^//,). 

4.  Modal  parameter  estimation  errors  due  to  invalid  estimation  of  number  of  degrees  of  freedom. 

The  identified  complex  modes  can  be  used  directly  in  the  applications  of  modal  modeling,  structural 
dynamic  modifications  .sensitivity  analysis,  and  validation  and  optimization  of  an  analytical  model. 
On  the  other  hand,  real  normal  modes  are  sometimes  more  desirable  in  those  applications  due  to  the 
facts  that  (i)  normal  modes  are  numerically  easier  to  handle  than  complex  modes  and  (ii)  analysts 
usually  compute  normal  modes  rather  complex  modes  in  the  finite  element  analysis.  Therefore,  if  the 
normal  modes  are  desired,  a  real-normalization  procedure  is  needed  to  be  implemented.  Some  of  the 
proposed  normalization  methods  are  summarized  in  the  following. 

2.7.2  Normalization  Using  Direct  Parameter  Estimation  Technique 

The  Direct  Parameter  Estimation  technique  ^  is  a  multi-degree-of-freedom  frequency  domain 
parameter  estimation  algorithm.  This  technique  manipulates  multiple  response  functions  from  a 
single  reference  location  to  obtain  global  least-square  estimates  of  the  modal  properties. 

Consider  the  equation  of  motion  in  the  frequency  domain, 

-w?[A/K*}  *  [Kpf.MF.}  (59) 

where: 

<4  =  frequency  i 
[M  ]  =  NpxNp  mass  matrix 
[C  ]  =  NpxNp  viscous  damping  matrix 
(A  ]  =  NpxNp  stiffness  matrix 
{A",}  =  vector  of  Np  responses  at  frequency  i 
{ F,}  =  vector  of  Np  forces  at  frequency  i 
Np  =  number  of  physical  degrees  of  freedom 

For  the  single  reference  case,  the  force  vector  would  be  made  up  of  all  zeroes  except  at  the  reference 
location.  Premultiplying  Equation  (59)  by  the  inverse  of  the  mass  matrix. 


(60) 


where: 


M 

> 

i 
♦ 


[a]  =  [Af  ]"*[C] 

m=ml[K] 

{w}  =  -[A/]-l{F} 


Reordering  Equation  (60): 


7hMA'.}  +  [^K^}+{v}  =  w?{a;} 


Dropping  vector  and  matrix  notation  for  simplicity: 


[<*  0  v]  -  X,  =  a%2*, 


D,=  Xi 
1 


and  E,  =  J}  X{ 


Consider  k  spectral  lines  at  w  =  wx  a^, _ ,w*, 


[*Pv]  [Dl\D2\....  |  Dk]  =  [Ex  |  £a  |  ....  |Et] 


Equation  (64)  can  be  solved  using  a  least  square  technique.  The  system  matrices  a  and  0  can  be  used 
to  solve  for  the  eigenvalues  and  eigenvectors  of  the  system  using  a  companion  matrix  approach. 


The  companion  matrix  for  the  general  case  of  nonproportional  damping: 


[7  A 


where  [1]  is  the  identify  matrix. 


Since  the  a  matrix  contains  all  the  damping  information,  and  the  0  matrix  contains  the  mass  and 
stiffness  information,  the  companion  matrix  can  also  be  defined  to  solve  for  the  undamped  natural 
frequencies  and  a  set  of  normal  modal  vectors  of  the  system.  For  this  case,  the  a  matrix  is  neglected 
and  the  resulting  companion  matrix  is  shown  below: 


[/#fl 


An  estimate  of  the  proportional  damping  can  be  obtained  using  the  left  and  right  eigenvectors, 
resulting  from  this  solution,  to  diagonalize  the  damping  matrix.  Using  the  diagonal  elements  and  the 
natural  frequency  of  the  system,  the  estimates  of  proportional  damping  can  be  obtained. 


rirnyw 
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Principal  Response  Analysis  t24'  is  applied  to  improve  the  numerical  efficiency  of  the  algorithm  by 
allowing  a  large  number  of  response  functions  to  be  used  in  the  analysis. 

2.7.3  Ibrahim’s  Time  Domain  Technique 

S.R.  Ibrahim  indicated  I25,26!  that  large  errors  may  result  from  simplified  normal  mode 
approximations  to  complex  modes. 


fasRi*  jli 


where: 

fa  =  Ith  element  of  a  complex  modal  vector 
R,  =  real  part  of  fa 
Ii  -  imaginary  part  of  fa 

The  approximated  normal  mode  element  fa  corresponding  to  fa  is: 


where: 

fa  =  ith  element  of  the  normal  modal  vector 

The  ±  sign  in  Equation  (66)  is  determined  by  the  phase  angle  of  the  the  complex  modal  coefficient 
fa.  For  example,  if  a  normal  modal  vector  has  zero  phase  angle,  then  fa  has  "+”  sign  if  the  phase  angle 
of  the  complex  fa  is  between  ±  90°. 

A  time  domain  technique  has  been  proposed  to  compute  a  set  of  normal  modes  from  the  identified 
complex  modes.  The  required  data  are  a  set  of  modal  parameters  such  as  may  be  identified  from  a 

modal  survey  test.  These  modal  parameters  are  namely  a  set  of  complex  modes  {V>}r,  r=l, . N  and  a 

set  of  corresponding  roots  Ar,  r=l, — ,N(  and  their  complex  conjugates).  The  modal  vectors  have  Np 
elements  where  Np>  >N.  To  compute  the  normal  modes,  two  approaches  were  suggested  in  the 
following  way: 

2.7.3.1  Approach  1:  Using  an  Oversized  Mathematical  Model 

From  the  given  modal  parameters,  displacement,  velocity,  and  acceleration  responses  are  formed 
according  to  the  equations. 


«o}=  mu'*  +{«i(0} 


V-.-V 
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V'  O  v*«  *  »  *  -  ‘r.  "  -  «  .  . 


{i(0}=  ♦{MO} 


{i(0}=  ZXMre*’*  ♦{MO} 


where  MO.MO.MO  are  added  random  noise  of  uniform  distribution.  These  responses  are  then 
used  in  the  state  vector  equation. 


Ji(OlJ  [0] 
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[K}\  -[[A/rl[C]] 
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M  =  [£Ky} 


where  {y}  is  now  system’s  state  vector  containing  the  displacements  and  velocity  responses.  By 
repeating  Equation  (71)  2NP  time  instants,  the  following  equation  is  satisfied: 


[y]  =  [£]tv] 


where  [y  ]  and  [y]  contain  responses  measured  at  the  2NP  time  instants.  From  Equation  (72)  the  [E] 
matrix  can  be  identified  as. 


[^]  =  tv][y]-] 


By  computing  matrix  [E],  the  [A/]'x[X]  matrix  gives  normal  modes  according  to  the  eigenvalue 
equation, 


[[A/Wl]  {*)  =<*?{*} 


Naturally,  without  any  noise,  the  matrix  [y]  is  singular  since  the  number  of  degrees  of  freedom  is 
larger  than  the  number  of  modes  present  in  the  responses.  A  small  amount  of  noise  makes  the 
inversion  of  [y]  possible  for  the  purpose  of  extracting  modal  information.  The  mechanism  on  which 
the  success  of  this  approach  is  based  is  explained  in  Reference 

2.7.32  Approach  2:  Using  Assumed  Modes 

The  given  set  of  complex  modal  parameters  satisfy  the  equation 


snsw 
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[[A/H*]  [A#m]{  ** 

» 


=  {-A?W}  fr=l . N) 


(75) 


Since  there  are  only  N  modes  and  the  system  has  Np  degrees  of  freedom.  Equation  (75)  cannot  be 
solved  for  [[A/]*1!*]  [A/  J-1  [C  ] j .  Let  us  assume  that  there  exists  a  set  of  vectors  Yk  and  a  set  of 
characteristic  roots  sk ,  k=  N+l,  N+2, — ,/Vp.  This  set  of  assumed  parameters  are  selected  such  that 


A,  P  sk 


YkP  [V»i  V>2  ••••  V’jvH'J} 


where  {a}  is  any  vector  of  coefficients.  This  implies  that  Yk  and  {\p}r  for  all  r  and  k  from  a  linearly 
independent  set  of  vectors.  In  such  case,  it  can  be  written  that 

[[A/m  [A*rl[C]j  |Jtyytj={-s2n}  (k=N+\,N+2 . Np)  (76) 

and  Equations  (75)  and  (76)  can  be  solved  for  [[A/]*1  [AC]  [A/]‘1[C]|  from  which  normal  modes  are 
computed  according  to  Equation  (74). 

2.7.4  Frequency  Domain  Normalization  Technique  using  orthogonality  constraints 

Zhang  and  Lallement  ^  proposed  a  numerical  method  to  determine  the  real  eigenvalues  and 
eigenvectors  of  the  associated  undamped  structure  from  the  identified  complex  eigensolutions. 

2.7.4.1  Recall  of  General  Relations 

Under  autonomous  conditions,  the  dynamic  behavior  of  a  mechanical  system  is  assumed  to  be 
represented  by  a  second  order  linear  differential  equation. 


[A/  K*(0}  +  [C  K*(0}  ♦  [*  K*  (0}  =  0 


(77) 


where  [M],[C],[K]  are  NpxNp  positive  definite  real  matrices. 

The  solution  of  Equation  (77)  is  of  the  form:  {x  (r)}  =  {X}e*.  This  leads  to  the  eigenvalue  problem. 


[ A2[Af  ]  ♦  A[C  ]  +  [K  ]]{</>}  =  0 


(78) 


Assuming  the  case  of  nonproportional  damping,  the  2 Np  eigensolutions  ( Ar,  (V'Jr)  are  then  found  by 
conjugate  pairs.  Introducing  the  spectral  matrix  fAj  and  the  modal  matrix  [4>]  defined  by 


r  _  A]fv*p 

“  Ta*. 


[^LvpjOJVp  -  [[^LvpXATp  I  [*  LvpXWp] 


the  first  N„  relations  of  Equation  (78)  are  collected  in  the  form: 


[Af][*Kx*p  rAJ  +  [X]['tUpXWp=0 


The  transformation  of  Equation  (77)  into  the  space  of  state  variable  y(t)  gives: 


where: 


■"HS1  P 


The  solutions  of  Equation  (82)  are  also  of  the  form:  {y(0}  ={>'}eA*.  This  leads  to  the  eigenvalue 
problem: 


[in*A[t/]]  {y}  =  0. 


Introducing  the  modal  matrix  where  [Y]  ==  [ _ ,{y}r, _ ],  and  it  can  be  shown  directly 

that  [  Y]  is  of  the  form: 


1v>-[(a*i  ia1Vi]  and  W--{A; 


and  that  the  matrices  [U],  [V],  rAj  and  [Y]  satisfy  the  orthogonality  relations: 


[Yf  [t/JM  •  r*j;  [Yr[V][Y]=  r/?i  TAj 


(86) 


where  ^  ^  =  diagonal{rT)  = 

eigenvector  {y}r. 


jNp^ 


J**p»»*p 


,  and  rT  characterizes  the  norm  chosen  for  the 


Considering  the  partition  into  submatrices  of  (VMY],  and  fAj,  the  two  following  equations  are 
directly  deduced  from  Equation  (86): 


rAj[*f  [M][V]  fAj  -  [*ffK][V]  =  fAj  (NpxNp) 


(87) 


rA*j  [¥*  f  [M  ][*]  rAj  -  [*f  [K  ][*]  =  0  (A rpxNp) 


(88) 


The  conservative  system  associated  with  Equation  (77)  is  finally  characterized  by: 


[MHx(0MKHx(0>  =  0 


(89) 


It  is  known  that  its  solutions  of  the  following  form:  {x  (f)}  =  {X)e±jut,  yield  the  eigenvalue  problem: 


[[A]-n?[M]]  {^}  =  0 


(90) 


Collecting  these  m  eigenvalues  (f^.,{V>}r)  in  the  spectral  matrix  ffy  =  diagonal^},  and  in  the  modal 
matrix  [$]  =  [  — — 1  ,  these  matrices  satisfy  the  orthonormality  relations: 


[*f[A/][*]  =  [/]  ;  [*m]I*]-raj 


(91) 


2.7 .4.2  Data  and  Hypothesis 


The  problem  considered  here  has  the  following  assumptions: 
1.  The  matrices  [M],[C],[K)  are  unknown. 


2.  From  dynamic  tests  carried  out  on  the  structure,  N  (N<  <Nf)  eigensoiutions  are  identified  by 
means  of  Nm  (N  <Nm<Np)  measurement  points  located  on  the  structure.  These  N 
eigensoiutions  are  characterized  by  the  spectral  submatrix:  fA^  =  diagonal{  A,},  and  the 
modal  matrix: 

3.  corresponding  submatrices  of  the  associated  conservative  system,  are  to  be 
determined  from  fAj,[^]. 

The  hypotheses  to  solve  the  unknowns  are  described  as  follows: 

Hypothesis  1:  It  is  assumed  that  the  identified  eigenvectors  can  be  represented  by  a  linear 
combination  of  the  corresponding  normal  eigenvectors  of  the  associated  conservative  systems: 


[ip]  =  [$]  [r],  where:  [r]  is  a  complex  NyN  matrix 


(92) 


Let: 


[ip]  =  [ipk  +  /[ip]^  (NmxN)  where:  re  =  real  part,  im  =  imaginary  part 


[r]  =  [rk+;[ru  (n*n) 


(93) 


fAj=  (NyN) 


Hypothesis  2:  The  matrices  [¥]  and  [<&]  are  assumed  to  have  rank  N.  Therefore,  matrix  [r]  is 
nonsingular.  The  matrix  [r]«  is  also  assumed  to  be  nonsingular. 

Hypothesis  1  leads  to  the  following  relationship: 


I 


[*]  =  [*L  [ri™1 

[rU  =  frjre  [IV]  where:  [W\  is  a  real  NxN  matrix. 

The  combination  of  Equations  (92)  and  (94b)  yields: 

[*u  =  i*u  m 


(94a) 

£ 

% 

(94b) 
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The  best  estimation  of  [W],  by  means  of  the  least  squares  method,  is  then 


According  to  Hypotheses  1  and  Equation  (94b),  it  can  be  then  written: 


[*]  =  [*]  in,  [[/Lv^+ym] 


(97) 


[*1  =  [*W  [Ttt 


(98a) 


The  matrix  [♦]„  being  known,  the  problem  is  now  to  determine  the  matrix  [r]«  first,  then  the  matrix 
raj- 


The  solution  procedure,  therefore,  is  briefly  described  as  follows: 

—  By  substituting  Equation  (92)  and  (97)  into  Equation  (87)  and  (88),  and  then  separating  the  real 
and  imaginary  parts,  expressions  of  the  symmetrical  matrices  [r]£  [1%  and  [r]«  rfVj  [r]„  can  be 
derived.  These  matrices  can  be  expressed  as  functions  of  matrices  [$],  fAj  and  [W], 

—  Calculate  matrices  (T]£  (rU  and  [r]£  ffy  [r)«. 

—  Matrices  [r]«  and  ify  are  determined  by  solving  the  following  eigenvalue  problems: 


[[fii|-MB2]HM  =  0.  {vr}T[fl3Kvr}  =  1,  r- 1,---. at 


(98b) 


where  [flj  =  [T]l  ^  [rU  and  [02]  =  [r£  [ V U- 
and 


Uj  =diagonal{\},  and  [V]  =  [  •  •  •  {vr}  •  •  • )  =  [1%* 


(98c) 


2.7.5  Normalization  Using  Improved  Computational  Model 

Natke  and  Rotert  I28'  proposed  to  obtain  a  set  of  normal  modes  by  improving  an  existing 
computational  model  using  measured  eigenquantities  (incomplete  and  erroneous)  in  a  subsystem 
formulation.  This  improved  computational  model  is  then  used  to  compute  the  eigenquantities  of  the 
associated  undamped  system.  Since  this  method  assumes  the  use  of  a  analytical  or  finite  element 
model,  detailed  discussion  of  this  technique  will  not  be  presented  here. 


2.7.6  Time  Domain  Technique  using  Principal  Response  Analysis 


a 

a 


Nl 


A  time-domain  technique  is  presented  to  compute  a  set  of  normal  modes  from  the  measured 
complex  modes.  The  number  of  degrees  of  freedom  which  is  equal  to  the  number  of  measurements, 
in  general,  is  much  larger  than  the  number  of  measured  modes.  By  using  the  proposed  method  ,  a 
large  number  of  physical  coordinates  are  reduced  to  a  smaller  number  of  principal  coordinates  using 
principal  response  analysis  technique.  From  the  given  modal  parameters,  free  decay  responses  are 
calculated  using  properly  scaled  complex  modal  vectors.  The  companion  matrix  for  the  general  case 
of  nonproportional  damping  is  also  derived  in  the  principal  coordinates.  Normal  modes  are  obtained 
through  eigenvalue  solutions  of  the  [M  and  transformed  back  to  the  physical  coordinates  to 
get  a  set  of  normalized  real  modes.  A  numerical  example  is  presented  to  support  the  outlined  theory. 


> 
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2.7.6. 1  Theory  and  Formulation:  computation  of  normal  modes 


Assume  there  exists  a  set  of  measured  modal  parameters  from  a  modal  test.  These  modal  parameters 

consists  of  N  complex  modes(  and  their  complex  conjugates),  Ar,  r=l . N  and  {ip}„  r=l, . N. 

Each  modal  vector  has  dimensions  Nm  and  in  general  Nm>  >N.  To  compute  a  set  of  normal  modes 
from  the  given  complex  modes,  a  time-domain  approach  is  formulated  as  follows: 

Free  decay  displacement,  velocity  and  acceleration  responses  can  be  expressed  as: 

aw  .  . 

{*(0)-  (99) 
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The  equations  of  motion  for  the  general  case  of  nonproportional  damping  is: 


Another  form  of  the  homogeneous  Equation  (102)  can  be  written  as: 
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{>(/)}  =  [£]  {>(/)} 


where  {y(t)}  is  now  the  system’s  state  vector  containing  the  displacement  and  velocity  responses.  In 
the  physical  coordinates,  dropping  vector  notation.  Equation  (104)  can  also  be  written  as: 


Mi) . y('™Jl  =  l£][y(M . rt'aO) 


(105) 


Therefore,  [E]  matrix  can  be  computed  as: 


I£]  =  tv]  trr1 


(106) 


where  [E],  [y]  and  [y]  are  all  2Nm-x2Nm  matrices. 

By  computing  the  matrix  [E],  [A/]_1[K]  matrix  gives  the  normal  modes  according  to  the  eigenvalue 
solution: 


[[AfWjl*]-**!*! 


(107) 


The  matrix  [y]  is  always  singular  due  to  the  number  of  degrees  of  freedom  Nm  is  larger  than  the 
number  of  measured  modes  N.  Therefore,  it  is  numerically  difficult  and  unstable  to  solve  [E]  matrix 
in  the  physical  coordinates.  Therefore,  it  is  proposed  to  compute  the  companion  matrix  [E]  in  the 
principal  coordinates  in  the  following  way: 

Define  [<X>]  =  modulus  matrix  of  complex  modal  vector  matrix  [»],  i.e.,  [$>]  is  real  modal  vector 
approximation  according  to  Equations  (65)  and  (66),  and 


m  =  [*]  m 


(108) 


where  [W|  is  a  NxN  complex  matrix  and  can  be  obtained  through  pseudo  inverse  technique. 

Using  singular  value  decomposition  technique  and  assume  [$]  has  full  rank  N  (which  is  true  if  an 
independent  set  of  modal  vectors  is  used),  [$]  can  be  decomposed  as: 


m-wtziisr 


(109) 


where: 


[P]  =  orthonormal  matrix  ( NmxN ) 

r  Ej  =  diagonal  matrix  consists  of  eigenvalues  of  (NxN) 

[5  ]  =  unitary  matrix  consists  of  eigenvectors  of  [$]T[$]  (NxN) 

By  using  matrix  [P]  as  the  transformation  matrix,  physical  coordinates  can  be  transformed  to  principal 
coordinates  as  follows: 


-28- 


{p(/)}  =  [Pf{x(0} 

substitute  Equations  (99)  and  (108)  into  (110),  we  get 

=  [pfmm  r^j 
=  [/>mr£j  [5rmr^j 
=  rEjisn^jr^j 


Similarly, 


{/>(0}  =  [/>]7WAV‘J  (113) 

Therefore,  free  decay  displacement,  velocity  and  acceleration  responses  can  be  evaluated  using 
Equations  (111)  through  (113)  in  the  principal  coordinates.  Once  {q(f)}  and  {<j(f)}  are  computed  at 
2N  time  points  ,  companion  matrix  [E]  (2N  by  2N)  can  be  computed  in  the  principal  coordinates 
accordingly  as: 


[£]  =  l9][9]-1 


where:  {<7  (0}  =  K  (<7(0} 


Since  [qj  is  a  full  rank  matrix,  it  is  always  invertible  and  the  computation  is  numerically  stable  and 
accurate. 

Based  on  the  result  of  matrix  [E],  the  N  by  N  [Af  ]_1  [/C  ]  matrix  gives  the  normal  mode  solutions  [$]  in 
the  principal  coordinates  according  to  Equation  (106). 

A  set  of  normal  modes, [$]  in  the  physical  coordinates  can  be  obtained  through  the  coordinate 
transformation: 

[*1  =  [/’1(*J  (H5) 

2.1.62  Discussion  and  Conclusion 

From  the  numerical  example  of  a  nonproportionally  damped  discrete  system  described  later  in 
Section  6.2,  it  indicates  that  a  set  of  real-normalized  modal  vectors  can  be  obtained  using  this 
proposed  method.  Nevertheless,  unless  a  complete  set  of  complex  modal  vectors  are  included  in  the 
analysis,  analysis  results  are  subject  to  modal  truncation  errors  to  certain  extent. 
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Finite  element  modes  can  be  added  to  the  measured  modal  data  set  outside  the  measurement  band 
to  improve  the  analysis  results.  The  analysis  modes  can  be  real  normal  modes  with  proportional 
damping  or  complex  modes  with  nonproportional  damping. 

For  a  very  lightly  damped  system  the  magnitude  of  the  M'lC  matrix  is  in  general,  much  smaller  than 
the  magnitude  of  the  M'lK  matrix.  Based  on  past  experiences,  the  real-normalized  modal  vectors  will 
be  similar  to  those  normal  modes  by  taking  the  magnitudes  of  the  complex  modes. 

The  true  free  decay  responses  of  a  system  are  dependent  on  the  initial  conditions  of  the  structure,  but 
the  modal  properties  of  a  system  are  independent  of  the  initial  conditions.  With  the  proposed 
method,  it  is  recommended  that  all  modal  vectors  included  in  the  analysis  band  be  scaled  to  the  same 
order  of  magnitude  before  computing  the  free  decay  responses  according  to  Equations  (99)  to  (101). 
And  by  doing  so,  the  diagonal  matrix  1  £j  in  Equation  (109)  will  always  have  rank  N  and  make  the 
coordinate  transformation  feasible  in  the  application  of  principal  response  analysis. 

2.8  Errors  in  the  Experimental  Model 

Since  the  structural  dynamic  modification  method  uses  the  modal  parameters  of  the  unmodified 
structure  to  predict  the  effects  of  the  modification  ,  it  is  apparent  that  the  accuracy  of  the  results  is 
dependent  on  the  validity  of  the  data  base  supplied.  For  this  reason,  it  is  important  that  all  errors 
associated  with  the  data  acquisition  and  processing  of  the  unmodified  structure  are  minimized. 

Among  the  sources  of  error  that  must  be  addressed  are  nonlinearity,  standard  FFT  errors  (aliasing 
and  leakage),  scaling  errors,  as  well  as  truncation  errors  of  the  modal  model  itself.  Only  the  most 
serious  of  these  errors  will  be  addressed  at  this  time. 

One  of  the  assumptions  in  the  experimental  modal  analysis  is  linearity  of  the  system.  Therefore,  all 
methods  discussed  here  are  based  on  this  assumption.  In  reality,  inaccurate  results  will  arise  when 
linear  system  coupling  algorithm  is  used  to  predict  the  dynamic  characteristics  of  structure(s)  coupled 
together  by  nonlinear  joints. 

Leakage  is  a  measurement  error  that  arises  from  the  processing  of  signals  that  are  not  periodic  in  the 
time  window  of  the  signal  analyzer.  Because  of  this  truncation  in  the  time  domain,  the  Fourier 
coefficients  of  the  sampled  signal  do  not  lie  on  the  A/ of  the  analyzer.  This  causes  energy  at  a  specific 
frequency  to  spread  out  into  adjacent  frequency  bands,  and  results  in  an  amplitude  distortion  at  the 
actual  frequency.  It  is  this  amplitude  error  that  causes  scaling  errors  in  the  modal  mass  and  stiffness 
estimates  associated  with  each  mode.  This  in  turn  affects  the  modification  process  by  scaling  the 
predicted  modes  of  the  modified  structure  by  an  amount  that  is  proportional  to  the  amplitude  error 
of  the  original  data. 

This  leakage  problem  can  be  reduced  by:  (1)  using  periodic  excitation,  or  transient  excitation  signals 
such  as  burst  random,  or,  (2)  using  the  Ht  frequency  response  function  procedure,  or,  (3)  using  cyclic 
averaging  technique  in  the  measurement  stage,  or,  (4)  taking  data  with  smaller  A/,  or,  (5)  using 
windows  on  the  time  domain  measurement,  such  as,  impact  with  exponential  window  applied  to  the 
response  signals.  If  exponential  window  is  used,  it  must  be  accounted  for  in  the  parameter  estimati  on 
because  it  adds  artificial  damping  to  the  structure.  Care  must  be  taken  not  to  overcompensate  for 
this  damping  allowing  the  poles  of  the  system  to  become  negatively  damped. 

Modification  errors  often  arise  from  using  a  modal  data  base  that  is  not  properly  scaled  relative  to 
the  system  of  units  used  by  the  modeling  software.  This  type  of  error  will  occur  if  the  modal 
parameters  are  estimated  using  frequency  response  functions  which  were  measured  using  improper 
transducer  calibrations.  This  improper  scaling  once  again  results  in  improper  estimates  of  modal 


mass  and  stiffness  used  by  the  modeling  software. 

Up  to  this  point,  it  has  been  shown  that  the  accuracy  of  the  results  is  dependent  on  the  amount  of 
error  in  the  data  base.  Another  concern  that  needs  to  be  addressed  is  the  validity  of  the  modal 
model.  Because  the  effects  of  a  structural  modification  are  calculated  in  modal  space,  if  an 
insufficient  number  of  modes  are  included  in  the  original  data  base,  there  will  be  a  limit  to  the 
number  of  modal  vectors  that  can  be  predicted.  This  phenomena  is  known  as  modal  truncation  and 
should  be  considered  in  choosing  a  frequency  range  for  the  analysis.  It  may  be  desirable  in  some 
cases  to  extend  the  frequency  range  of  data  acquisition  above  the  actual  frequency  range  of  interest 
for  the  data  base  to  include  a  few  extra  modes.  This  extended  frequency  range  will  improve  the 
calculation  of  out  of  band  residuals,  and  may  help  for  the  case  where  these  out  of  band  modes  are 
shifted  into  the  frequency  range  of  interest  by  the  modification.  Care  must  be  taken  in  extending  the 
frequency  range,  to  prevent  an  excessive  loss  of  frequency  resolution. 

Another  concern  in  the  development  of  the  modal  model  is  the  number  of  degrees  of  freedom  to 
include  in  the  analysis.  By  definition,  the  number  of  degrees  of  freedom  must  be  equal  to  or  greater 
than  the  number  of  modes  in  the  frequency  range  of  interest.  Realistically,  the  number  of  DOF 
should  be  much  larger  than  the  number  of  modes  of  interest  to  accurately  define  the  individual 
modal  vectors.  Since  modal  coefficients  exist  only  at  points  where  data  has  been  taken,  it  is  possible 
to  miss  nodal  lines  of  the  structure  if  too  few  points  are  included  in  the  analysis.  This  generally 
becomes  a  more  significant  problem  for  higher  order  mode  shapes. 

Once  the  original  data  base  has  been  established,  the  modal  parameters  can  be  estimated  using  any  of 
several  existing  parameter  estimation  algorithms.  All  of  these  algorithms  attempt  to  yield  a  best 
estimate  of  the  actual  parameters.  Because  there  will  always  be  some  degree  of  experimental  error  in 
the  data,  the  resulting  estimates  of  modal  parameters  will  be  subject  to  error.  In  order  to  minimize 
this  error,  it  is  advantageous  to  use  some  sort  of  least  squares  implementation  to  yield  a  best  estimate 
of  modal  parameters. 

The  use  of  SDOF  versus  MDOF  parameter  estimation  algorithms  is  determined  by  the  modal  density 
of  the  structure  being  analyzed.  If  a  SDOF  method  is  used  for  a  structure  with  closely  coupled 
modes,  poor  estimates  of  modal  mass  and  stiffness  are  obtained,  and  the  modification  routine  will 
yield  poor  results. 

The  various  errors  mentioned  in  the  previous  paragraphs  are  commonly  committed,  and  easily 
overlooked  when  performing  a  modal  test.  This  is  not  intended  to  be  an  exhaustive  list  of  errors 
affecting  modal  modeling,  but  an  indication  of  the  types  of  things  that  must  be  kept  in  mind  when 
establishing  a  valid  modal  data  base.  Without  good  estimates  of  the  original  structures  modal 
parameters,  there  can  be  no  serious  attempt  at  accurately  predicting  the  characteristics  of  the 
modified  structure. 

2.9  Validation  of  Experimental  Modal  Model 

As  mentioned  in  the  previous  section,  the  accuracy  of  the  modal  modeling  or  structural  modification 
results  is  dependent  on  the  validity  of  the  modal  model  supplied.  There  will  always  be  some  degrees 
)f  experimental  and  modal  parameter  estimation  errors  in  the  data  base.  Therefore,  it  is  important 
for  the  users  to  qualitatively,  and  if  possible,  quantitatively,  examine  the  validity  and  errors  of  the 
modal  model  before  it  is  used  to  predict  the  system  modeling  or  modification  results.  Although 
perfect  results  should  not  be  desired  in  the  application  of  modal  modeling  technique,  it  is  important 
to  realize  that  any  modal  model  obtained  experimentally  is  far  from  being  perfect.  Therefore,  it  is 
suggested  that  the  experimental  modal  model  be  validated,  or  optimized  ,  before  the  model  is  input 
to  any  modal  modeling  algorithm.  In  this  section,  some  of  the  validation  methods  are  briefly 


summarized  as  follows: 


2.9.1  Frequency  Response  Function  Synthesis 

Synthesizing  a  frequency  response  function  (not  used  in  the  estimation  of  modal  parameters)  using 
extracted  modal  parameters  and  compared  with  a  measured  frequency  response  function  at  the 
synthesized  measurement  degree  of  freedom  is,  in  general,  a  common  practice  during  the  modal 
parameter  estimation  process.  If  a  good  match  exists  between  these  two  sets  of  frequency  response 
functions,  then  it  is  a  good  indication  that  the  extracted  modal  parameters  are  agreeable  with  the 
measurement  data.  But  this  doe  s  not  guarantee  that  there  is  no  error  exist  in  the  modal  data  base. 

2.9.2  MAC 

Modal  Assurance  Criteria  (MAC)  ^  is  commonly  used  to  check  the  consistency  of  the  extracted 
modal  parameters,  when  more  than  one  estimate  of  each  mode  is  available. 

2.9.3  Detection  of  Mode  Overcomplexity 

This  method  qualifies  each  mode  by  a  number  called  the  Mode  Overcomplexity  Value  (M.O.V.)  and 
the  global  Modal  Model  by  the  Mode  Overcomplexity  Ratio  (M.O.R.)  The  basic  idea  of  the 
Mode  Overcomplexity  test  is  that,  for  good  modal  models  with  complex  modes,  the  frequency 
sensitivity  for  an  added  mass  change  should  be  negative.  If  it  happens  that  the  sensitivity  is  positive, 
it  is  caused  by  either  an  incorrect  scale  factor  (modal  mass)  or  by  the  fact  that  the  phase  angle  of  the 
complex  modes  compared  to  the  normal  mode  phase  angle  exceeds  a  certain  limit;  in  other  words,  it 
is  due  to  an  overcomplexity  of  the  mode  shape. 

The  MOV  is  defined  as  the  ratio  of  the  number  of  positive  frequency  sensitivities  over  the  number  of 
all  the  frequency  sensitivities  for  a  particular  mode.  To  give  more  weight  to  points  with  a  high  modal 
displacement  compared  to  points  with  a  small  modal  displacement,  a  weighted  sum  is  introduced  to 
give  a  more  general  evaluation  of  the  modal  model.  The  value  of  MOV  is  between  1  and  0,  the 
bigger  the  value  is,  the  modal  model  is  more  overcomplex. 

m  m 

The  MOR  is  defined  as  the  ratio  of  over  (  1  -  )  which  gives  a  one  figure  assessment 

*=1  1=1 

of  the  modal  model  with  respect  to  its  overcomplexity.  The  MOR  ranges  from  zero  to  infinity.  A 
low  MOR  value  indicates  good  modal  data,  while  a  large  MOR  indicates  a  scale  factor  problem  or  a 
overcomplexity  problem. 

2.9.4  Mass  Additive/Removal  technique 

This  technique  employs  a  mass  additive  or  removal  procedure  to  verify  or  validate  experimental 
modal  model  in  the  application  of  modal  modeling  technique.  Modal  model  can  be  obtained  from 
either  the  original  structure,  or,  mass-loaded  structure  W.  If  modal  model  is  obtained  from  the 
original  structural  configuration,  then,  comparisons  can  be  made  between  the  analytically  predicted 
dynamic  characteristics  of  the  mass-loaded  structure  from  the  modal  modeling  algorithm,  and  the 
test  results  (such  as  modified  resonant  frequencies)  obtained  from  the  physically  modified  mass- 
loaded  structure.  If  there  is  no  good  agreement  between  these  two  sets  of  results  for  the  mass-loaded 
structure,  then  this  is  an  indication  that  global  or  local  scaling  errors,  or  overcomplexity  of  some 
measured  complex  modes,  existing  in  the  experimental  modal  model.  If  high  quality  data  are  desired 


in  predicting  the  system  dynamics  of  the  altered  structure  or  combined  structure(s),  then  the 
previously  determined  experimental  modal  model  needs  to  be  validated,  if  possible,  or,  a  new  set  of 
data  needs  to  be  recollected  before  any  modeling  application  is  attempted. 

For  the  second  case,  i.e.,  if  modal  model  is  obtained  from  the  physically  mass-loaded  structure,  then 
comparisons  can  be  made  between  the  analytically  predicted  and  experimentally  measured  dynamic 
characteristics  of  the  mass-removed  structure.  Similarly,  if  there  is  no  good  agreement  between  these 
two  sets  of  data,  this  indicates  some  errors  existing  in  the  original  modal  model.  In  Reference  [2], 
using  approximated  real  modes  from  the  measured  complex  modes,  a  modal  scaling  procedure  can  be 
used  to  correct  the  global  scaling  errors  in  the  experimental  modal  model. 

The  number  of  masses  and  the  size/weight  of  each  additive  mass,  can  be  added  to  the  structure  is 
dependent  on  the  total  mass  and  size  of  the  structure(s).  In  general,  the  following  rules  can  be  used 
as  guidelines  in  considering  the  number  and  size(s)  of  the  additive  mass(es): 

•  The  added  mass(es)  can  be  considered  rigid  in  the  frequency  range  of  interest. 

•  With  small  amount  of  mass(es)  added  or  removed  to  or  from  the  structure,  the  mode  shapes  can 
be  considered  unchanged  before  and  after  the  modification. 

•  Sensitivity  of  the  change  of  system  dynamics  is  dependent  on  the  location(s)  of  the  added 
mass(es).  In  other  words,  if  there  is  only  one  mass  added  to  a  large  structure,  then  some  of  the 
modes  may  not  be  sensitive  enough  to  alter  their  frequencies  due  to  the  fact  that  the  added  single 
mass  is  near  the  nodal  points  of  such  modes. 

•  Rotational  degrees  freedom,  if  permitted,  can  be  extracted  from  the  rigid  body  motion  of  the 
lumped  mass(es).  This  information  is  very  useful  if  the  mass  mounting  point(s)  is(are)  the 
connection  or  coupling  point(s)  of  the  structure(s). 

2.9.5  Improvement  of  Norms  of  Modal  Vectors 

Zhang  and  Lallement  proposed  a  method  to  improve  the  norms  of  the  measured  modal  vectors 
and  then  calculate  the  generalized  modal  masses  of  the  original  structure.  This  method  will  correct 
modal  scaling  errors  in  the  modal  data  base.  This  method  requires  a  set  of  modal  data  from  the  initial 
structure  and  a  set  of  data  from  the  mass  loaded  (perturbated)  structure. 

2.9.6  Conclusion 

There  are  many  other  techniques  currently  being  studied  in  order  to  evaluate  and  improve  the 
accuracy  of  the  experimental  modal  model.  Some  involve  the  use  of  a  finite  element  model.  The 
validation  methods  described  in  this  section  serve  under  the  assumption  that  the  experimental  modal 
model  needs  to  be  verified,  and  furthermore,  validated  through  an  acceptable  engineering  procedure. 
This  procedure,  which  is  needed  to  guarantee  the  accuracy  of  the  modal  modeling  results  within  an 
acceptable  limit,  has  been  proven  to  be  a  missing  link  between  the  modal  data  base  and  the  analytical 
model  modeling  technique. 

2.10  Summary 

In  summary,  modal  modeling  has  been  discussed  from  its  inception  by  Kron  through  present 
day  research  on  beam  modifications  in  the  modal  domain.  Modal  modeling  is  a  technique  that  is 
very  quick  because  the  generalized  coordinates  have  a  reduced  number  of  degrees  of  freedom. 
Therefore,  many  modifications  can  be  investigated  in  a  short  time.  Earlier,  One  this  technique 


was  presented  mainly  as  a  trouble  shooting  technique.  In  fact,  researchers  ^  have  found  this 
method  to  be  three  to  six  times  faster  than  analytical  techniques.  This  ratio  increases  with  the  size  of 
the  problem.  The  speed  of  this  technique  and  its  interactive  implementation  make  it  well-suited  for 
on-site  problem  solving  and  initial  design  cycle  work. 

Many  limitations  are  apparent  from  the  development.  First,  if  experimental  data  is  used,  the 
frequency  response  functions  must  be  carefully  calibrated.  This  technique  is  extremely  sensitive 
to  experimental  errors  in  general.  Data  must  be  carefully  acquired  to  avoid  bias  errors  such  as 
leakage  and  aliasing.  Errors  made  in  the  estimation  of  the  frequency  response  functions  translate 
into  errors  in  the  modal  model  and  modal  matrix  [$].  Modal  parameter  estimation  is  extremely 
critical  in  modal  modeling.  Parameter  estimation  is  a  two-stage  process  that  estimates  eigenvalues 
which  are  used  to  compute  the  modal  model  (see  Equation  (17))  and  the  modal  vectors  which 
make  up  the  transformation  matrix. 

Recall  that  a  convenient  form  of  the  model  is  for  unity  modal  mass  or  unity  [V]  (see  Equation 
(52))  scaling.  This  results  in  equations  of  the  form  of  those  in  Equations  (18)  through  (20). 
Examination  of  the  modal  stiffness  and  damping  matrix  reveals  that  the  estimate  of  the  damping 
ratio  f  is  involved  in  both  matrices.  Unfortunately,  damping  is  a  difficult  parameter  to  estimate. 
This  is  one  of  the  major  limitations  of  the  accuracy  of  an  experimental  modal  model.  Fortunately,  if 
great  care  is  taken  in  the  measurements,  the  magnitudes  of  this  error  are  not  great  enough  to  cause 
more  variation  than  found  in  normal  experimental  error. 

Another  source  of  error  is  truncation.  Truncation  errors  occur  in  two  forms:  geometry  and  modal. 
Geometry  truncation  is  a  problem  that  occurs  when  not  enough  physical  coordinates  are  defined 
to  adequately  describe  the  dynamics  of  the  structure.  Higher  order  mode  shape  patterns  are  not 
properly  defined  unless  enough  points  are  defined  along  the  shape  to  describe  it.  A  good  rule  of 
thumb  is  to  apply  Shannon’s  sampling  theorem  to  the  highest  order  mode  expected.  Geometry 
truncation  also  occurs  when  all  pertinent  translational  and  rotational  degrees  of  freedom  are  not 
measured.  If  a  structure  exhibits  motion  in  all  translational  degrees  of  freedom  and  only  one  is 
measured,  the  associated  error  is  defined  as  geometry  truncation.  In  general,  the  number  of  data 
points  should  be  much  greater  than  the  number  of  modes  of  interest  to  avoid  geometry  truncation. 

Modal  truncation  refers  to  the  number  of  modes  included  in  the  data  set.  From  Equation  (38),  the 
modified  mode  shapes  are  a  linear  combination  of  the  original  mode  shapes.  Therefore,  the  rank 
of  the  original  modal  matrix  limits  the  possible  dynamic  changes  that  can  be  calculated.  The  lower 
limit  of  the  number  of  modes  required  for  even  simple  structures  is  six  l3 z\  to  have  sufficient  rank  to 
accurately  predict  the  results  for  the  first  few  modes.  A  good  rule  of  thumb  is  to  include  several 
modes  beyond  the  frequency  range  of  interest  to  insure  the  validity  of  the  results  within  the 
frequency  range  of  interest.  Another  serious  modal  truncation  error  occurs  when  the  rigid  body 
modes  i  .c  a  free-free  structure  are  not  included  when  that  structure  is  tied  to  ground. 

Based  on  the  preceding  discussion  it  is  apparent  that  the  use  of  modal  modeling  programs  with 
experimental  data  requires  carefully  acquired  data  and  good  parameter  estimation  results.  These 
problems  can  be  overcome  by  carefully  designing  the  modal  test  and  using  the  proper  parameter 
estimation  algorithms  for  the  given  data  l33'.  This  technique  works  equally  well  with  analytical  data 
and  has  been  implemented  in  this  manner  by  Structural  Measurement  Systems  ^34l 

The  issue  of  complex  versus  real  modes  has  been  debated  greatly  in  recent  years.  To  be 
completely  accurate  the  complex  form  of  the  modal  modeling  technique  should  be  used  when 
nonproportional  or  heavy  damping  exists  in  a  structure.  Using  a  real  normal  mode  in  this  case  will 
cause  erroneous  results  One  compromise  is  to  use  a  normalized  set  of  real  modes  derived  from 
the  measured  complex  modes  described  in  Section  2.7. 


The  use  of  beam  modifications  greatly  increases  the  capability  of  modal  modeling.  Simple 
scalar  modifications  and  lumped  masses  are  limiting  and  unrealistic.  Beam  modifications  require 
rotational  information  at  the  modification  points.  This  information  is  not  readily  available  but  can 
be  obtained  with  some  effort  experimentally  or  analytically.  Once  rotational  information  is 
readily  available  from  experimental  sources,  modal  modification  will  become  a  more  powerful  trouble 
shooting  tool. 


3.  SENSITIVITY  ANALYSIS 


3.1  Sensitivity  Analysis  Overview 

Sensitivity  Analysis  is  an  approximate  technique  that  determines  the  rate  of  change  of 
eigenvalues  and  eigenvectors  using  a  Taylor  expansion  of  the  derivatives.  This  technique  was 
developed  by  Fox  and  Kappor  t3”  and  Garg  ^  ^  initially  in  the  late  Sixties  and  early  Seventies. 
Van  Belle  and  VanHonacker  I20,37!  further  developed  its  use  with  mechanical  structures  and 
implemented  it  for  use  directly  on  modal  parameters.  This  technique  is  approximate  because  only 
one  term  (differential)  or  two  terms  (difference)  of  the  series  expansion  are  used  to  approximate 
the  derivative. 

Sensitivity  Analysis  is  useful  in  two  ways.  First,  if  a  certain  type  of  modification  of  a  structure  is 
required.  Sensitivity  Analysis  determines  the  best  location  to  make  effective  structural  changes. 
Sensitivity  Analysis  also  is  used  to  predict  the  amount  of  change  by  linearly  interpolating  the  amount 
of  change  from  the  sensitivity  to  achieve  the  desired  dynamic  behavior.  This  last  method  is  very 
time  consuming,  especially  when  using  difference  sensitivities  to  maintain  accuracy. 

3.2  Formulation  of  Sensitivity  Analysis 

Van  Belle  initially  developed  the  expression  for  differential  and  difference  sensitivities  using  the 
theory  of  adjoint  structures.  A  design  problem  is  reduced  to  a  classical  optimization  problem  and 
Tellogen’s  Theorem  is  applied.  Wang  and  De  Landsheer  ^  independently  derived  the  sensitivity  of 
the  dynamic  transfer  function  by  direct  algebraic  manipulation  and  partial  differentiation.  This 
results  in  a  general  expression  for  the  derivative  of  the  transfer  function  matrix.  Starting  with  the 
basic  second  order  differential  equation  of  a  mechanical  structure: 

[MKi>[CKiM*K*}  ={/(')}  (116) 

where  M,  C  and  K  are  the  structure’s  mass,  damping  and  stiffness  matrices  of  order  n  by  n.  Laplace 
transforms  are  used  to  arrive  at  a  set  of  linear  equations. 

[J3(s)K*}={F(s)}  (117) 


where, 

[fl(.r)]  =  [[A/jiMCM*]] 
or, 


[H][F}={x}  (118) 

where, 

[W]=  [[MlsMCJH*]]'1 

From  Equations  (117)  and  (118),  an  identity  equation  can  be  formulated  as: 

[B  ][//]  =  [#/][ll]  =  [/J  (119) 


For  a  modified  system  with  an  added  mass,  damper  or  stiffener,  a  similar  identity  equation  can  be 
derived  : 


[BL[WU  =  ["MBJ™  =  [/] 


(120) 


[WU  =  [W]  MAW] 

[BU  =  [fl]  +  [AS] 

where  [A//]  and  [AS]  are  formed  by  applying  Equation  (117)  and  (118)  with  added  mass(es), 
damper(s)  or  stiffener(s).  Substitution  is  made  to  Equation  (120)  yielding : 

[tf  ]»  [Bln  =  [H]  [B]  +  [W]  [Afl]  ♦  [A//]  [J 5]  ♦  [AW] [AB]  =  [/]  (121) 

By  the  identity  Equation  (119),  Equation  (121)  can  be  simplified  as: 

[AW]  [[B]  +  [AS]]  =  -[//]  [AS] 

By  postmultipling  the  equation  above  by  [H]  and  then  taking  the  inverse  of  ([/]  +  [AB]  [W]),  an 
expression  for  [AW]  can  be  obtained  as: 

[AW]  =  -  [W  ]  [AS  ]  [W  ]  [[/  ]  ♦  [Afl  ]  [W  j'1  (122) 

Recall  the  Macalaurin  series  and  Taylor  series  as : 

(l+JC^-l+Jt-y-...  (123) 

A/  =  /-/(xo)=/'(^o)Ar  +  ^y^-(Ax)2  (124) 

Applying  Equation  (123),  Equation  (122)  can  be  expanded  as  a  matrix  polynomial  form: 


[AW]  =  -  [W][AB][W]  +  [W][Afl][W][AB][W]  ♦ . . .  (125) 

The  first  and  second  order  derivative  of  the  transfer  function  matrix  [H]  can  then  be  obtained 
according  to  Equation  (124): 

[W]'  =^='[«]Iff1[W]  (126) 

dp  bp 


[Wf^^W^IWl^V]  027) 

dp2  A/J  Ap 

The  transfer  function  matrix  [H(s)J  can  be  conveniently  expressed  as  a  summation  of  modal  vectors 
in  partial  fraction  form  by  expanding  Equation  (119)  as  : 

,"<,)1\?1  <128) 


where. 


=  K  •  { 0}r +n  -  { t/’Jr  . 


r  =  1,  2,  3 . N 


Differentiating  Equation  (128)  with  a  general  variable  p  (could  be  mass,  damping  or  stiffness),  the 
derivative  of  the  transfer  function  matrix  can  be  found  as : 

IK 2  (TSr  *  uM 

With  the  substitution  of  Equation  (128)  into  Equation  (126),  the  derivative  of  the  transfer  function 
matrix  can  be  clearly  shown  as  a  summation  of  power  series  of  - — - — -  through  the  use  of  partial 
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fraction  expansion : 
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By  factoring  out  the  same  terms,  the  derivative  of  the  transfer  function  matrix  can  be  simplified  as 

i"0"i\iJ<7V,*(^x>H  ° 


where. 
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A  comparison  between  Equations  (129)  and  (130)  related  to  Equation  (126)  yields : 


A  =  -{'l>}r{'l>}rT  ~ 
dp 


diWrii’V) 


Expanding  Equation  (131),  the  eigenvalue  sensitivity  can  be  obtained  through  the  following 
procedure. 


{  ^}r  {  i>)rT  -r~  =  -  {  ^}r  {  V’}/ 

dp 


where  a„  is  a  scalar  and. 
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Hence,  the  eigenvalue  sensitivity  is  : 


=  -a„  =  -{4,}rT 


Also,  expanding  Equation  (132),  the  eigenvector  sensitivity  can  be  acquired. 
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where, 


bf  =  {^1  AMI 


ds  dp 
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Investigating  the  equation  above,  a  further  detailed  expression  of  eigenvector  sensitivity  about  each 
modal  coordinate  can  be  found  as  : 
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A  comparison  between  the  next  two  equations  yields  the  eigenvector  sensitivity  expression, 
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(135) 


Following  the  same  derivation  process,  starting  from  Equation  (127),  the  second  order  eigenvalue 
sensitivity  equation  is : 
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3.3  Analysis 


The  differential  sensitivity  analysis  technique  utilizes  the  term  on  right  hand  side  of  Equation  (126). 
Different  expressions  are  developed  for  the  mass,  stiffness  and  damping  cases.  For  the  calculation  of 
the  sensitivity  of  an  eigenvalue,  A*,  only  the  corresponding  eigenvector  is  necessary.  Calculation  of 


eigenvalue  derivations  do  not  require  complete  information  on  the  dynamics  of  the  structure  ^38l 

Finite  difference  or  difference  sensitivities  use  the  terms  on  the  right  hand  side  of  Equations  (126) 
and  (127)  trying  to  improve  the  approximation  of  the  differential  sensitivity.  Nevertheless,  a  better 
approximation  is  obtained  only  when  the  change  of  the  structure  parameter  is  small.  If  the 
magnitude  of  change  is  increased  beyond  a  certain  value,  the  result  will  be  even  worse.  Equation 

(136)  shows  that  a  term  - — - — -  is  involved  in  the  second  order  derivative  of  an  eigenvalue.  If 

C\fc-\n) 

there  are  two  close  modes,  the  term  will  become  large  so  the  second  order  derivative  contribution 
dominates  the  approximation.  Further  more,  if  two  adjacent  modes  are  very  close  to  each  other,  the 
term  will  diverge  so  the  result  will  be  unacceptable.  Hence,  care  must  be  taken  when  a  set  of  modal 
data  shows  repeated  eigenvalues. 

From  this  discussion,  it  is  seen  that  the  expressions  used  to  compute  differential  or  difference 
sensitivity  from  modal  parameters  are  in  the  form  of  the  transfer  function  matrix.  Because  only 
one  or  two  terms  are  used  from  the  Taylor  expansion,  the  technique  is  an  approximate  one.  Since 
only  modal  parameters  are  necessary,  this  technique  is  equally  applicable  to  experimental  or 
analytical  data  Currently  it  is  implemented  with  experimental  data  t38l 

VanHonacker  ^38'  has  shown  this  method  to  be  accurate  for  only  small  incremental  changes.  The 
differential  method  is  far  less  accurate  than  the  difference  method.  For  small  changes  of  mass, 
stiffness,  or  damping  the  differential  technique  will  accurately  predict  the  eigenvalue  shift.  For 
more  significant  parameter  changes,  the  difference  technique  is  recommended.  Therefore, 
sensitivity  analysis  has  only  limited  application  in  the  prediction  of  the  effects  of  structural 
modifications. 

Sensitivity  is  extremely  useful  as  a  preprocessor  to  Modal  Modeling  or  Finite  Element  Analysis 
techniques.  The  sensitivities  of  a  structure  can  be  computed  rapidly  from  the  modal  parameters 
to  determine  the  optimal  location  at  which  to  investigate  a  modification.  Furthermore,  the 
sensitivity  value  is  useful  in  determining  how  much  of  a  modification  is  required.  Therefore, 
Sensitivity  Analysis  is  a  valuable  tool  in  the  optimization  of  a  design. 

This  technique  has  several  limitations.  First,  the  results  are  only  as  good  as  the  modal  parameters 
used  in  the  calculations.  Therefore,  all  of  the  experimental  errors  and  parameter  estimation 
limitations  which  hinder  other  modeling  methods  apply  to  sensitivity  techniques  as  well.  Notably,  a 
limited  number  of  modes  are  available  from  zero  to  the  maximum  frequency  measured.  Although 
not  currently  implemented  with  analytical  data,  any  inaccuracies  in  an  analytical  model  would 
similarly  deteriorate  the  calculations  when  used  with  modal  data.  In  the  experimental  case,  geometry 
truncation  errors  are  significant  due  to  the  exclusion  of  rotational  degrees  of  freedom. 

As  a  preprocessor  to  other  modeling  techniques,  sensitivity  has  advantages.  The  computations  are 
fast  and  stable,  especially  when  compared  with  a  complete  eigensolution.  It  is  intuitive  in  nature 
because  it  provides  rates  of  change  that  allow  the  selection  of  the  best  type  and  location  of 
modification  as  well  as  a  comparison  of  different  modifications.  This  provides  a  large  amount 
of  information  that  offers  much  insight  into  the  dynamic  behavior  of  a  structure. 


4.  IMPEDANCE  MODELING 


4.1  Impedance  Modeling  Overview 

The  general  impedance  method  was  first  introduced  by  Klosterman  and  Lemon  (391  in  1969. 
Due  to  the  state  of  measurement  equipment  at  that  time  the  method  was  not  pursued  further. 
As  the  ability  of  Fourier  analyzers  to  accurately  measure  frequency  response  functions  improved 
in  the  late  Seventies,  the  interest  in  General  Impedance  Techniques  was  renewed.  Two 
techniques  are  developed  in  this  chapter  using  experimental  data.  Both  methods  use  measured 
frequency  response  functions  or  synthesized  frequency  response  functions. 

The  general  impedance  technique  is  formulated  in  two  ways:  frequency  response  and  dynamic 
stiffness.  The  dynamic  stiffness  approach  was  initially  developed  by  Klosterman  ^  and 
implemented  by  Structural  Dynamics  Research  Corporation  as  SABBA  (Structural  Analysis  using 
the  Building  Block  Approach)  ^41l  This  technique  is  primarily  used  to  couple  together 
structures  to  predict  the  total  dynamic  characteristics  using  the  concept  of  superposition.  Thus,  the 
phrase  building  block  approach  was  applied  to  this  technique. 

The  frequency  response  method  was  published  and  implemented  by  Crowley  and  Klosterman  ^  at 
the  Structural  Dynamics  Research  Corporation  in  1984  and  referred  to  as  SMURF  (Structural 
Modifications  Using  Response  Functions).  It  is  primarily  a  trouble  shooting  technique.  It 
operates  on  frequency  response  functions;  therefore,  no  modal  model  is  necessary  to  investigate 
structural  modifications.  Problems  may  be  solved  by  acquiring  several  frequency  response 
functions  and  investigating  the  effects  of  modifications.  The  modified  frequency  response 

function  is  computed  as  a  function  of  frequency  by  simple  block  operation  using  the  original 
frequency  response  function  and  a  frequency  representation  of  the  structural  modification. 

Consider  once  again  the  multiple  degree  of  freedom  linear  second  order  differential  equation: 

[MKi'MCKiM*K*}={/}  (137) 

which  describes  the  equations  of  motion  of  a  structure.  Laplace  transforms  or  assumption  of  the 
solution  {x}  =  { X}e *  yields: 

[[M]s2+[C]s+[K]]  =  F(s)  (138) 

where  the  initial  conditions  are  assumed  to  be  zero.  The  substitution: 

[B(s)]  =  [[M]rMCM*]] 
is  made  to  Equation  (138)  to  yield: 

(«(s)l{*(s)}={F(r)}  (139) 

where  (B(s)]  is  the  system  matrix.  The  system  transfer  function  [H(s)J  is  defined  as  [^(s)]"1  allowing 
Equation  (140)  to  be  written  as: 


I 


‘"(s)1'lw  (140) 

In  practice,  the  transfer  function  matrix  is  not  measured.  Experimentally,  the  Fast  Fourier 
Transform  measures  only  the  frequency  response  function  at  s=/w  in  the  complex  plane  yielding: 

["Ml-  7m  (14,) 

This  expression  is  expanded  using  partial  fractions  and  including  inertial  restraint  and  residual 
flexibility  for  mass  and  stiffness  as: 

[//(u,)]  =  [/?/]+  £  +  77 - V-  ♦  [  Rp]  (142) 
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Equation  (142)  is  the  expression  that  forms  the  basis  of  impedance  methods.  Measured 
frequency  response  functions  are  represented  analytically  by  this  equation.  If  the  modal 
parameters  exist  from  either  Finite  Element  Analysis  or  an  experimental  modal  model,  the 
frequency  response  function  can  be  synthesized  from  the  model  by  substituting  in  the 
eigenvalue  and  eigenvector  information.  The  use  of  residual  terms  increases  the  validity  of  the 
analytical  expression. 

4.2  Dynamic  Stiffness  Method 

Historically,  the  dynamic  stiffness  method  was  implemented  first  by  Klosterman  and  Lemon  as  a 
building  block  approach.  This  approach  is  used  to  predict  the  response  of  two  components 
coupled  together  or  the  effects  of  a  structural  change,  when  the  dynamic  characteristics  are  available 
in  dynamic  stiffness  form. 

{F}  =  [//m)  (143) 
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[H  I"1  is  the  inverse  of  the  frequency  response  function  matrix.  The  dynamic  stiffness  matrix  may 
consist  of  experimental,  analytical,  or  general  matrix  components  in  their  dynamic  stiffness 
form. 

A  system  with  two  points  p  and  q  has  inputs  at  both  points.  Frequency  response  functions  are 
linear  for  linear,  time  invariant  systems.  Therefore,  the  principle  of  superposition  is  used  to  write 
the  steady  state  response. 

{Fp}  =  [MppV{Xp}*[Hm]-1{X,}  (144) 

{F,}  =  [W J-1  {*„}  ♦  [HJ-1  {X,} 


Now  consider  a  technique  to  couple  two  systems,  1  and  2,  at  points  p  and  s.  The  equations  for 
sinusoidal  input  on  the  second  structure  at  points  s  and  t  are: 

{Z7.}  =  [//-I'1  {*«}  ♦  l H A1  {*<}  045) 


{Ft}  =  [/U1  W  ♦  W*Yl  W 


Now,  equilibrium  and  compatibility  conditions  are  enforced  upon  the  system 


FP*F,  =  0 


(146) 


XV=X, 


to  combine  the  two  sets  of  equations  into  a  system  matrix. 
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(147) 


This  set  of  equations  represents  a  set  of  complex  linear  simultaneous  equations  of  the  form 
F=KX. 

This  method  can  be  used  with  several  sources  of  frequency  response  function  data.  The  first  is 
direct  computation  of  mass  and  stiffness  matrices  to  satisfy  Equation  (143).  The  second  is  to  utilize 
the  impedance  properties  of  many  components.  For  example,  mass,  stiffness,  beam,  or  bearings 
are  modeled  using  a  library  of  impedance  representations  in  SDRC’s  SABBA  modeling 
program.  The  third  is  the  direct  measurement  of  frequency  response  functions  and  inverting  them 
to  assemble  the  dynamic  stiffness  matrix.  Last  is  the  computation  of  frequency  response  function 
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data  from  the  Equation  (141)  from  a  modal  model. 

The  last  two  methods  require  an  inversion  of  the  matrix  [H].  Care  must  be  taken  to  prevent  this 
inversion  from  becoming  ill-conditioned.  The  use  of  frequency  response  function  data  was  found  to 
be  undesirable  by  Klosterman  initially,  due  to  instability  in  the  numerical  calculations  caused  by 
poor  measurement  quality.  Hence,  the  general  impedance  technique  is  implemented  using 


synthesized  frequency  response  functions  from  modal  data.  The  possible  errors  associated  with 
using  a  modal  model  are  discussed  in  Section  2.8  .  Considering  the  possible  errors  associated 
with  a  modal  model,  this  technique  is  still  preferable  to  direct  frequency  response  function 
calculations. 

The  modal  method  is  implemented  in  two  fashions,  the  unconstrained  modal  method  and  the 
constrained  modal  method.  Both  methods  are  compatible  with  modal  models  or  Finite  Element 
Analysis.  These  methods  are  used  when  a  structure  is  free  or  constrained  respectively  during  test 
or  analysis. 

Consider  the  degrees  of  freedom  of  the  connection  point  to  be  designated  by  qm  in  modal 
coordinates,  and  [4>]  is  the  eigenvector  abbreviated  to  include  only  the  coordinates  at  the 
connection  points.  The  modal  model  of  a  component  is  expressed  as: 


% -J*  fA/j  ■  [qm]  =  {Fm}  ♦  [$f  {F»} 


where:  {Fm}  -  forces  due  to  external  connection  of 
another  component 

{Fn}  -  additional  external  forces 

This  set  of  equations  could  represent  an  experimental  modal  component. 
Now  consider  a  finite  element  component  in  standard  form 


K*  K*  ?  M, 

Kn,  Krn  M„ 


r-  f  Fj 

mi  j^mj  |*Fm  +  F\ 


where:  {F2}  is  external  forces  on  the  second  component 
{Ai}  is  remaining  coordinates 
{^m}  is  the  coordinates  of  the  connection  points 


The  two  systems,  modal  and  analytical,  are  combined  by  enforcing  the  compatibility  and  equilibrium 
conditions  yielding: 


(150) 
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This  matrix  can  be  condensed  into  the  general  form  of: 


[[K]-u?[Af]]{*}={F} 


051) 


by  substitution. 

This  is  the  standard  form  used  by  finite  element  programs.  Recall  that  the  dynamic  matrix  is 
formed  from  frequency  response  function  information  and  must  be  solved  frequency  by  frequency 
using  a  determinant  search  algorithm  to  determine  the  new  model  properties.  In  its  initial 
implementation  the  series  of  equations  were  used  to  solve  for  the  modified  frequency  response  due 
to  a  structural  modification  or  coupling  with  another  structure.  The  unconstrained  modal 
method  requires  a  large  number  of  degrees  of  ireedom  depending  upon  the  severity  of  the 
constraints  placed  upon  the  component.  This  led  to  the  development  of  the  constrained  modal 
method. 

The  constrained  modal  method  is  used  when  the  test  or  analysis  is  conducted  with  the  structure 
constrained  in  some  manner.  This  is  desirable  in  several  situations.  Rigid  body  modes  are 
eliminated  which  is  desirable  in  light  of  the  fact  that  they  are  difficult  to  determine  in 
experimental  tests.  A  constrained  mode  model  closely  represents  the  final  structure.  Klosterman 
has  shown  f40*  that  the  damping  characteristics  are  frequency  dependent;  therefore,  a  minimal 
shift  in  the  dominant  modes  will  result  in  a  more  accurate  representation  of  damping.  Fewer 
modes  or  degrees  of  freedom  are  necessary  if  the  structure  is  analyzed  or  tested  with  constraint 
similar  to  the  final  composite  structure. 

Klosterman  has  proposed  a  method  to  derive  a  representation  when  the  connection  coordinates 
{A’m}  =  0.  This  result  is: 


[*F  Jf{<5 
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where:  m  =  fixed  connection  coordinates 
i  =  remaining  coordinates 
{</,}  =  modal  coordinate  of  the  remaining  points 

[R  ]  =  -[#f  fKj  =  forces  applied  at  the  connections  to  restrain  each  normal  mode 

/w  =  [*™»j  -  [*r 

When  the  coordinates  { qm }  are  not  redundant,  [ K ]  =  0  and  Equation  (152)  is  the  best  form  to 
use.  All  of  the  information  necessary  to  use  with  Equation  (152)  is  available  from  any  modal  analysis 
if  all  of  the  forces  of  constraint  are  available.  This  point  is  the  largest  disadvantage  to  this 


technique.  A  more  detailed  development  is  found  in  Reference  4. 

The  dynamic  stiffness  approach  has  evolved  from  the  level  where  single  frequency  response 
functions  were  computed  for  modified  structures,  to  where  complex  beam  and  matrix  constraints  can 
be  placed  on  a  structure  and  the  resulting  modal  parameters  are  found  using  a  frequency 
dependent  determinant  search  method.  This  technique  has  been  implemented  on  super  mini¬ 
computers  ^  and  is  available  for  on-test-site  trouble  shooting  and  design  cycle  applications  .  It 
relies  on  modal  data;  therefore,  any  limitations  noted  on  modal  parameter  identification  in  Section 
2.8  are  relevant  to  this  technique. 

4.3  Frequency  Response  Method 

The  frequency  response  method  is  a  pure  trouble  shooting  technique  that  operates  on  measured 
data  or  synthesized  data  It  is  implemented  ^  to  be  used  for  on-site  structural  modifications  or 
verification  of  a  modal  model.  Only  single  constraint  structural  modifications  are  practical. 
This  technique  is  an  extension  of  the  technique  originally  proposed  by  Klosterman.  It  became 
feasible  because  the  numerical  problems  associated  with  using  measured  data  were  reduced  as 
measurement  techniques  improved  in  recent  years. 

Frequency  response  testing  results  in  equations  of  the  general  form  : 


{*}  =  [//]{F} 
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In  the  following  section,  equations  will  be  developed  using  simple  block  manipulations  of 
frequency  response  functions  to  couple  two  structures,  add  mass  to  a  structure,  add  a  stiffener  to 
a  structure,  and  constrain  a  point  on  a  structure  to  ground. 

First,  the  coupling  of  two  structures  is  treated.  The  forcing  vector  is  split  into  external  {Fe}  and 
internal  {F,}  forces,  yielding: 


{*}  =  [//]{Fe}  +  [tf]{F,} 


(154) 


which  is  manipulated  into: 

[!/M//j{j£|  =  l"]{Fe}  (155) 

The  uncoupled  equations  of  two  structures,  1  and  2,  become: 
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Now  the  equilibrium  and  compatibility  equations  are  applied  at  the  connections  to  compute  the 
frequency  response  of  the  modified  structure.  The  bar  denotes  the  combined  structures  1  and 


2,  p  refers  to  the  response  point  on  1,  q  refers  to  the  input  point  on  1,  and  k  is  the  connection  point . 
The  equation  to  compute  the  new  frequency  response  is: 


#pgl  ffitgi 
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(157) 


Now  consider  structure  1  with  point  m  tied  to  ground  such  that  Xm  *  0  The  new  frequency 
response  between  p  and  q  is  found  by  writing  the  response  of  point  m  due  to  input  at  m  and  q  as: 


Xm=HmmFm  +  Hm 


since  Xm  -  0 


Fm  =  Ft  (159) 

The  response  of  point  p  is  also  written  in  terms  of  the  reaction  force  at  m  and  input  at  q. 

Xp  =  Hpm  Fm  *  Ff  (160) 

Equation  (158)  is  substituted  into  (159)  to  obtain  the  constrained  response  at  point  p. 

Xp  =  (H„-HpnH,LHmi)Ft  (161) 

By  simply  dividing  by  Ft  the  frequency  response  function  is  found. 

/C»  (162) 

H „  is  computed  frequency  by  frequency  from  the  measured  frequency  response  functions.  Note 

that  if  a  structure  is  to  be  constrained  at  a  point  m,  a  driving  point  measurement  at  m  and  cross 
measurements  between  m,  p,  and  q  are  necessary  to  compute  the  modified  frequency  response 
function. 

Fixing  a  point  to  ground  is  a  special  case  of  the  more  general  one  of  constraining  a  single  point 
of  a  structure.  Crowley  has  shown  ^  that  the  constraint  is  modeled  by  finding  a  stationary 
point  of  the  constrained  quadratic  functional  by  minimizing  the  strain  energy: 

j{X}aml{X}-{X}B{F}  (163) 


subject  to  the  constraint  f/  j7”  {X}  =  0 


where:  J  is  a  matrix  of  constraint  coefficients 


Employing  Lagrange  multipliers  {/?}  is  equivalent  to  solving  a  linear  equation  system. 


Algebraic  manipulation  results  in  the  general  form  for  constrained  response: 

[H]  -  [[/]  -  [H]  [J]  [[Jf  [H]  [J  j'1  [/]*]  [//]  (165) 

Equation  (165)  is  used  to  develop  general  equations  for  addition  of  mass  at  a  point,  and  a  stiffener 
addition  between  two  points.  The  cases  of  constraining  one  point  to  ground  and  coupling  two 
structures  are  already  developed. 

To  compute  the  effect  of  addition  of  mass  at  point  q  to  Equation  (165)  takes  the  form: 


(H„)  [/*„  +  W-n]-1  (//,*) 


(166) 


where:  H ^  is  the  analytical  driving  point  of  the  mass  addition  at  point  q. 

The  equation  for  the  addition  of  a  stiffener  between  q  and  m  for  the  modified  response  at  in 
general  form  is: 


(167) 


The  expression  for  adding  stiffness  at  points  q  and  m  is  computed  by  adding  the  analytical 
stiffness  to  H„  and  Hm 

The  frequency  response  method  is  a  simple  technique  to  utilize  in  trouble  shooting  situations. 
Measurements  are  required  at  and  between  the  driving  point  and  ail  points  of  constraint  as  well 
as  the  desired  response  point.  Generally,  driving  point  information  is  needed  at  the  driving  point 
and  points  of  constraint.  Cross  measurements  are  also  needed  between  each  of  these  points. 
Crowley  suggests  only  one  constraint  be  applied  at  a  time  to  avoid  matrix  inversions  l42’ .  Based  on 
the  current  implementation,  it  appears  that  an  impact  testing  technique  is  best  suited  to  quickly 
acquire  the  necessary  frequency  response  functions  to  investigate  modifications  using  the  frequency 
response  method. 


4.4  Analysis 

The  general  impedance  technique  is  a  method  that  employs  the  use  of  frequency  response  data  to 
investigate  coupling  of  structures  and  structural  modifications.  The  dynamic  stiffness  approach  is 
the  more  powerful  of  the  two.  The  advantage  of  this  technique  is  that  a  large  array  of 
investigations  may  be  conducted.  Also,  the  method  is  not  as  cumbersome  as  the  component 
modal  synthesis  technique  (Section  5)  because  it  has  a  reduced  number  of  degrees  of  freedom. 
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Due  to  the  dynamic  stiffness  formulation,  modal  models  of  experimental  or  analytical  origin  may 
be  combined  or  modified  using  an  array  of  mass,  stiffness,  damping,  beam,  or  matrix  elements 
represented  in  impedance  form.  This  brings  more  analytical  capability  directly  to  an  engineer  in 
a  mini-computer  environment  as  implemented  currently  l41)  when  compared  with  modal  modeling 
techniques. 

Due  to  numerical  problems,  Klosterman  recommends  use  of  synthesized  frequency  response 
functions  to  build  the  dynamic  stiffness  matrix.  This  introduces  errors  made  in  modal  parameter 
estimation,  but  reduces  numerical  problems  associated  with  noisy  frequency  response  functions 
measurements  because  the  parameter  estimation  process  fits  a  smooth  curve  through  the  measured 
data.  Modeling  of  this  type  requires  carefully  acquired,  and  properly  calibrated,  data  to  obtain 
the  best  modal  model  possible. 

One  of  the  major  problems  with  the  dynamic  stiffness  approach  is  the  determination  of  the 
[//  J-1  matrices.  The  inverse  of  the  matrix  must  exist.  This  problem  forces  the  use  of  modal  data 
because  the  inversion  process  is  numerically  unstable  ^43l  The  number  of  modes  must  be  much 
greater  than  the  number  of  constraint  points  to  insure  the  existence  of  the  inverse. 

The  most  serious  limitation  of  the  dynamic  stiffness  approach  is  the  computational  speed  and 
stability.  When  implemented  initially,  only  individual  frequency  response  functions  were  computed 
for  the  resultant  structure.  This  implementation  is  efficient,  but  the  computations  are  somewhat 
unstable  unless  synthesized  frequency  response  functions  are  used.  This  led  to  the  application  of 
a  determinant  search  algorithm  to  compute  the  resultant  eigenvalues  and  eigenvectors.  This 
algorithm  is  not  computationally  efficient  because  the  equations  are  solved  frequency  by 
frequency.  This  fact  has  led  to  more  widespread  use  of  component  modal  synthesis  techniques  that 
are  described  in  Chapter  5. 

The  frequency  response  method  is  implemented  for  single  constraint  situations  to  avoid  the 
matrix  inversion  problem.  This  makes  it  useful  for  trouble  shooting  situations.  Since  measured 
frequency  response  functions  are  used  in  the  calculations,  no  modal  model  is  necessary.  Therefore, 
if  the  necessary  frequency  response  functions  exist,  modifications  may  be  made  directly  to  obtain  the 
modified  frequency  response.  This  makes  the  frequency  response  technique  the  fastest  trouble 
shooting  technique,  but  only  modified  frequency  response  functions  are  computed  not  modal  data. 
This  technique  is  slower  than  modal  modeling  when  modal  data  is  desired. 

Data  is  required  for  driving  points  and  cross  measurements  at  the  constraints,  response  point,  and 
driving  point  to  compute  the  modified  frequency  response.  Therefore,  the  impact  testing  or 
multiple  reference  techniques  are  most  convenient  for  the  frequency  response  method. 
Synthesized  frequency  response  functions  can  be  used  in  this  technique  when  the  desired  frequency 
responses  are  not  available,  but  the  advantage  of  avoiding  modal  parameter  estimation  is  lost. 
Modal  modeling  is  a  better  alternative  at  this  point  because  the  entire  set  of  modified  modal 
parameters  is  computed.  Only  individual  modified  frequency  response  functions  result  from  the 
frequency  response  method. 

The  frequency  response  technique  is  sensitive  to  measurement  errors  such  as  leakage,  aliasing,  and 
random  noise.  Wang  l44'  has  found  the  errors  largest  at  anti-resonance.  This  is  due  to  the  fact  that 
the  signal  to  noise  ratio  is  poor  at  anti-resonance.  Therefore,  as  the  magnitude  of  modification  or 
constraint  increases,  the  more  the  accuracy  of  the  calculations  will  deteriorate  because  frequencies 
will  shift  closer  to  or  past  anti-resonances. 


5.  COMPONENT  DYNAMIC  SYNTHESIS 


5.1  Introduction 

Component  dynamic  synthesis  is  an  analytical  procedure  for  modeling  dynamical  behavior  of  complex 
structures  in  terms  of  the  properties  of  its  components  or  substructure.  The  procedure  involves 
explicitly  every  component  in  the  structure  the  advantages  of  which  are  many-fold:  analysis  and 
design  of  different  components  of  the  structure  can  proceed  independently,  component  properties 
can  be  obtained  from  tests  and/or  analysis,  and  the  size  of  the  built-up  structure  analysis  problem  can 
be  reduced  to  manageable  proportions. 

Component  synthesis  with  static  condensation  t45,46l  has  long  been  used  for  improving  efficiency  of 
static  analysis.  In  this  method,  known  as  substructuring  technique,  unique  or  functionally  distinct 
components  of  a  structural  system  are  analyzed  separately,  condensed,  and  then  combined  to  form  a 
reduced  model.  This  reduced  model,  having  fewer  degrees  of  freedom,  is  generally  more  economical 
to  analyze  than  the  original  structural  model.  The  static  condensation  is  an  exact  reduction 
procedure. 

In  dynamic  analysis,  exact  reduction  of  an  individual  component  is  dependent  upon  the  natural 
frequencies  of  the  total  structural  system  which  are  yet  unknown  at  the  component  level.  Frequency 
independent  or  iterative  reduction  methods  must  therefore  be  used,  which  introduce  approximations. 
The  various  reduction  methods  are  collectively  known  as  component  dynamic  synthesis  or  modal 
synthesis  (CMS). 

The  objectives  of  this  section  are  to  review  the  state  of  the  art  in  component  dynamic  synthesis  and  to 
develop  and  implement  an  improved  dynamics  synthesis  procedure. 

5.1.1  Dynamic  Synthesis  Methods 

In  order  to  review  the  existing  dynamic  synthesis  procedures,  it  is  necessary  to  define  certain 
frequently-used  terms.  A  component  or  substructure  is  one  which  is  connected  to  one  or  more 
adjacent  components  by  redundant  interfaces.  Discrete  points  on  the  connection  interface  are  called 
boundary  points  and  the  remainder  are  called  interior  points.  The  following  classes  of  modes  are 
commonly  used  as  basis  components  in  component  dynamic  model  definition.  Details  of  these  mode 
sets  is  given  later  in  this  report. 

1.  Normal  Modes:  These  are  free  vibration  eigenmodes  of  an  elastic  structure  that  result  in  a 
diagonal  generalized  mass  and  stiffness  matrix.  The  normal  modes  are  qualified  as  free  or  fixed 
interface  modes,  depending  on  whether  the  connection  interfaces  are  held  free  or  fixed.  Loaded 
interface  normal  modes  simulate  intermediate  fixity  of  the  interfaces. 

2.  Constraint  Modes:  These  are  static  deflection  shapes  resulting  from  unit  displacements  imposed 
on  one  connection  degree  of  freedom  and  zero  displacements  on  all  the  remaining  degrees  of 
freedom. 

3.  Attachment  Modes:  The  attachment  modes  are  static  deflection  shapes  defined  by  imposing  a  unit 
force  on  one  connection  degree  of  freedom  while  the  remaining  connection  degrees-of-freedom 
(DOFs)  are  force  free.  If  the  structure  is  unrestrained,  this  mode  set  will  consist  of  inertia  relief 
modes.  Attachment  modes  are  also  static  modes. 

4.  Rigid-body  Modes:  These  are  displacement  shapes  corresponding  to  rigid  body  degrees  of 
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freedom.  They  may  be  considered  a  subset  of  normal  modes  corresponding  to  null  eigenvalues  or 
else  defined  directly  by  geometrical  consideration. 


5.  Admissible  Shape  Functions:  These  are  any  general  distributed  coordinates  or  space  functions, 
linear  combinations  of  which  simultaneously  approximate  the  displacement  of  all  points  of  an  elastic 
structure.  The  only  requirements  are  that  the  admissible  functions  satisfy  geometry  boundary 
conditions  of  the  component  over  which  they  are  defined,  and  satisfy  certain  differentiability 
conditions.  These  are  analogous  to  finite  element  shape  functions. 

Static  condensation  or  Guyan  Reduction  ^  is  the  simplest  of  ail  component  dynamic  synthesis 
techniques.  The  approach  is  a  direct  extension  of  static  condensation.  The  transformation  matrix  of 
static  constraint  modes  which  is  used  to  reduce  the  order  of  the  stiffness  matrix  in  static  analysis  is 
also  used  to  reduce  the  order  of  the  structure  mass  matrix.  The  kinetic  energy  of  the  interior  nodes  is 
represented  by  only  static  mode  shapes.  Drawbacks  of  this  approach  are  obvious.  The  static  modes 
are  not  the  best  Ritz  modes  for  component  dynamic  representation. 

The  concept  of  Component  Modal  Synthesis  (CMS)  was  first  proposed  by  Hurty  ^47l  Component 
members  were  represented  by  admissible  functions  (low-  order  polynomials)  to  develop  a  reduced 
order  model.  This  procedure  is  essentially  the  application  of  the  Rayleigh-Ritz  procedure  at  the 
component  level.  Hurty  extended  the  method  to  include  discrete  finite  element  models^.  This 
method  proposed  that  the  connect  DOF  of  a  component  were  fixed  or  had  a  zero  displacement. 
Hurty  then  partitioned  the  modes  of  the  structure  into  rigid  body  modes,  constraint  modes,  and 
normal  modes.  The  constraint  and  rigid  body  modes  were  found  by  applying  unit  static  load  to  each 
of  the  connection  points  individually  to  obtain  static  deformation  shapes  of  the  structure.  These 
modes  were  added  to  the  constrained  normal  modes  to  form  a  truncated  mode  set  used  in  the 
synthesis  of  the  entire  structure.  A  simplification  of  Hurty's  fixed  interface  method  was  presented  by 
Craig  and  Bampton  ^49l  Substructure  component  modes  were  divided  into  only  two  groups: 
constraint  modes  and  normal  modes.  This  resulted  in  a  procedure  which  is  conceptually  simpler, 
easier  to  implement  in  analysis  software.  Bamford  t50l  further  increased  the  accuracy  of  the  method 
by  adding  attachment  modes  which  improve  the  convergence  of  the  method.  The  attachment  modes 
are  the  displaced  configurations  of  a  component  when  a  unit  force  is  applied  to  one  boundary  degree 
of  freedom  while  all  other  boundary  DOF  remain  free  of  loads. 


Goldman  introduced  the  free  interface  method,  employing  only  rigid  body  modes  and  free-free 
normal  modes  in  substructure  dynamic  representation.  This  technique  eliminates  the  computation  of 
static  constraint  modes,  but  their  advantage  is  negated  by  the  poor  accuracy  of  the  method.  Hou 
presented  a  variation  of  Goldman’s  free-interface  method  in  which  no  distinction  is  made  between 
rigid  body  modes  and  free-free  normal  modes.  Hou’s  approach  also  includes  an  error  analysis 
procedure  to  evaluate  convergence. 
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Gladwell  developed  "branch  mode  analysis"  by  combining  free  interface  and  fixed  interface 
analysis  to  reduce  the  order  of  the  stiffness  and  mass  matrices  for  individual  substructures.  The 
reduction  procedure  depends  upon  the  topological  arrangement  of  the  substructures  in  the  model. 
Thus,  reduction  of  any  one  substructure  requires  knowledge  of  the  arrangement  of  all  substructures 
in  the  model. 

Bajan,  et  al.  developed  an  iterative  form  of  the  fixed  interface  method.  He  showed  that 
significant  improvements  in  synthesis  accuracy  can  be  achieved  by  repeating  the  reduction,  based  on 
updated  estimates  of  system  frequencies  and  mode  shapes. 

Benfield  and  Hruda  ^  introduced  inertia  and  stiffness  loading  of  component  interfaces  to  account 
for  adjacent  substructures.  The  use  of  loaded  interface  modes  is  shown  to  have  superior  convergence 
characteristics. 
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Motivated  by  the  need  to  use  experimental  test  data,  MacNeal  ^  introduced  the  use  of  hybrid 
modes  and  inertia  relief  modes  for  component  mode  synthesis.  Hybrid  modes  are  substructure 
normal  modes  computed  with  a  combination  of  fixed  and  free  boundary  conditions.  Inertia  relief 
attachment  modes  are  attachment  modes  for  components  with  rigid  body  freedoms.  MacNeal  also 
included  residual  inertia  and  flexibility  to  approximate  the  static  contribution  of  the  truncated  higher 
order  modes  of  a  component.  Rubin  ^  extended  the  residual  flexibility  approach  for  free  interface 
method  by  introducing  higher  order  corrections  to  account  for  the  truncated  modes.  Klosterman  I58' 
more  frilly  developed  the  combined  experimental  and  analytical  method  introduced  by  MacNeal. 
Hintz  discussed  the  implications  of  truncating  various  mode  sets  and  developed  guidelines  for 
retaining  accuracy  with  a  reduced  size  model 

Many  authors  have  compared  the  techniques  discussed.  No  method  clearly  appears  to  be  superior  to 
the  other.  The  constrained  interface  method  of  Craig  and  Bampton  and  Hurty  is  expected  to  be  the 
most  accurate  when  the  connect  degrees  of  freedom  have  little  motion.  The  free  interface  method 
with  the  use  of  residuals  as  proposed  by  Rubin  appears  to  be  more  accurate  than  the  constrained 
approach. 

Recent  research  has  centered  on  comparisons  of  the  various  methods.  Baker  l60!,  for  example, 
compares  the  constrained  and  free-free  approach  using  experimental  techniques  and  also  investigates 
using  mass  additive  techniques  and  measured  rotational  DOF  ^45l  This  investigation  was  motivated 
by  a  need  to  find  the  best  method  for  rigidly  connected  flexible  structures.  In  this  connection,  the 
constrained  method  produced  the  best  result.  Klosterman  ^  has  shown  the  free-free  method  to  be 
accurate  for  relatively  stiff  structures  connected  with  flexible  elements.  This  supports  Rubin’s 
conclusion  that  the  free-free  method  is  at  least  as  accurate  when  residual  effects  are  accounted 
for.  These  conclusions  are  intuitive  because  the  type  of  boundary  condition  imposed  in  the  analysis 
that  best  represents  the  boundary  of  the  assembled  structure  provides  the  best  accuracy  in  the  modal 
synthesis. 

Meirovitch  and  Hale  have  developed  a  generalized  synthesis  procedure  by  broadening  the 
definition  of  the  admissible  functions  proposed  by  Hurty  l4  l  This  technique  is  applicable  to  both 
continuous  and  discrete  structural  models.  The  geometric  compatibility  conditions  at  connection 
interfaces  are  approximately  enforced  by  the  method  of  weighted  residuals. 

The  method  due  to  Klosterman  ^581  has  been  implemented  in  an  interactive  computer  code  SYSTAN 
and  that  due  to  Herting  ^  is  available  in  NASTRAN.  The  latter  is  the  most  general  of  the 
modal  synthesis  techniques.  It  allows  retention  of  an  arbitrary  set  of  component  normal  modes, 
inertia  relief  modes,  and  all  geometric  coordinates  at  connection  boundaries.  Both  the  fixed- 
interface  method  of  Craig  and  Bampton,  and  the  MacNeal ’s  residual  flexibility  method,  are  special 
cases  of  the  general  method.  Other  analyses  presented  in  the  literature  based  on  modal  synthesis 
techniques  are  not  incorporated  into  general  structural  analysis  codes.  In  general,  there  is  a  lack  of 
sophistication  in  available  software. 

5.1.2  Damping  Synthesis  Methods 

The  methods  of  dynamic  synthesis  are  particularly  useful  and  sometimes  the  only  alternative  available 
in  damping  prediction  for  built-up  structures.  Most  frequently,  damping  has  been  synthesized  in  the 
manner  analogous  to  stiffness  and  mass  synthesis  with  the  assumption  of  proportional  damping. 
Hasselman  and  Kaplan  ^  used  complex  modes  of  components  with  nonproportional  damping. 
Obtaining  damping  matrices  in  the  presence  of  general  energy  dissipating  mechanisms  in  a  complex 
structure  is  one  of  the  complicating  factors,  however.  In  such  cases  an  average  damping  behavior  can 
be  obtained  from  tests  in  the  form  of  cyclic  energy  dissipated  versus  peak  stored  energy  correlation  or 
damping  law.  Kana  et  al.  ^  synthesized  system  damping  based  on  substructure  stored  energy  at  the 


system  modal  frequency.  Soni  I66*  developed  a  substructure  damping  synthesis  method  applicable  to 
cases  where  substructure  damping  varies  greatly  and  irregularly  from  mode  to  mode.  The  procedure 
has  been  validated  in  experimental  studies  l67l  Jezequel  employed  fixed  interface  component 
modes  together  with  mass  loaded  interface  modes  replacing  the  static  constraint  mode  in  his  damping 
synthesis  method.  Mass  loading  results  in  an  improved  representation  of  interface  flexibility  and 
dissipation;  however,  the  use  of  constrained  interface  modes  make  it  difficult  to  implement  it  in 
experimental  testing. 

The  subject  of  component  dynamic  synthesis  has  received  increasing  attention  in  recent  years. 
Reference  ^691  presents  several  detailed  reviews,  applications,  and  case  histories,  with  particular 
emphasis  on  experimental  characterization  of  component  dynamics. 

5.1.3  A  Comparison  of  the  Synthesis  Methods 

While  differing  in  their  detailed  treatment,  all  the  synthesis  methods  have  the  following  objectives: 
(1)  to  efficiently  predict  the  dynamics  of  a  structure  within  required  accuracy  for  a  minimum  number 
of  DOFs;  (2)  to  analyze  the  components  totally  independent  of  other  components  so  the  design 
process  is  uncoupled,  and  (3)  to  use  component  properties  derived  from  tests  and/or  analyses.  The 
various  methods  discussed  in  the  preceding  paragraphs  only  partially  satisfy  the  three  basic 
requirements.  Common  to  all  modal  synthesis  methods  discussed  in  the  preceding  is  the  complexity 
of  the  matrix  manipulations  involved  in  setting  up  the  coupling  and  assembly  procedure  to  obtain  the 
final  reduced  equation  system. 

The  major  limitation  on  the  use  of  the  existing  modal  synthesis  methods  is  their  lack  of  compatibility 
with  practical  experimental  procedures.  Although  various  types  of  component  dynamic 
representations  have  been  proposed,  only  those  requiring  normal  modal  properties  are  practical. 
Test  derived  modal  representation  is,  in  general,  incomplete;  the  component  normal  modes  obtained 
assuming  any  type  of  support  conditions  at  the  interfaces  are,  in  general,  different  from  those 
occurring  when  the  components  are  acting  within  the  compound  of  the  total  system.  Since  only  a 
limited  set  of  modal  data  is  obtainable,  the  interface  flexibility  is  not  adequately  modeled.  Depending 
upon  the  synthesis  method  used,  additional  information  is  therefore  usually  required  to  approximate 
the  effect  of  interface  condition  or  modal  truncation. 

Fixed  interface  mode  synthesis  methods  employ  static  deflection  shapes.  In  an  experimental  setup, 
constraining  interface  degrees  of  freedom  proves  impractical,  particularly  when  large  dimensions  or  a 
large  number  of  connection  points  are  involved.  Also,  damping  data  associated  with  static  modes  is 
unavailable.  For  these  reasons  the  free  interface  based  modal  synthesis  methods  are  best  suited  for 
achieving  test  compatibility.  These  methods  also  lend  themselves  to  accuracy  improvements  via  the 
artifice  of  interface  loading  or  by  augmenting  the  normal  modal  data  with  residual  flexibility  and 
inertia  effects  of  truncated  modes. 

5.1.4  Objectives  and  Scope 

The  objective  of  this  work  is  to  investigate  and  develop  component  dynamic  synthesis  procedures  and 
associated  computer  software  which  (a)  combines  component  dynamic  characteristics  obtained  from 
modal  tests  or  analyses  or  both;  (b)  accounts  for  the  effects  of  differences  in  interface  boundary 
constraints  of  the  component  structures  in  the  modal  test  and  in  the  comparison  of  total  structure; 
and  c)  reduces  inaccuracies  due  to  modal  truncation. 

In  view  of  the  above  objectives  and  the  assessment  of  existing  synthesis  methods  made  in  the 
preceding  section,  the  free  interface  modal  synthesis  methods  are  studied  further  in  this  work.  For 


completeness,  the  fixed  interface  and  the  discrete  element  representations  are  also  considered.  For 
certain  components  the  use  of  constrained  interface  conditions  may  be  unavoidable.  Structural 
components  such  as  panels,  stringers,  simple  masses,  vibration  control  devices,  etc.  are  conveniently 
input  to  the  synthesis  procedure  via  discrete  elements. 

A  principal  feature  of  the  work  developed  here  is  the  component  dynamic  model  reduction 
procedure  that  leads  to  an  exact  and  numerically  stable  synthesis.  In  order  to  affect  component 
coupling,  neither  the  specification  of  external  coupling  springs  nor  an  user-specified  selection  of 
independent  coordinate  is  required.  Existing  synthesis  procedures  suffer  from  these  drawbacks. 
Component  dynamic  models  considered  include  free-free  normal  modes  with  or  without  interface 
loading,  up  to  second  order  stiffness  and  inertia  connections  accounting  for  the  effect  of  modal 
truncation,  fixed  interface  modes,  and  also  the  physical  coordinate  components.  The  modal 
reduction  procedure  involves  interior  boundary  coordinate  transformations  which  explicitly  retain 
connection  interface  displacement  coordinates  in  the  reduced  component  dynamic  representations. 
Interior  coordinates  may  include  physical,  modal,  or  any  admissible  coordinates.  Components  in  this 
reduced  form  are  termed  "superelement"  because  they  are  a  generalization  of  the  conventional  finite 
elements  of  structural  mechanics.  The  problem  of  component  dynamic  synthesis  is  then  reduced  to 
the  assembly  of  the  superelement.  The  direct  stiffness  approach  and  all  subsequent  processing 
operations  of  the  finite  element  method  are  then  applicable. 

The  remainder  of  this  section  presents  the  formulation  implementation  aspects  of  the  superelement 
component  dynamic  synthesis  method.  Section  5.2  contains  a  discussion  of  various  types  of 
component  dynamic  models,  a  synopsis  of  coupling  methods  followed  by  a  detailed  derivation  of  the 
superelement  component  dynamic  equations.  The  derivation  of  component  dynamic  models  from 
measurements  is  described  in  Section  5.3.  A  verification  study  of  the  developed  technique  is 
presented  in  Section  6.  The  computer  software  implementation  of  the  superelement  synthesis 
method  is  described  in  the  Appendix  B. 

5.2  Superelement  Method 

In  order  to  develop  an  improved  component  dynamics  synthesis  procedure,  there  are  two  key  issues 
to  focus  on:  the  modeling  of  component  dynamics  and  the  coupling  of  component  coordinates.  As 
seen  in  the  review,  no  one  method  of  component  dynamic  modeling  is  shown  to  be  superior  to  any 
other.  The  methods  of  synthesis  developed  in  the  literature  use  one  or  other  type  of  component 
representation.  Aerospace  structures  involve  a  wide  variety  of  components  and  any  single  component 
dynamic  modeling  method  may  not  be  uniformly  suitable  to  all  the  components.  With  this  in  mind,  a 
generalized  synthesis  method  is  develop  which  permits  different  types  of  component  models  and  an 
associated  coupling  procedure.  In  Section  5.2.1  various  component  dynamic  models  are  described. 
In  Section  5.2.2  coupling  procedures  are  described,  illustrating  the  accuracy  and  computational 
efficiency  of  direct  matrix  assembly  procedures  to  affect  component  coupling.  A  basic  requirement  of 
the  superelement  method  is  discussed.  In  order  to  implement  the  direct  matrix  assembly  procedure 
for  components  described  in  any  general  coordinates,  the  necessary  transformations  are  derived  in 
Section  5.2.3.  The  coupled  systems  equations  are  derived  in  Section  5.2.4. 

5.2.1  Component  Dynamic  Model 

The  equation  of  motion  of  a  component  undergoing  motion  within  the  compound  of  a  total 
structural  system  may  be  expressed  as 
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[M]{i}  +  [C]{i}+[[X]+;[^{x}={/} 


(168) 


where  M,  C,  K  and  X  are  respectively  the  mass,  viscous  damping,  stiffness,  and  hysteretic  damping 
matrices  and  {x}  the  column  matrix  of  displacement  due  to  internal  forces  {F}  at  the  interface  degrees 
of  freedom  connected  with  adjacent  components.  External  forces  are  assumed  to  be  applied  only  at 
connect  degrees  of  freedom.  In  the  following  derivations  we  ignore  damping.  Further,  we  use  the 
same  notation  for  the  time  dependent  quantities  and  their  amplitudes,  that  is,  (x(t)}  =  {x}  eJ*,  {F}  = 
{Fje^,  to  minimize  the  number  of  symbols.  For  later  discussion  the  system  of  Equation  (168)  is 
partitioned  according  to  boundary  {x}fl  and  remainder  (x}0  degrees  of  freedom  as 


The  above  represents  a  total  of  n  (=  nB  +  n0  )  number  of  equations  where  nB  and  n0  are, 
respectively,  the  number  of  boundary  and  other  degrees  of  freedom.  The  boundary  degrees  of 
freedom  include  the  connection  interface  DOFs  as  well  as  any  other  DOF  that  need  to  be  carried  to 
system  synthesis  in  physical  coordinates.  Thus  the  grid  point  DOFs  associated  with  supports,  external 
excitations,  response  monitoring  stations,  concentrated  weights  and  attachment  to  other  components 
for  higher  level  synthesis  are  also  termed  boundary  DOFs. 

In  order  to  reduce  the  size  of  the  above  system  equation  and  also  to  obtain  the  model  from 
experimental  tests  the  displacements  {*}  of  the  component  are  represented  by  a  linear  combination 
of  a  finite  number  of  displacement  functions  as 


{*}  =  [*]{<!}  (170) 


where  {q}  is  the  vector  of  generalized  coordinates  and  [$]  is  the  matrix  of  generalized  displacement 
functions.  Using  the  above  transformation  in  Equation  (168)  leads  to  the  equation  of  motion  in 
generalized  coordinates  as 


[**],{«}  ♦[x]<{?}  =  [$r{F}  (171) 


where  : 

[M],  =  [fc]7*  [M]  [$]  :  Generalized  mass  matrix 

[X  ]f  =  [fcf  [X]  [$]  :  Generalized  stiffness  matrix 

Equation  (171)  defines  the  component  dynamic  model  in  generalized  coordinates.  This  definition  is 
very  general  in  the  sense  that  no  restriction  is  placed  on  the  coordinate  {q}  other  than  the 
requirements  that  the  columns  of  this  transformation  matrix  [$>]  be  linearly  independent  admissible 
shape  functions  whose  linear  combination  adequately  represents  the  deformed  configuration  of 
the  elastic  component.  In  the  following,  however,  in  view  of  the  objectives  of  this  work,  attention  is 
focused  on  the  use  of  component  normal  modal  coordinates.  Some  auxiliary  displacement  functions 
arc  also  used  in  order  to  improve  the  accuracy  of  component  dynamic  representation.  The  normal 
mode  and  auxiliary  displacement  functions  for  different  interface  fixity  conditions  are  described  next. 


5.2.1. 1  Free  Interface  Modal  Models 


The  generalized  displacement  functions  in  Equation  (170)  are  obtained  from  the  solution  of  the 
eigenproblem 

(172) 


where  the  boundary  forces  {F}B  in  Equation  (169)  are  set  to  zero.  The  displacement  functions 
include  rigid  body  modes  if  any,  and  the  flexural  modes  of  the  component.  The  generalized  mass  and 
stiffness  matrices  are  diagonal  with  coefficients  given  as 

Mr={«nM]{4 

*,={*}?■ [*]{*},  (173) 


The  equation  of  motion  of  the  modal  component  is 


rMj{q}+  %{q}={f] 


(174) 


In  practice,  only  a  partial  set  of  normal  modes  is  available  which  fails  to  adequately  represent  the 
component  dynamic  behavior.  The  two  models,  loaded  interface  modal  model  and  the  residual 
flexibility  modal  model  discussed  in  the  following,  approximate  the  effects  of  high  frequency 
truncated  modes  and  improve  the  accuracy  of  the  component  dynamic  representation. 

The  loaded  interface  modal  model  is  obtained  from  a  modified  eigenproblem  of  the  component.  The 
modifications  account  for  the  stiffness  and  inertia  loading  of  the  component  interfaces.  Let 
and  [A M]BB  be  the  connection  point  stiffness  and  mass  loading  matrices  respectively.  The  modified 
component  eigenproblem  becomes 

[K  loo  [K  ] OB 

[K  Ibo  \bb  +  [A/Cjosj 

The  computed  eigenvectors  {^}r  diagonalize  the  modified  stiffness  and  mass  matrices  and  satisfy  the 
orthogonality  relations  of  the  types  given  by  Equations  (173). 

The  residual  flexibility  model  accounting  for  the  effects  of  modal  truncation  is  derived  using  the 
residual  flexibility  attachment  modes  used  in  a  number  of  previous  works  [56'57.58,63,69,70]  jn  tjie 
following  a  derivation  of  the  model  including  inertia  as  well  as  flexibility  effects  is  given. 

Partitioning  the  complete  modal  matrix  [4>]  into  orthogonal  subsets  [$>]*  and  [$];j ,  where  [4>k  is  a  n 
by  mK  matrix  of  kept  normal  modes  and  is  a  n  by  mD  matrix  of  omitted  high  frequency  modes,  n 
and  mK  being  the  total  number  of  component  degrees  of  freedom  and  the  number  of  retained 
modes,  respectively,  mD  is  the  number  of  omitted  modes  (  mD  =  n  -  mK  ).  The  physical 
displacement.  Equation  (170),  can  be  written  as 


-A? 


\M  loo  [M  ] OB 

[M]bo  [[mw  +  [amw| 


(175) 
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{x}  =  [$1k-{<7}a-  +  [*fo{<7}D 


(176) 


where. 


[*]=[l*k  I  [*h\ 

Substituting  Equation  (176)  into  Equation  (171)  gives  the  equations  of  motion  for  generalized 
coordinates {q}  as 

{$}*♦  rAJ/f{<7}*  =  [*]£{/}  (177) 


{<i)D+  rAj0{<7}D  =  [*]£{/}  (178) 

where  modes  are  orthonormalized  with  respect  to  the  mass  matrix,  and  and  lAjD  are  defined  as 
per  Equations  (173)  corresponding  to  the  kept  and  detailed  modes  respectively.  For  harmonic 
motion  of  the  component  at  frequency  w,  solving  Equation  (178)  for  the  generalized  coordinate  {<7} 
and  substituting  the  result  into  Equation  (176)  leads  to 


{*}  =  [$!*{<?}*: + 


{/] 


(179) 


The  second  series  in  the  above  expression  represents  the  dynamic  flexibility  or  simply  the  residual 
flexibility  of  the  modes  [<&]D.  Approximation  for  the  residual  effects  of  the  truncated  modes  is  arrived 
at  indirectly  since  the  modes  and  frequencies  in  the  second  series  are  not  known.  Following  Rubin 
the  static  contribution  of  the  truncated  modes  is  obtained  by  neglecting  inertia  in  Equation 
(178)  and  solving  for{^}D  : 

{<1)d  =  ^  ml{f)  (180) 


Substitution  of  the  above  into  Equation  (176)  gives  component  displacements  as 


{*}  =  [♦]*{*}*  +  [*h  TAj*'1  [*]£{/} 


(181) 


Using  the  definition  of  the  generalized  stiffness 
rAj=[*™[*l 

the  total  static  flexibility  is  given  as 

i*r 


for  constrained  components,  and  [ G]  =  \A'fr[G]\A ]  for  free  components  where  [A]  is  a  matrix 
operator  which  filters  out  rigid  body  modes: 

Ml  =  [/]-[M][*MA/ ]*[*]£ 

[<&!/!  and  [M\r  being  respectively  the  matrix  of  rigid  body  modes  and  corresponding  mass  matrix. 
[G]  is  the  flexibility  of  the  component  subjected  to  arbitrary  statically  determinate  constraints.  Using 
the  partitioned  form  of  the  modal  matrix  [<&]  the  above  equation  becomes: 


[G ]  =  [*Lc  fAjK-1  [*]£  ♦  [*k  ^  mi 


Therefore  the  component  displacements  can  be  written  in  terms  of  the  known  modes  [$]*  and  the 
static  flexibility  as 


W  =  WkMk  +  [{G]-  [*k  ^  [*]S|  {/} 


{*}-[*k{?}*  +  [G]*{/} 


where, 


[GJ*-[[G ]-[*]*  fc*1 


represents  the  residual  flexibility  of  the  truncated  modes. 

Correction  for  the  inertial  effects  of  the  truncated  modes  is  obtained  using  Equation  (180)  to 
approximate  the  inertia  term  in  Equation  (178).  Thus, 


{Q}d  =  V  [t/]^[G  ]*]{/} 


Correspondingly, 


The  approximation  for  the  flexibility  and  inertia  of  the  truncated  modes  is  clearly  seen  in  the 
displacement  transformation  relations  given  by  Equations  (183)  and  (186).  These  transformations. 


however,  are  not  suitable  for  component  coordination  reduction  because  of  the  force  term  {F}.  A 
Rayleigh-  Ritz  representation  for  Equation  (183)  is  readily  obtained  by  considering  the  residual 
displacement  term  in  Equation  (181) 

{x}R  =  [G]R{f) 

This  equation  defines  residual  displacement  of  the  component  due  to  interface  force  {F}.  Comparing 
it  with  the  definition  of  an  attachment  mode,  this  equation  leads  to  the  definition  of  a  residual 
flexibility  attachment  mode  for  a  unit  force  imposed  at  the  interface  coordinate.  The  residual 
displacement  vector  is  the  column  of  residual  flexibility  matrix  corresponding  to  the  interface 
coordinate  degree  of  freedom.  Component  displacement  {x}*  can  therefore  be  expressed  as  a 
superposition  of  residual  flexibility  attachment  modes  [G  V- 


{*}*  =  \Gh{q}B  (187> 


where  {q}B  is  the  associated  generalized  coordinate  and  [G  la  is  the  (rtxnB  )  partition  of  the  (nx/i) 
residual  flexibility  vector  [G  ]#■  That  is, 

[Gk'[[c lew.,.]  <188) 

Rayleigh-Ritz  transformation  for  the  component  displacement  then  becomes 

{*}  =  [T]A<1} 


where, 

!  [G]*|  ,  (<7}=j^j  (189) 

Equation  (182)  governs  component  dynamic  equilibrium  in  the  generalized  coordinates  {#}  given  by 
Equation  (189)  with  the  modal  parameters  defined  as 


„  riAf,  [0] 

[Q]  (19°) 

n  =  [T]J[K][T]^  ^  091) 


where  1 M and  are  respectively  the  generalized  mass  and  stiffness  matrices  of  the  kept  normal 
modes,  and  [H  \BD  and  [G]sb  are  the  boundary  partitions  of  the  residual  mass  and  flexibility 
matrices. 


rA#k-[*S[^][*lr 

(192) 

rAk  =  [*]5Ml*k 

[cl»-[C]S[«][cl, 

The  following  orthogonality  properties  of  the  modes  are  used  in  the  above  derivation: 

[$]£  [Af  ]  [G  ]*  =  [G  ]£  [M  ]  [$]*  =  0 

(*ff[^HCla-[C]5^][*k»0  (193) 

K;|J[/n[GU  =  [*]orAJi5[^  =  [Gk 

It  should  be  noted  that  in  contrast  to  Equation  (183),  through  the  use  of  residual  flexibility 
attachment  modes,  the  Raieigh-Ritz  transformation  of  Equation  (189)  requires  the  measurement  of 
residual  flexibilities  at  connection  DOFs  only.  Furthermore,  the  transformation  of  Equation  (189)  is 
amenable  to  superelement  assembly  procedure  as  is  shown  later  in  this  section.  An  estimate  of  the 
improvement  in  accuracy  of  component  displacement  representation  with  the  introduction  of  the 
residual  flexibility  effects  is  obtained  as  follows.  Comparing  Equations  (179)  and  (184),  the  error  in 
displacement  given  by  the  two  equations  is  less  than  or  equal  to 


II 

(H-^) 


where  uK  is  the  frequency  of  the  highest  retained  mode  of  the  component. 


5.2.1.2  Constrained  Interface  Modal  Model 


A  constrained  interface  modal  model  is  derived  using  the  displacement  transformation  given  as 


[0]  [/] 


where  the  columns  of  the  [<!>]#  and  [4>]c  matrices  are  respectively  the  constrained  interface  normal 
modes  and  the  static  constraint  deflection  shapes  of  the  component  obtained  as  follows.  For  forced 
harmonic  motion  of  the  component  at  system  frequency  A,  the  component  equilibrium  Equation 
(169)  can  be  written  as 


[K  ]oo  [K  J OB 
[k\bo  \db 


[Af  ]oo  [Af  Jo# 
[Af  ]  bo  [A 1\bb 


Since  the  interior  DOFs  are  force  free,  their  motion  can  be  written  in  terms  of  the  boundary  DOFs  as 


Mo  =  '  loo  '  A2  [A/  looj  j[K  lop  •  A2  [M  JobJ{x  )b 


0%) 


The  static  constraint  modes  are  derived  by  ignoring  inertia  forces  on  the  interior  coordinates.  Thus, 


{*}o=[*lcW*  (197) 


where 

[*lc  =  ■  [K  loo  [K  1  OB 

[$)o  is  the  constraint  mode  transformation. 

Physically  a  static  constraint  mode  is  the  displacement  configuration  of  the  component  resulting  from 
a  unit  displacement  applied  to  one  boundary  DOF  while  all  other  boundary  DOFS  are  held  fixed. 

Inertia  of  the  interior  DOFs  is  accounted  for  by  considering  the  motion  of  the  component  with  all 
interface  DOFs  held  fixed  and  solving  the  eigenproblem 

[[/C]oo-A2[A/]oo]{jc}o={0}  (198) 

The  solution  of  this  eigenproblem  yields  the  matrix  of  fixed  interface  normal  modes  [$]0  having  the 
same  order  as  [ K]00  and  [M]00.  The  computed  eigenfrequencies  \ 2  are  those  of  the  component 
with  its  boundaries  fixed. 

The  complete  set  of  component  normal  modes  [<&]0  plus  the  static  constraint  modes  [tfk  provide 
the  means  to  transform  the  displacement  vector  {*}  from  geometric  coordinates  to  an  equivalent  set 
of  generalized  coordinates  {q }  as  given  in  Equation  (194),  where  [4>k  is  a  subset  of  [4>]0  chosen  so 
as  to  yield  a  reduced  order  model.  The  reduced  order  stiffness  and  mass  matrices  in  generalized 
coordinates  are  obtained  as 

(199) 


where 


rAk«[*ff[*][*k-  fA?j 

[m  hoc  =  [®]£  W  ]  [®k  =  f/w,j 

(200) 

[k] cc  =  [*]5  [X]  [4>Jc  =  [AW  +  [*ko  [*]c 

[m  ]cc  =  [®]c  [A/]  [$]c  =  [M]bB  +  [$]c  [A/]oo  [$]c  ♦  [®]c  [M] ob  ♦  [M\bo  [®]c 
[mkc  =  [M\bo  [$k  ♦  [®]c  [M]oo  [®k 
[m  lac  =  [m  ]kc 

Two  characteristics  should  be  noted  concerning  the  generalized  coordinates  of  the  constrained 
interface  model:  (1)  the  generalized  coordinates  associated  with  the  constraint  modes  are  simply  the 
boundary  displacement  coordinates;  and  (2)  the  generalized  coordinates  associated  with  their  normal 
modes  of  the  component  are  unique  to  the  component  and  do  not  require  compatibility  with  other 
components.  It  is  in  this  regard  that  the  fixed  interface  model  mimics  the  general  characteristics  of  a 
displacement  based  finite  element.  Assimilation  of  this  and  the  other  models,  namely  the  interface 
modal  models,  finite  element  models  in  the  system  synthesis  is  discussed  in  the  next  section. 

Procedures  for  constructing  the  component  dynamic  models  discussed  in  the  preceding  from 
experimental  test  data  are  summarized  later  in  this  Section  5.  Detailed  discussions  may  be  found 
elsewhere  l57-5®, 60,69, 72] 

5.2.2  Component  Coordinate  Coupling 

The  assembly  of  the  component  dynamics  models  and  the  enforcement  of  the  compatibility 
constraints  at  their  connection  interfaces  to  form  coupled  system  equations  is  the  key  step  in  a 
synthesis  procedure.  The  accuracy  and  efficiency  of  synthesis  depends  on  the  manner  in  which  these 
interface  constraints  are  enforced,  which  in  turn  depends  upon  the  type  of  basis  coordinates  used  in 
the  component  dynamics  representation.  In  the  following,  several  coupling  procedures  are  discussed. 
Observe  that  the  direct  matrix  approach  is  the  only  exact  and  most  efficient  coupling  procedure. 
Finally,  observations  are  made  regarding  the  form  of  the  component  dynamics  model  needed  to 
affect  synthesis  using  direct  matrix  assembly  approach. 

In  the  following,  the  coupling  of  two  components  described  in  physical  coordinates  (Equation  (168)), 
generalized  coordinates  (Equation  (171)),  and  mixed  coordinates  are  considered.  In  discussing  the 
generalized  coordinate  model  full  generality  of  the  basis  coordinates  is  permitted,  that  is,  the  columns 
of  the  modal  matrix  [$>]  in  Equation  (170)  are  not  required  to  be  orthogonal  vibration  modal  vectors, 
the  generalized  stiffness  and  mass  matrices  are  in  general  fully  populated. 

The  problem  of  component  coordinate  coupling  can  be  approached  in  basically  two  ways  depending 
on  how  the  compatibility  condition  at  connection  interfaces  is  enforced.  In  the  direct  or  implicit 
method,  the  component  matrices  are  partitioned  into  interior  and  boundary  DOFs  and  then 
assembled  in  the  manner  of  direct  matrix  assembly  approach  used  in  conventional  finite  element 
method.  Thus  the  coupling  system  governing  equations  are 
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where  the  vectors  and  matrices  are  partitioned  as 
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The  superscripts  identify  the  component,  and  the  subscripts  B  and  O  represent  respectively  the 
boundary,  and  the  other  than  boundary  degrees  of  freedom.  The  interface  displacement  and  force 
compatibility  conditions 


{X}1B  =  {X}1 

{/U  =  -{/}! 


(203) 


are  implicitly  satisfied  in  obtaining  the  coupled  system  equations  (Equation  (201)).  Note  that  the 
presence  of  boundary  DOFs  {X}%  in  the  component  representation  permitted  the  use  of  direct 
matrix  assembly.  In  the  event  the  components  are  described  in  terms  of  generalized  coordinates  the 
constraints  in  Equation  (203)  need  to  be  enforced  explicitly  as  shown  below. 

Consider  the  coupling  of  the  components  described  by  Equation  (171).  The  governing  equations  are 


[A/]}  [Of 
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[*]}  [0] 
[0]  [X]? 


(204) 


or 


[M]q{q}  ♦  [*], {*}  =  {/}*, 

together  with  the  interface  compatibility  constraints 

MK<?}={0}  (205) 

which  are  generalized  coordinate  analog  of  Equation  (203),  [A]  is  the  matrix  of  connection  point 
modal  coefficients,  and  {q}  is  the  uncoupled  system  coordinate  vector.  The  constraint  Equation  (205) 
implies  a  reduction  of  the  number  of  system  DOFs  by  the  number  of  constraints.  This  reduction  is 
denoted  by 


{q}  =  [T){q} 


(206) 
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where  m  is  coordinate  coupling  or  synthesis  transformation  and  {<?}  is  the  coupled  system 
coordinate  vector.  The  manner  in  which  the  constraints  Equation  (206)  are  inserted  into  the  system 
Equation  (205)  significantly  affects  the  accuracy  and  efficiency  of  the  synthesis  procedure.  Possible 
alternates  are  as  follows. 


5.2.2. 1  Lagrange  Multiplier  Method 


The  constraints  are  introduced  by  extending  the  system  of  Equation  (204)  through  the  use  of 

nodeless  variable  of  Lagrange  multiplier  %  ,  i  =1,2 . NB  ,  NB  being  the  total  number  of 

connection  DOF.  The  overall  system  equations  become 


[M],  [0| 
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Contrary  to  the  reduction  in  the  number  of  coordinates,  the  system  size  is  extended  by  the  number  of 
compatibility  constraints.  Also,  it  leads  to  some  extraneous  system  frequencies  and  mode  shapes. 


The  system  size  reduction  can  be  achieved  by  applying  Guyan  condensation  to  the  system  Equation 
(207).  This  involves  assumptions  concerning  the  inertia  and  mass  of  the  condensed  out  DOFs  and 
also  the  reduction  is  nonunique.  An  inappropriate  selection  of  master  DOFs  may  lead  to 
meaningless  results. 


5.2.2.2  Elimination  of  Coordinates 


The  rectangular  matrix  [A]  of  Equation  (205)  is  partitioned  into  an  invertible,  nonsingular,  square 
matrix  [A  ]D  and  the  remainder  [A  1/  ,  the  coordinate  vector  {q}  is  also  partitioned  likewise.  Then, 


{q)D  =  -[A]^[A}I{q}l 


and  therefore 


\[-[AU[A] 


{g}  =  [Tm 


which  leads  to  system  equations  in  {q}  coordinates  alone  at  the  expense  of  computations  implied  in 
Equation  (208).  Note  that  the  coefficients  of  matrix  [A]  are  the  connection  point  modal  coefficients, 
finding  an  invertible  partition  of  which  may  present  computational  difficulties. 


v.v.v.v 


/  / 


S.2.2.3  Orthogonal  Complement  Method 


The  coordinate  reduction  transformation  [T]  is  generated  by  means  of  an  orthogonal  complement  of 
[A].  The  orthogonal  complement  [0\  of  [A]  consists  of  the  eigenvectors  corresponding  to  the  zero 
eigenvalues  of  the  following  problem: 

[Af\A){q)-\{q)  (209) 

The  uncoupled  system  coordinates  may  then  be  represented  in  terms  of  the  coupled  coordinates  {<7} 
through  the  relation  {q}  =  [0\{q}-  The  size  of  the  above  eigenproblem  is  equal  to  the  total  number  of 
generalized  coordinates  in  the  system.  Further,  due  to  roundoff  errors  the  zero  eigenvalues  may  not 
be  discernible  from  small  eigenvalues. 

5.2.2.4  Elastic  Coupling  Element 

Between  each  DOF  to  be  coupled  an  artificial  coupling  spring  is  introduced  which  leads  to  coupled 
system  equation  as 

[M ],  {9}  ♦[[KWAKj  {?}  =  {/},  (210) 

where  [AK]  is  additional  stiffness  due  to  coupling  spring  given  as 

[AK]  =  [*]£[*  M*lc  (211) 


where 


r*i  Nfc  I°J 
(  b "  10]  l*lfc 


is  the  matrix  of  modal  coefficients  corresponding  to  coupling  DOFs.  [$]fc  and  [$]fc  are  the  row 
partition  corresponding  to  the  connection  DOFs  of  the  modal  matrices  for  the  components  A  and  B 
respectively  and  \K  \c  is  the  stiffness  of  the  coupling  spring.  The  nature  of  the  coupled  system  so 
obtained  may  be  altered  drastically  by  an  inappropriate  choice  of  the  coupling  spring  stiffness. 


In  view  of  the  above,  it  is  obvious  that  the  direct  matrix  assembly  procedure  of  coupling  component 
coordinates  is  the  most  efficient  and  exact.  All  other  methods  introduce  added  computations  and/or 
approximations.  In  order  to  be  able  to  use  the  direct  matrix  assembly  procedure,  it  is  necessary  to 
have  the  connection  DOFs  explicitly  represented  in  the  component  dynamics  model.  On  the  other 
hand,  the  dynamic  representation  given  in  Equation  (170)  is  more  practical  and  obtainable  from 
experimental  tests.  The  superelement  dynamic  reduction  procedure  developed  in  the  following 
section  makes  use  of  the  dynamic  model  given  by  Equation  (170)  and  reduces  it  to  a  form  amenable 
to  direct  matrix  assembly  procedure.  All  the  component  dynamic  models  described  in  Section  5.2.1 
are  reducible  to  this  form,  thus  leading  to  a  unified  treatment  for  any  general  component 
representation. 


5.2.3  Superelement  Model  Reduction 

In  this  section  the  mathematical  transformations  are  derived  for  reducing  the  component  dynamic 
models  to  a  form  involving  physical  DOFs  of  the  connection  interface  nodes  plus  some  additional 
DOFs  related  to  the  interior  nodes.  The  component  dynamic  models  in  this  reduced  form  possess  a 
structure  similar  to  that  of  the  displacement  finite  element  method  and  are  therefore  termed 
superelements.  The  superelement  components  can  be  assembled  to  other  similarly  reduced 
components  by  the  direct  stiffness  method,  leading  to  an  efficient  and  exact  component  dynamics 
synthesis.  Four  different  types  of  component  dynamic  models  are  considered  in  the  following:  (1) 
finite  element  model;  (2)  fixed  interface  modal  model  with  static  constraint  modes;  (3)  free  interface 
modal  model  with  residual  flexibility  attachment  modes;  and  (4)  free  interface  model  with  normal 
modes  alone.  The  superelement  reduction  procedure  can  also  be  applied  to  the  component  dynamic 
models  involving  other  than  normal  modes  of  the  component  as  long  as  certain  requirements  for 
matrix  partitioning  and  invertibility  are  satisfied,  the  model  type  4  above  includes  this  model  as  well. 
A  system  may  involve  any  combination  of  the  above  types  of  component  dynamics  model  since  the 
models  are  reduced  to  a  common  form  before  assembly.  Furthermore,  the  component  dynamic 
models  chosen  above  are  such  that  the  component  independence  is  preserved.  Thus  the 
superelement  method  leads  to  a  very  general  synthesis  procedure. 

The  finite  element  model  described  in  Equation  (168)  or  in  partitioned  form  by  Equation  (169) 
naturally  contains  the  connection  interface  DOFS  and  is  already  in  the  required  superelement  form. 

The  constrained  interface  modal  model  described  by  Equations  (194)  through  (201)  is  also  in  the 
required  superelement  form  since  the  generalized  coordinates  associated  with  the  constraint  modes 
are  the  physical  DOFs  of  the  connection  interface  nodes  as  seen  in  Equation  (197).  It  should  be 
noted  that  as  the  number  of  modes  in  [$1*  is  reduced,  the  transformation  relation.  Equation  (195), 
shrinks  to  just  the  static  constraint  modes  and  the  fixed  interface  model  degenerates  to  a  Guyan 
reduced  component.  Likewise,  as  the  number  of  kept  normal  modes  in  [<&]*  is  increased,  the 
Equation  (195)  approaches  an  exact  transformation.  Selection  of  these  normal  modes  accounts  for 
the  accuracy  of  the  constrained  interface  superelement. 

Next,  consider  the  free  interface  modal  model  with  residual  flexibility  attachment  modes.  The 
displacement  transformation  of  Equation  (189)  can  be  written  in  partitioned  form  as 

fro)  _  [*ko  \G  ]ob  J<7k| 

[$ko  [G  \bd  [QbJ 

To  obtain  the  superelement  reduction  transformation,  solve  the  lower  partition  of  the  equation  for 
{q}D.  Thus, 


{<7 )b  -  '[G]bb[^\kb{<I}k  +  [ G]bb{x)b  (213) 


and  therefore. 
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The  above  transformation  reduces  the  residual  flexibility  attachment  model  to  the  superelement 
form.  The  net  transformation  can  be  obtained  by  combining  Equations  (212)  and  (214),  leading  to 


where 


m  = 
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The  generalized  coordinates  {q}K  are  now  the  participation  factors  of  the  modified  modes 
(  [$1ko  -  [G  Job  [G  Jbb  [OJkb  )  The  internal  partitions  of  the  modified  normal  modes  are  in  effect  free 
interface  normal  mode  partitions  minus  a  set  of  constraint  modes  internal  partitions  and  therefore 
represent  modes  constrained  at  boundary  degrees  of  freedom.  Note  the  similarity  of  the 
transformation  of  Equation  (216)  with  that  of  Equation  (194).  It  can  be  shown  that  the 

partition  of  the  transformation  is  in  fact  the  matrix  of  residual  constraint  modes  ^731, 

which  represents  the  static  deflection  of  the  component  due  to  the  residual  portion  of  the  stiffness 
matrix  resulting  from  a  unit  displacement  at  each  of  the  boundary  degrees  of  freedom.  The 
superelement  equations  of  motion  are  obtained  by  transforming  the  Equation  (168)  leading  to 
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where 
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[7  Jbb  -  [G  Jbb  [H  Jbb  [G  Jbb- 


The  modification  of  the  component  generalized  mass  and  stiffness  matrices  by  the  contributions  from 
the  residual  flexibility  of  the  deleted  modes  is  clearly  seen  in  the  above  equation.  To  recapitulate, 
[M  \kx  and  are  generalized  mass  and  stiffness  matrices  respectively.  fA^  and  [H  ]BB  are  the 
generalized  stiffness  and  mass  associated  with  the  residual  flexibility  attachment  modes  obtained  from 
transformations  given  in  Equation  (192). 

In  the  event  only  a  truncated  set  of  free  interface  normal  modes  is  used  to  represent  component 
displacements,  the  superelement  reduction  transformation  is  obtained  from  a  partial  inversion  of  the 
modal  matrix  as  described  in  the  following.  Consider  the  free  interface  normal  mode  transformation 
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where  [$ko  and  [$>1,®  are  (n0  xmK  )  and  ( nBxmK  )size  partitions  of  the  component  modal  matrix 
corresponding  to  the  interior  and  boundary  degrees  of  freedom  {x}0  and  ,  respectively. 
Consider  the  lower  portion  of  Equation  (219), 

{x  }B  =  [$]kb  {q}k  (220) 

and  partition  [$k»  into  a  nonsingular  invertible  square  matrix  [$]kbi  and  remainder  [$k»2-  Thus, 

{x}b  =  (221) 

which  requires  that  the  number  of  kept  modes  be  greater  than  the  number  of  boundary  degrees  of 
freedom.  Solving  the  above  equation  for  {q}xi  gives 

{q}ki  -  -  [$1kbi  [$]kb2 {<7}/C2  +  [$1/cbi  {x}b  (222) 
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The  generalized  coordinates  {q}K2  are  called  reduced  modal  coordinates  since  they  are  associated 
with  modified  modes  as  seen  in  the  above  equation.  Combining  Equations  (219)  and  (224),  a  net 
transformation  from  physical  coordinates  {Y}  to  free  interface  normal  mode  superelement 
coordinates  can  be  written  as 
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where  [$1koi  and  [<&lif02  are  partitions  of  [$]ko  created  using  the  partitioning  information 
generated  in  Equation  (221).  Analogous  to  Ritz  transformation  for  constrained  interface  modal 
model  Equation  (194)  and  the  transformation  for  residual  flexibility  modal  model  Equation  (217), 
the  left  and  right  partitions  of  the  transformation  of  Equation  (225)  may  be  interpreted  as  expressing 
modal  matrices  of  modified  normal  modes  and  boundary  constraint  modes,  respectively.  The 
component  equations  of  motion  in  superelement  coordinates  ({g}2  ,  {x}B  )  are  obtained  from  pre- 
and  post-multiplication  of  Equation  (169)  with  the  transformation  [T]  of  Equation  (225). 
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or  alternately  using  Equation  (221)  and  transforming  the  modal  equilibrium  Equation  (171). 
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where  full  generalized  mass  and  stiffness  matrices  are  used  to  include  the  case  where  the  columns  of 
the  modal  matrix  in  Equation  (170)  are  not  orthogonal.  (Mass  and  stiffness  loaded  component 
modes  are  nonorthogonal  with  respect  to  the  original  unloaded  mass  and  stiffness  matrices,  for 
example.) 

It  is  to  be  emphasized  again  that  the  reduced  modal  representation  given  in  Equations  (224)  through 
(227)  is  obtainable  from  any  given  free  interface  modal  representation  provided  the  partitioning 
conditions  leading  to  Equation  (222)  can  be  met.  The  columns  of  modal  matrix  in  Equation  (219) 
need  not  be  normal  modes.  The  only  necessary  condition  is  that  they  be  linearly  independent, 
adequate  in  number,  and  be  from  a  complete  set,  i.e.,  a  linear  combination  of  them  should  be  capable 
of  representing  the  deformation  shapes  of  the  component  undergoing  motion  within  the  compound 
of  the  built-up  system. 

5.2.4  Superelement  Synthesis 

The  superelement  reduction  procedures  described  in  the  preceding  section  yield  component  models 
with  a  common  characteristic:  the  boundary  degrees  of  freedom  are  explicitly  retained  and  the 
internal  degrees  of  freedom  are  transformed  to  reduced  generalized  coordinates.  The  equations  of 
motion  of  a  superelement  component  can  be  expressed  as 


(228) 


or 


[M  )•'*>  {p  }<“>  +  [C  )(a>  {p  },a>  ♦  [K  l<a>  [p  }<°»  =  {/}'“> 


Each  component  substructure  of  a  given  built-up  system  can  be  expressed  in  the  above  format. 
Analogous  to  the  matrix  equation  of  a  displacement  based  finite  element  model,  the  generalized 
coordinates  may  be  interpreted  as  internal  DOFs.  Assuming  that  the  global  structure  as  partitioned 
into  N  components,  the  assembled  global  matrices  [M],  [C],  [K]  relating  the  global  displacement 
vector  {x}  and  the  global  force  vector  {F}  can  now  be  synthesized  by  the  direct  stifbiess  approach  from 
the  component  matrices  of  Equation  (228).  The  direct  stiffness  approach  postulates  equal  nodal 
displacement  and  nodal  equilibrium  of  forces  at  the  connection  interfaces  of  adjacent  components, 
thus 
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m 

No  compatibility  conditions  exist  among  the  reduced  generalized  coordinates  and  remain 

intact.  The  coupled  system  coordinate  vector  thus  becomes 


{^}r  =  |{9}(l)T.{9}(2)T'---.{<7}(W)T{^}^ 


so  that  the  global  mass,  damping,  and  stiffness  matrices  are  given  by 
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where  [y9](a)  is  the  portion  of  the  component  to  global  transformation  matrix  [/9]  corresponding  to  the 
a"*  component. 
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{ F)  being  the  vector  of  uncoupled  superelement  coordinates. 


m7 =•[{<? },1)T{*}*(1,T  {</}(2)T{*}*(2)T  •  ••{<?  },Ar)T{*}*w‘ 


and  {P]n[a)  representing  the  compatibility  conditions  for  the  degrees  of  freedom  at  the  boundaries  of 
the  component  a. 

The  built-up  system  mass,  damping,  and  stiffness  matrices  are  of  the  following  form: 
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where  the  coefficients  a{i (“>  are  the  mass,  damping,  or  stiffness  coefficients  matrices  given  in 
Equation  (228),  and  [A]  is  the  corresponding  global  matrix.  The  diagonal  submatrices  [a  ]„*“*  in 
Equation  (234)  correspond  to  the  reduced  generalized  coordinates  of  the  components  and  are 
uncoupled  from  other  components. 

[«li»-  E  i/?b(a)T  (a  Wu)  L9b(a) 

«=i 

The  submatrices  are  symmetric  and,  in  general,  fully  populated.  The  system  governing  equations  can 
be  expressed  as 


[mKxMCK*M*K*}  =  {/} 


(235) 


which  itself  is  in  a  superelement  form  since  the  system  basis  coordinate  {x}  contains  the  boundary 
degrees  of  freedom  explicitly.  This  is  particularly  useful  from  the  standpoint  of  coupling  the  system  of 
Equation  (235)  to  a  higher  level  superelement. 

Next,  the  recovery  of  component  displacements  from  the  system  displacement  vector  {x}  is 
considered.  The  process  is  essentially  back  substitution.  The  superelement  displacement 
({q}HT{*}0HT)  is  recovered  from  the  system  displacement  vector  {x}  using  Equations  (229)  and 
(230).  The  use  of  superelement  transformations  and  component  transformations  described  in 
Sections  5.2.3  and  5.2.1  leads  to  the  recovery  of  component  displacements.  The  software 
implementation  of  the  superelement  synthesis  method  is  described  in  the  Appendix  B. 


5.3  Summary  and  Conclusions 

This  Section  documents  analytical  investigations  and  software  developments  aimed  at  providing  an 
improved  procedure  for  component  dynamic  synthesis  from  test  and/or  analysis  derived  component 
dynamic  models. 

A  set  of  consistent  Ritz  transformations  was  derived  that  lead  to  an  exact,  efficient  and  unified 
procedure  for  coupling  component  dynamic  models.  A  broad  class  of  test  and/or  analysis  derived 
component  dynamic  models  were  considered  in  this  work.  These  dynamic  models  are  compatible 
with  the  state  of  the  art  experimental  modal  testing  as  well  as  analytical  procedures  and  permit 
improvement  of  synthesis  accuracy  through  the  inclusion  of  flexibility,  inertia  and  damping 
corrections  of  truncated  high  frequency  modes  in  the  component  dynamic  models.  The  synthesis 
procedure  developed  in  this  work  may  be  considered  as  a  generalization  of  the  Craig-Bampton 
method  to  include  free  interface  normal  modes,  residual  flexibility  attachment  modes,  loaded 
interface  normal  modes,  and  any  general  type  of  component  modes  or  admissible  shapes  that 
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adequately  represent  the  dynamics  of  a  component.  Several  existing  methods  such  as  MacNeal’s 
method,  and  Rubin’s  method  are  shown  to  be  the  particular  cases  of  the  generalization  presented  in 
this  work.  As  a  result  of  this  generalization,  the  different  components  of  a  built-up  system  may  be 
characterized  using  any  type  of  dynamic  model  which  is  most  convenient  regardless  of  the  manner  in 
which  the  other  components  are  characterized.  Furthermore  the  component  coupling  is  exact  and 
computationally  efficient;  no  artificial  coupling  element  or  the  user  specification  of  independent 
coordinates  is  required.  The  developed  procedure  is  implemented  in  a  stand  alone  computer 
program  COMSYN  in  Appendix  C. 

The  component  dynamics  reduction  method  developed  in  this  work  transforms  the  component 
dynamic  coordinates  to  the  superelement  coordinates  containing  the  physical  coordinates  of  the 
connection  interfaces  as  well  as  any  desired  noninterface  points.  As  a  result  of  this  reduction 
procedure  the  component  subsystems  take  the  form  of  a  finite  element.  It  is  therefore  possible  to 
obtain  system  synthesis  even  with  nonlinear  components.  The  addition  of  the  necessary  data 
handling  and  solution  algorithms  to  treat  nonlinear  components  will  greatly  enhance  the  capabilities 
of  COMSYN.  It  is  recommended  as  a  further  work. 


-73- 


6.  PRACTICAL  EXAMPLES 


Three  examples  are  presented  in  this  section  to  demonstrate  the  practical  application  of  the  theories 
outlined  in  previous  chapters. 

6.1  Example  1:  Structural  Dynamics  Modification  Using  System  Modeling 
Techniques 

This  section  presents  the  conditions  and  results  of  the  experimental  testing  and  modeling  of  a  H 
frame  completed  to  form  a  practical  analysis  and  comparison  of  the  techniques  investigated  in 
Sections  2  and  4.  Results  of  finite  element  analysis,  modal  modeling,  component  mode  synthesis, 
and,  impedance  modeling  using  frequency  response  functions  are  summarized  in  Tables  2  through  6. 
All  models  are  compared  to  experimental  results  as  a  basis  to  evaluate  and  compare  the  technique. 

6.1.1  Development  of  the  Experimental  Modal  Model 

The  test  structure  used  in  the  modal  analysis  is  an  H  frame  structure  (see  Figure  1)  constructed  of 
two  by  six  inch  rectangular  steel  tubing.  Bolted  to  each  leg  of  the  structure  is  a  ten  by  eleven  by  one 
inch  aluminum  end  plate.  The  combined  weight  of  the  structure  with  end  plates  is  approximately  250 
pounds.  This  structure  was  used  in  the  analysis  primarily  because  its  welded  steel  construction  closely 
approximates  a  true  linear  system,  a  basic  assumption  with  all  modeling  routines.  In  addition  ,  the 
simple  geometry  of  the  structure  is  easily  described  in  the  analytical  model. 

Preliminary  analysis  indicated  20  -  25  modes  existing  in  the  frequency  range  from  0-512  Hertz. 
Figure  2  shows  a  typical  driving  point  frequency  response  function  of  the  unmodified  structure. 
Based  on  the  modal  density  of  the  structure,  resolution  of  1  Hertz  at  this  frequency  range  was 
considered  sufficient  to  accurately  identify  each  of  the  modes. 

One  important  characteristic  of  this  structure  noted  during  the  preliminary  analysis  was  that  the 
vertical  and  lateral  modes  of  the  structure  were  highly  uncoupled.  Because  of  this  ,  three  complete 
multi-point,  burst  random  modal  surveys  were  performed  on  the  structure  using  dual  inputs  in  the 
vertical,  lateral  ,  and  skewed  directions.  Burst  random  excitation  was  used  to  minimize  the  leakage 
errors  associated  with  the  lightly  damped  structure. 

The  modal  parameters  from  the  skewed  input  test  were  extracted  using  the  polyreference  time 
domain  algorithm  74l  Information  from  this  skewed  test  generated  one  complete  modal  model  since 
it  contained  both  the  lateral  and  vertical  modes.  The  modal  parameters  from  the  lateral  and  vertical 
input  test  were  calculated  using  the  polyreference  time  domain  algorithm  and  merged  to  generate  a 
second  model  of  the  structure.  This  model  was  used  primarily  for  reference  in  the  analysis.  Table  1 
lists  the  frequency  ,  damping,  and  modal  vectors  descriptions  for  the  19  modes  of  the  baseline  data 


6.1.2  Rotational  Degrees  of  Freedom 

Because  the  aluminum  end  plates  were  known  to  remain  rigid  over  the  entire  frequency  range  of 
interest,  a  rigid  body  program  was  used  to  generate  a  least  squares  fit  of  the  rigid  plate  motion  using 
all  available  measurements.  In  addition  to  generating  enhanced  modal  vectors  of  the  rigid  bodies,  the 
program  was  used  to  determine  rotational  degrees  of  freedom  associated  with  each  rigid  body.  These 
rotational  degrees  of  freedom  were  then  stored  as  new  measurement  degrees  of  freedom  located  at 
the  center  of  gravity  of  each  plate.  When  incorporated  into  the  modal  model,  these  rotational 


degrees  of  freedom  improve  the  results  of  the  predicted  mode  shapes  by  adding  rotational  constraints  \ 

to  the  model.  Ignoring  these  rotational  degrees  of  freedom  tends  to  result  in  high  estimates  of  the 

predicted  natural  frequencies  (especially  for  those  modes  with  high  rotations  at  the  end  plates)  due  to 

lack  of  rotational  mass  moment  of  inertia  in  the  model.  Table  2  shows  the  results  of  mass 

modification  (the  addition  of  28  pound  steel  masses  bolted  to  the  end  plates  in  the  long  ends  of  the  i 

H  frame)  using  the  Snyder  method  both  with  and  without  the  rotational  information  included  in  the 

model.  i 

t 

t 

6.1.3  Hardware  Modifications 

i 

i 

The  remainder  of  this  section  presents  the  results  of  the  modeling  phase  of  these  practical 

applications.  The  procedure  and  detail  of  the  modifications  are  discussed  to  highlight  the  capabilities 

and  limitations  of  each  technique.  The  results  of  the  various  modeling  programs  are  compared  with  J 

the  results  of  the  experimental  tests  of  the  modified  structure.  In  addition,  a  finite  element  model 

was  constructed  using  SMS  FESDEC  finite  element  code.  This  model  was  tuned  to  match  the  1 

dynamic  characteristics  from  the  baseline  experimental  data.  Each  physical  modification  to  the 

structure  was  also  performed  in  the  analytical  model  and  these  results  are  included  in  the  analysis.  A  I 

percent  error  is  calculated  for  each  eigenvalue  by  dividing  the  difference  between  the  experimental 

retest  frequency  and  the  predicted  frequency  by  the  experimental  retest  frequency.  A  more  thorough 

comparison  would  include  a  comparison  of  the  damping  values  and  the  computed  correlation 

coefficients  between  the  experimental  and  predicted  eigenvectors. 

Four  modifications  were  implemented  on  the  H  frame  to  evaluate  the  effectiveness  of  the  various 
modeling  techniques.  All  of  these  modifications  to  the  structure  were  made  at  the  rigid  end  plates  to 
take  advantage  of  the  calculated  rotational  degrees  of  freedom  and  improved  modal  coefficients.  The 
experimental  tests  of  the  modified  structure  were  multi-reference  impact  tests.  Two  triaxial 
accelerometers  were  mounted  to  the  structure  and  the  frequency  response  matrix  was  generated  by 
roving  the  impact  to  generate  the  modified  data  base.  The  modal  parameters  of  this  modified  data 
base  were  extracted  using  the  polyreference  time  domain  algorithm. 

The  first  of  the  modifications  investigated  was  a  long  thin  steel  rod  with  negligible  bending  stiffness 
connected  across  the  two  end  plates  on  the  long  end  of  the  H  frame.  This  modification  was  modeled 
as  a  simple  scalar  stiffness,  or  spring  element  in  the  analytical  models.  Table  3  summarizes  the  results 
of  the  stiffness  modeling  results  (using  SDM,  SYSTAN  and  SMURF)  compared  with  the 
experimental  re-test  and  finite  element  results. 

The  second  modification  investigated  was  the  addition  of  28  pound  steel  masses  bolted  to  the  two 
end  plates  on  the  long  end  of  the  H  frame.  These  masses  were  modeled  as  lumped  mass  elements 
with  the  correct  inertia  properties  about  each  of  the  three  principle  axis.  These  lumped  mass 
elements  were  applied  to  the  model  at  the  center  of  the  end  plates  to  use  the  rotational  information 
previously  calculated  for  these  points.  Table  4  summarizes  the  results  of  the  mass  modification 
(using  SDM  and  SYSTAN).  Note  that  the  frequency  response  method  of  general  impedance 
modeling  (SMURF)  was  not  investigated  for  the  mass  modification,  since  only  a  single  constraint 
equation  can  be  modeled  at  one  time. 

The  final  two  modifications  involved  elements  that  required  a  more  detailed  description  of  their 
properties.  The  first  of  these  was  a  two  and  one  half  inch  by  one  half  inch  beam  having  mass  and 
stiffness  properties  in  all  directions  connected  across  the  long  end  of  the  H  frame.  The  last 
modification  involved  the  connection  of  a  ten  by  one  quarter  inch  steel  plate  across  the  two  end 
plates  on  the  long  end  of  the  H  frame.  Because  these  two  modifications  required  matrix  descriptions 
of  their  mass  and  stiffness  properties,  only  the  Component  Mode  Synthesis  (SYSTAN)  and  finite 
element  solutions  were  evaluated.  The  results  of  these  modifications  are  summarized  in  Tables  5  and 
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6.1.4  Discussion  of  Modification  Results 


In  general,  the  SYSTAN  results  compared  best  with  the  experimental  test  of  the  modified 
structure(s).  The  average  error  for  each  of  the  modifications  were:  stiffness  -  Z8  %,  mass  -  1.1  %, 
beam  -  3.0  %,  and  the  plate  -  3.6  %.  The  SMURF  technique  yielded  excellent  results  for  the  single 
modification  investigated.  The  average  error  was  1.7%.  AH  of  the  techniques  predicted  the 
modifications  within  the  realm  of  experimental  accuracy.  As  the  complexity  of  the  modification 
increased  the  results  deteriorated  somewhat.  This  application  has  demonstrated  that  all  of  the 
methods  used  produce  accurate  results  when  used  properly  with  good  modal  data  and  procedure 
used  to  obtain  accurate  results.  More  rigorous  applications  will  determine  the  upper  bounds  of  the 
capabilities  of  each  of  the  techniques. 


0 .0  HZ  (  0) 


500.0  Hz  (512) 


Figure  2.  Typical  Driving  Point  Frequency  Response  Function 
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TABLE  1.  H  frame  Baseline  Test  Results 


Mode 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


Frequency 

14.8 

24.1 

36.6 

55.3 

77.4 

158.5 

162.6 

169.5 
191.0 
193.9 

209.1 

297.5 

309.1 

312.6 
417.8 


Damping  (%) 
07 
0.1 
0.2 
0.1 
0.1 
0.1 
0.2 
0.2 
0.4 
0.1 
0.4 
0.3 
0.4 
0.1 
0.1 


16 

434.1 

0.1 

17 

438.1 

0.1 

_ Description _ 

1st  lat.  -  rigid  legs  anti. 

1st  torsion  -  legs  rigid 
1st  shearing 

2nd  lat.  -  1st  full  frame  anti. 
3rd  lat.  -  1st  full  frame  sym. 
4th  lat.  -  2nd  full  frame  anti. 
1st  vert.  -  1st  full  frame  sym. 
2nd  vert.  -  1st  full  frame  anti. 
1st  twist  long  leg  -  anti. 

5th  lat.  -  2nd  full  frame  sym. 
2nd  twist  long  legs  -  sym. 

1st  twist  short  legs  -  anti. 

2nd  twist  short  legs  -  anti. 

X  bar  fore  &  aft  bending 
6th  lat.  -  3rd  full  frame  sym. 


3rd  vert.  -2nd  full  frame  sym 


7th  lat.  -  3rd  full  frame  anti. 


TABLE  2.  The  Effects  of  RDOF  on  Modeling  Results  (Adding  Masses) 


Mode 

Exp. 

Snyders  Method  (SDM) 

Description 

<4- 

<4  w/RDOF 

<4  w/o  RDOF 

1 

10.50 

10.85 

10.87 

1st  lat-anti  rigid 

2 

19.40 

19.95 

19.95 

1st  torsion 

3 

30.20 

29.71 

29.79 

1st  shearing 

4 

52.80 

53.27 

53.35 

2nd  lat-1  anti 

5 

72.60 

72.02 

72.44 

3rd  lat-l  sym 

6 

122.40 

113.13 

189.62 

1  tw  long-anti 

7 

126.80 

117.77 

204.81 

2  tw  long-sym 

8 

139.40 

142.81 

148.60 

4th  lat-2  anti 

9 

146.10 

144.04 

144.51 

1st  vertical 

10 

159.70 

159.11 

159.38 

2nd  torsion 

11 

174.20 

179.14 

184.28 

5  Iat-2  sym 

12 

281.80 

287.22 

295.80 

1  tw  short-anti 

13 

299.90 

302.42 

308.60 

2  tw  short-sym 

14 

370.20 

387.12 

307.70 

6  lat-3  sym 

15 

382.50 

397.32 

413.61 

7  lat-3  anti 

16 

395.80 

400.62 

411.39 

3rd  ver-2  sym 

17 

343.40 

457.11 

467.64 

4  ver-2  anti 

18 

476.50 

478.96 

485.75 

8  lat-4  anti 
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TABLE  3.  Stiffness  Modification  Modeling  Results  (Steel  Thin  Rod) 


ode 

Exp. 

FESDEC 

SDM 

SYSTAN 

SMURF 

Description 

<A 

<*v 

%  error 

<4 

%error 

“V 

%  error 

<4 

%  error 

1 

24.2 

23.7 

2.1 

23.8 

1.7 

24.1 

0.4 

25.6 

5.8 

1st  torsion 

2 

36.6 

37.3 

1.9 

36.1 

1.4 

36.5 

03 

36.6 

03 

1st  shearing 

3 

49.9 

50.6 

1.4 

48.2 

3.4 

50.9 

2.0 

49.9 

2.6 

1  lat-1  anti 

4 

77.5 

78.7 

1.5 

64.5 

17.0 

77.4 

0.1 

77.5 

0.8 

2  lat-1  sym 

5 

93.4 

94.5 

1.2 

91.7 

1.8 

96.0 

2.8 

93.4 

6.8 

3  lat-2  anti 

6 

161.6 

167.8 

3.8 

159.8 

1.1 

162.0 

0.2 

161.6 

1.4 

4  lat-  2  anti 

7 

164.3 

169.5 

3.2 

162.7 

0.7 

164.0 

0.2 

164.3 

0.1 

1  ver-1  sym 

8 

168.1 

176.1 

4.8 

171.8 

2.2 

169.2 

0.7 

186.1 

0.9 

2  ver-1  anti 

9 

191.8 

198.3 

3.4 

193.7 

1.0 

194.0 

1.1 

191.8 

15 

5  lat-2  sym 

10 

211.8 

216.0 

2.0 

201.1 

6.1 

209.0 

1.3 

211.8 

0.9 

1  tw  l-sym 

11 

276.3 

289.9 

4.9 

292.4 

5.8 

267.0 

3.6 

276.3 

2.2 

1  tw  s-anti 

12 

308.9 

304.4 

1.2 

306.1 

0.6 

308.0 

3.2 

308.0 

3.2 

2  tw  s-sym 

13 

312.3 

336.4 

7.7 

312.0 

0.1 

312.0 

0.1 

312.3 

0.2 

X-bar  bend 

14 

417.7 

418.1 

0.1 

417.1 

0.1 

417.0 

0.2 

417.7 

0.2 

6  lat-3  sym 

15 

434.0 

442.6 

4.9 

430.9 

0.7 

429.0 

1.2 

434.0 

0.0 

3  ver-2  sym 

16 

453.1 

452.9 

0.4 

438.0 

3.3 

435.0 

4.0 

453.0 

3.5 

7  lat-3  anti 

17 

480.1 

490.5 

2.2 

475.2 

1.0 

480.0 

0.0 

480.1 

0.0 

4  ver-2  anti 

18 

489.6 

491.9 

0.5 

485.8 

0.8 

482.0 

1.6 

489.6 

0.0 

8  lat-  4  anti 

TABLE  4.  Mass  Modification  Modeling  Results  (Adding  Steel  Plates) 


Mode 

Exp. 

FESDEC 

SDM 

SYSTAN 

SMURF 

Description 

«V 

“V 

%  error 

<4 

%error 

<4 

%  error 

<4 

%  error 

1 

10.6 

10.6 

0.0 

10.9 

0.3 

10.6 

0.0 

N/A 

1  lat-anti  rigid 

2 

19.4 

19.7 

1.6 

19.9 

2.8 

19.4 

0.0 

1st  torsion 

3 

30.0 

31.5 

4.3 

29.7 

1.7 

29.2 

33 

1st  shearing 

4 

52.9 

53.6 

1.5 

53.3 

0.9 

53.3 

0.9 

2nd  lat-1  anti 

5 

72.7 

74.9 

3.2 

72.0 

0.8 

71.5 

13 

3  lat-1  sym 

6 

122.3 

130.3 

6.5 

113.1 

7.6 

113.0 

7.6 

1  tw  long-anti 

7 

126.8 

132.1 

4.2 

117.8 

7.1 

109.6 

14.0 

2  tw  long-sym 

8 

139.4 

148.3 

6.4 

142.8 

2.4 

139.5 

0.1 

4  lat-2  anti 

9 

146.2 

157.6 

7.4 

144.0 

1.4 

141.5 

3.1 

1  ver-1  sym 

10 

159.7 

154.9 

3.0 

159.1 

0.4 

155.7 

_  2-5  J 

2  ver- 1  anti 

11 

174.2 

183.2 

5.2 

178.1 

2.8 

164.3 

WEM 

5  lat-  2  sym 

12 

281.8 

290.4 

3.1 

281.2 

1.9 

277.4 

1.6 

1  tw  short-anti 

13 

299.7 

303.6 

1.2 

302.4 

0.8 

300.5 

0.2 

2  tw  short-sym 

14 

302.3 

328.1 

8.5 

305.5 

1.1 

294.9 

2.4 

X-bar  bending 

15 

370.1 

415.3 

12.0 

387.7 

4.6 

383.8 

3.7 

6  lat-3  sym 

16 

383.5 

435.5 

14.0 

397.3 

1.2 

399.2 

43 

7  lat-3  anti 

17 

395.8 

408.7 

3.3 

400.6 

1.2 

395.2 

0.2 

3  ver-2  sym 

18 

434.4 

461.5 

6.2 

457.1 

7.5 

441.6 

1.7 

4  ver-2  anti 

19 

476.4 

- 

- 

478.9 

0.5 

477.3 

0.2 

8  lat-4  anti 
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Mode 


FESDEC 


<4 

%  error 

38.0 

0.0 

40.6 

3.4 

51.7 

1.5 

77.7 

0.9 

126.9 

1.6 

160.0 

4.8 

161.9 

4.5 

179.1 

3.5 

192.5 

3.6 

208.7 

6.0 

291.4 

1.6 

304.0 

2.0 

333.8 

8.1 

373.1 

1.7 

426.6 

4.0 

429.4 

0.6 

459.5 

2.2 

460.4 

13.0 

487.8 

3.6 

543.0 

14.0 

SDM 


%error 


SYSTAN 

<4 

%  error 

38.4 

1.1 

37.6 

10.5 

52.2 

0.4 

78.1 

1.4 

SMURF 


Description 


1st  shearing 


1st  torsion 


1st  lat-short  anti 


2nd  lat-lst  sym 


1st  beam  bending 


2nd  beam  bending 


1st  ver-lst  sym 


2nd  ver-lst  anti 


3rd  lat-2nd  sym 


1st  tw  long-anti 


1st  tw  short-anti 


2nd  tw  short-sym 


X-bar  bending 


3rd  beam  bending 


2nd  tw  long-sym 


3rd  ver-2nd  sym 


4th  lat-3rd  sym 


5th  lat-2nd  anti 


4th  ver-2nd  anti 


6th  lat-3rd  anti 
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6.2  Example  2:  Normalization  of  Complex  Modes  using  PRA  Time  Domain 
Technique 

A  numerical  example  to  determine  a  set  of  normal  modes  from  a  given  set  of  complex  modes  using 
the  time  domain  principal  response  analysis  technique  is  demonstrated  as  follows: 

6.2.1  Numerical  Example 

To  demonstrate  the  validity  of  the  theory  outlined  in  Section  2.8.6,  a  four  degree-of-freedom  discrete 
system  is  used  here  to  simulate  a  test  structure.  Mass,  stiffness  and  damping  matrices  are  assumed  as 
follows: 


[3  0  0  01 

[200  -60  -80  -401 

3  -0.9  -0.6  -1 

0  2  0  0 

[*]= 

-60  340  -120  -120 

[C]= 

-0.9  3  -0.8  -0.5 

0  0  10 
0  0  0  2 

-80  -120  800  800 
-40  -50  -200  -200 

-0.6  -0.8  3  -0.6 
-1  -0.5  -0.6  2.5 

Complex  eigenvalues  and  eigenvectors  are  solved  and  listed  in  Tables  7  through  10.  Real-normalized 
modes  using  several  subsets  of  the  identified  complex  modes  are  also  presented  in  Tables  7  through 
10. 

From  this  example  of  a  nonproportionally  damped  discrete  system,  it  indicates  that  a  set  of  real- 
normalized  modal  vectors  can  be  obtained  using  the  proposed  method.  Nevertheless,  unless  all 
modes  are  included  in  the  analysis,  the  results  are  subject  to  modal  truncation  errors. 


TABLE  7.  First  Vibration  Mode 


First  Mode 

Complex  Mode 

f=  1.1598  Hz 

S  =  4.79  % 

Normal  Mode 

f  =  1.1604  Hz 

magnitude 

analytical 

include 

include 

include 

include 

phase. 

normal  mode 

all  4  modes 

modes  1,2,3 

modes  1,2,4 

modes  1,2 

0. 

100. 

100. 

100. 

100. 

3.66 

37.067 

37.067 

37.050 

37.051 

3.52 

18.825 

18.826 

18.812 

18.830 

18.809 

8.079 

7.33 

8.058 

8.058 

8.071 

8.063 

8.071 

Second  Mode 

Complex  Mode 

Normal  Mode 

f  =  2.0472  Hz 

f 

=  2.04502  Hz 

?  =  6.061  % 

analytical 

include 

include 

include 

include 

magnitude 

phase,0 

normal  mode 

all  4  modes 

modes  1,2,3 

modes  2,3,4 

modes  1,2 

26.473 

-173.34 

-26.262 

-26.263 

-26.280 

-26.475 

-26.227 

100. 

0. 

100. 

100. 

100. 

100. 

100. 

18.047 

2.06 

18.047 

18.047 

18.057 

18.066 

18.075 

7.807 

3.88 

7.794 

7.794 

7.801 

7.831 

7.818 

TABLE  9.  Third  Vibration  mode 


Third  Mode 

Complex  Mode 

Normal  Mode 

f  =  3.8228  Hz 

f  =  3.8237  Hz 

(  =  3.134  % 

analytical 

inc'ide 

include 

include 

include 

magnitude 

phase. 

normal  mode 

all  4  modes 

modes  1,2,3 

modes  2,3,4 

modes  3,4 

6.210 

-167.34 

-6.027 

-6.027 

-5.977 

-6.228 

-6.226 

17.322 

-174.47 

-17.236 

-17.236 

-17.220 

-17.349 

-17.368 

77.950 

-6.86 

78.318 

78.318 

77.980 

78344 

78.381 

100. 

0. 

100. 

100. 

100. 

100. 

100. 

TABLE  10.  Fourth  Vibration  Mode 


Fourth  Mode 

Complex  Mode 

Normal  Mode 

f  =  4.7423  Hz 

f 

=  4.75128  Hz 

S  =  5.0067  % 

analytical 

include 

include 

include 

include 

magnitude 

phase,0 

normal  mode 

all  4  modes 

modes  1,2,4 

modes  2,3,4 

modes  3,4 

2.378 

-2.413 

-2.363 

-2.367 

6.882 

-6.813 

-6.768 

-6.853 

100. 

■f 

100. 

100. 

100. 

100. 

100. 

40.241 

^ 1 

-40.552 

-40.552 

— 

-40.220 

-40.547 

-40.564 
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6.3  Example  3:  Superelement  Dynamics  Synthesis  Technique 

This  section  presents  the  results  of  a  numerical  study  conducted  in  order  to  demonstrate  the 
superelement  dynamics  synthesis  procedure  developed  in  Section  5. 

Figure  3  shows  a  discrete  spring-mass  damper  system  and  its  partitioned  components.  The 
component  A  represents  a  simple  model  of  an  aircraft,  while  that  labeled  B  represents  a  model  of  a 
store  to  be  attached  at  the  tip  of  the  wing.  The  following  types  of  dynamic  models  are  employed  in 
the  study. 

Component  A 

1.  Free  interface  normal  modes  followed  by  a  transformation  to  superelement  coordinates. 

2.  Free  interface  normal  modes  plus  a  residual  flexibility  attachment  mode  followed  by  a 
transformation  to  superelement  coordinates. 

Component  B 

1.  Free  interface  normal  modes  followed  by  superelement  coordinate  transformation. 

2.  Physical  coordinate  model. 


Discrete  Spring-Mass  Damper  System. 

^  I -  ^  ku  *SA  ■»  H  "a 

*1*  I -  ^ - ett  H 

3  i—  j—  ?”’■  I —  if —  1— 

The  Separated  Components  of  the  Unconstrained  System. 


mlA  =m%A  =  0.1 
klA  =  60.000 
m  lJ3  =  0.2 


=  36 

A:*,  =  120,000 

TT\  2JJ  -  0.4 


mx «  =  100 
ifcM  -  100.000 

ki£  =  2,000 


=  40 

=  50,000 


Figure  3.  Discrete  Spnng-Mass  System 
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For  comparison  purposes  the  problem  is  also  solved  using  existing  MacNeal  Method  and  Rubin 
Method  of  synthesis.  The  results  of  the  synthesis  are  shown  in  Table  11.  The  exact  results  are 
obtained  by  solving  the  eigenprobiem  of  the  built-up  system  without  partitioning. 


TABLE  11.  Direct  and  Synthesized  System  Natural  Frequencies  (rad/sec) 


Mode  No 

Exact 

Superelement 

Component 

Synthesis 

MacNeal 

Method 

Rubin 

Method 

Case  I 

Case  II 

Case  II 

Case  II 

1 

0 

0 

0 

0 

0 

2 

52.417 

52.417 

52.418 

52.419 

52.418 

3 

69.348 

69348 

69.829 

69.886 

69.829 

4 

72.818 

72.818 

291.697 

504.67 

291.697 

5 

418.029 

418.03 

6 

775.678 

775.678 

Case  I:  (1  rigid  body  +  4  elastic  modes)  of  Component  A 

Case  II:  (1  rigid  body  +  1  elastic  mode  +  1  residual  flexibility  mode)  of  Component  of  B 


The  superelement  synthesis  using  complete  mode  sets  of  the  components  leads  to  exact  system 
synthesis  as  expected.  Superelement  method  using  severely  truncated  component  mode  set  along 
with  residual  flexibility  attachment  modes  predicts  system  frequency  with  an  accuracy  better  than  that 
of  the  MacNeal  method  and  equal  to  that  of  the  Rubin’s  second  order  method. 

It  must  be  noted  that  in  this  study  different  types  of  component  dynamics  models  are  synthesized  by 
virtue  of  the  superelement  formulation,  and  that  the  accuracy  of  system  synthesis  is  entirely  governed 
by  the  type  of  component  dynamic  models  and  not  by  the  coupling  procedure.  No  further  numerical 
exercises  are  therefore  undertaken  to  study  the  accuracy  of  synthesis  with  respect  to  component 
modes. 


-88- 


A.- 


TTscrk.  wx'jyv^v^rvvv'jv,jv\r^^jteV^W7\3rw?wrwnK  *m  tiv^wy^vx^vttw^jvuvvvwuvvvwiv'.^tx^  k_*  v_*  \i*v*u  ^  ir jyv^’ A^^wrr^TTTTv 


7.  CONCLUSIONS 

Experimental  modal  analysis  developed  in  the  past  decade  can  provide  a  valid  data  base  used  in  the 
application  of  system  modeling  techniques.  The  success  of  applying  system  modeling  techniques  in 
improving  the  engineering  quality  of  the  industrial  products  through  a  design  cycle,  is  dependent  on 
the  quality  of  the  experimental  data,  and  the  accuracy  of  the  system  modeling  algorithm  used  to 
predict  the  altered  system  dynamics  of  a  structure  or  combined  structure(s). 

Generally  speaking,  all  system  modeling  techniques,  which  include  modal  or  impedance  modeling 
method,  sensitivity  analysis,  and  the  component  mode  synthesis  method,  can  predict  satisfactory 
results  if  a  complete  and  perfect  experimental  model  can  be  obtained  from  testing  and  used  as  data 
base  for  system  modeling  predictions.  In  reality,  there  exists  many  uncertainties  and  difficulties  in 
obtaining  a  complete  modal  or  impedance  model  representing  a  physical  structure.  Difficulties  in 
simulating  actual  boundary  conditions  in  the  testing  laboratory,  lack  of  rotational  degrees  of  freedom 
measurement,  incomplete  modal  model  due  to  limited  testing  frequency  range,  nonlinearities  existing 
in  the  structure  under  test,  scaling  errors,  and  mode  overcomplexity,  could  seriously  affect  the  quality 
and  completeness  of  the  experimentally-derived  modal  or  impedance  model.  These  deficiencies  in 
obtaining  a  reliable  and  complete  experimental  model  make  the  system  modeling  technique  a  much 
less  powerful  tool  in  the  application  of  engineering  design.  In  other  words,  currently,  the  weakness  of 
applying  system  modeling  techniques  comes  from  those  limitations  and  uncertainties  to  obtain  a 
desired  modal  or  impedance  model  of  physical  structure(s). 

At  the  present  time,  many  efforts  have  been  dedicated  by  researchers  to  overcome  those  deficiencies 
in  obtaining  a  desired  experimental  modal  or  impedance  model,  such  as  the  development  of 
rotational  transducers.  Further  research  and  practices  are  still  needed  to  develop  a  well-defined 
engineering  procedure  and  criterion  to  make  the  use  of  the  system  modeling  techniques  a  more 
powerful  and  reliable  tool  in  engineering  practices. 
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braces  enclose  column  vector  expressions 

row  vector  expressions 

brackets  enclose  matrix  expressions 

complex  conjugate  transpose,  or  Hermitian  transpose,  of  a  matrix 

inverse  of  a  matrix 
generalized  inverse  (pseudoinverse) 
size  of  a  matrix:  q  rows,  p  columns 
diagonal  matrix 
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complex  conjugate 
Fourier  transform 
inverse  Fourier  transform 
Hilbert  transform 
inverse  Hilbert  transform 
natural  logrithm 
Laplace  transform 
inverse  Laplace  transform 

complex  number:  real  part  "Re",  imaginary  part  "Im” 
first  derivative  with  respect  to  time  of  dependent  variable  x 
second  derivative  with  respect  to  time  of  dependent  variable  x 
mean  value  of  y 
estimated  value  of  y 

summation  oM,  5,  from  i  =  1  to  n 

partial  derivative  with  respect  to  independent  variable  "t” 
determinant  of  a  matrix 
Euclidian  norm 


Roman  Alphabet 


Apqr 

residue  for  response  location  p,  reference  location  q,  of  mode  r 
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damping 

COH 

ordinary  coherence  function! 

COH * 

ordinary  coherence  function  between  any  signal  i  and  any  signal  kt 

COHn 

conditioned  partial  coherence! 

e 

base  e  (2.71828...) 

F 

input  force 

r, 

spectrum  ofqth  reference! 

GFF 

auto  power  spectrum  of  reference! 

GFF„q 

auto  power  spectrum  of  reference  qf 

gff ■* 

cross  power  spectrum  of  reference  i  and  reference  k! 
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reference  power  spectrum  matrix  augmented  with  the  response/reference  cross 

power  spectrum  vector  for  use  in  Gauss  elimination 

cross  power  spectrum  of  response/referencet 

auto  power  spectrum  of  responset 

auto  power  spectrum  of  response  pt 

impulse  response  function* 

impulse  response  function  for  response  location  p,  reference  location  q  * 
transfer  function* 

frequency  response  function,  when  no  ambiguity  exist,  H  is  used  instead  of  H  (w)* 
frequency  response  function  for  response  location  p,  reference  location  q,  when  no 
ambiguity  exist,  //M  is  used  instead  of  H„iu)f 

frequency  response  function  estimate  with  noise  assumed  on  the  response,  when  no 
ambiguity  exist,  Hx  is  used  instead  of  H±  (w) t 

frequency  response  function  estimate  with  noise  assumed  on  the  reference,  when  no 
ambiguity  exist,  H2  is  used  instead  of  H2i^ 

scaled  frequency  response  function  estimate,  when  no  ambiguity  exist,  Hs  is  used 
instead  of  f/s(w)f 

frequency  response  function  estimate  with  noise  assumed  on  both  reference  and 
response,  when  no  ambiguity  exist,  Hv  is  used  instead  of  //v(w) t 
identity  matrix 

stiffness 

modal  participation  factor 
mass 

modal  mass  for  mode  r 
multiple  coherence  functiont 
number  of  modes 
number  of  references  (inputs) 
number  of  responses  (outputs) 
output,  or  response  point  (subscript) 
input,  or  reference  point  (subscript) 
mode  number  (subscript) 
residual  inertia 
residual  flexibility 
Laplace  domain  variable 
independent  variable  of  time  (sec) 
discrete  value  of  time  (sec) 
tk  =kAt 
sample  period 

displacement  in  physical  coordinates 
response 

spectrum  of p*  response* 

Z  domain  variable 


Greek  Alphabet 

S(t)  Dirac  impulse  function 

A /  discrete  interval  of  frequency  (Hertz  or  cycles/sec) 

At  discrete  interval  of  sample  time  (sec) 

e  small  number 

rj  noise  on  the  output 

Ar  r*  complex  eigenvalue,  or  system  pole 

I 

i 


-> ->  ->  -j>  ~j>  ->'  v-' -j’rsrj'r*wss.~J,.~S-~-~rf"7-r; 


Ar  =<rr  +  juT 

diagonal  matrix  of  poles  in  Laplace  domain 

noise  on  the  input 

variable  of  frequency  (rad/sec) 

imaginary  part  of  the  system  pole,  or  damped  natural  frequency,  for  mode  r 
(rad/sec) _ 

< 4  =  firV/  1  -?r2 

undamped  natural  frequency  (rad/sec) 

+  “V2 

scaled  p*  response  of  normal  modal  vector  for  mode  r 
scaled  normal  modal  vector  for  mode  r 
scaled  normal  modal  vector  matrix 
scaled  eigenvector 

scaled  ptt  response  of  a  complex  modal  vector  for  mode  r 
scaled  complex  modal  vector  for  mode  r 
scaled  complex  modal  vector  matrix 
variable  of  damping  (rad/sec) 

real  part  of  the  system  pole,  or  damping  factor,  for  mode  r 

damping  ratio 

damping  ratio  for  mode  r 


vector  implied  by  definition  of  function 


APPENDIX  A:  TEST  DERIVED  COMPONENT  DYNAMICS  MODELS 


The  appendix  served  as  a  supplement  to  the  theory  developed  in  Section  5,  describes  the  construction 
of  the  fixed  interface  modal  model  and  the  residual  flexibility  modal  model  from  experimental  test 
data  Also  given  is  the  derivation  of  free  interface  modal  model  using  measured  loaded  interface 
normal  modes. 

A.1  Fixed  Interface  Modal  Model 

Consider  the  component  Equation  195.  The  component  is  tested  with  =0  and  the  constraint 
modal  matrix  [$]  results  for  which: 

[*F[Klool*]«  r*fJ  (A-l) 

[*F[M]00[*]=  % 

Multiplying  by  [<&F  >n  the  upper  partition  of  Equation  195  and  substituting  the  transformation 
{A”}  =  [<&]{<?}  results  in 

[[KlS«[«j  KW  ]W  W  ’2 

The  modal  matrix  [4>],  modal  mass  rAf,Jt  and  modal  stiffness  rAfJ  are  readily  available  from  constraint 
modal  test,  but  the  stiffiiess  [K\bB  and  (A ]OB  must  be  derived  from  additional  static  testing.  From 
Equation  (A-2)  for  (x}B  =  {0} 

1  ob  [$]  {  <7  }  =  {PlB  ( A-3) 


The  matrix  on  the  left  hand  side  of  Equation  (A-3)  is  made  up  of  reaction  forces  for  each  mode.  This 
expression  is  expanded  out  into: 


[K\Ib  [*]{<?}  =  Wi.{Kk  •  • ..{/?} J  =  [R]  (A-4) 


The  matrix  [R]  is  the  modal  reaction  force  matrix  and  is  measured  directly  by  placing  a  force 
transducer  at  the  connect  DOF  during  the  modal  test. 

The  other  unknown  [K)BB  is  found  by  substituting  the  modal  reaction  force  matrix  [R]  and 
rearranging  Equation  (A-2)  into: 


r*,;1  [-rV  [R  l7] 
[(*1  rV]  [^W-[R]r^j[R]' 


(A-5) 


the  term  A2  represents  the  forces  on  the  unconstrained  DOFs  of  the  component.  For  static  case 
A-0  and  the  static  stiffness  of  the  boundary  DOF  is  obtained  from  the  lower  partition  of  Equation 


-5 


[*W  =  [*W-[*]r*«jW 


If  the  component  restraints  are  statically  non-redundant,  that  is,  just  sufficient  to  prevent  rigid  body 
motion,  then  there  is  no  stiffness  among  the  boundary  DOF  and  [K]sB  '  =0. 


[*]»-[*]  By1  [*F 


the  boundary  stiffness  can  be  obtained  from  the  modal  stiffness  and  modal  reaction  forces  alone. 

For  redundantly  restrained  components,  [K\bB  '  is  not  zero  and  therefore  some  additional  static  test 
or  analysis  is  required.  Consider  the  partitions  of  { x}B  nonredundant  subset  {x}c  and  remainder 

{*}/■ 


With  the  coordinates  {x}c  restrained,  the  static  deflection  due  to  some  applied  force  {f)j  of  the 
remainder  {x)f  coordinates  can  be  expressed  as  the  sum  of  the  rigid  body  motion  and  elastic 
deformation  of  the  component: 


{x),  =  [G){f),+  [T]{x}c 


where  [G]  is  the  static  flexibility  and  [T]  is  the  rigid  body  transformation.  Also  from  static 
equilibrium,  the  forces  applied  at  {x)j  must  be  balanced  by  forces  at  {x}c;  therefore 


[TV{f),  =  {f)c 


Rearranging  Equations  (A-8)  and  (A-9)  results  in  the  following  equations 


[G]*1  [-[G]-1^ 

[-[7T[cri]  [-[rncr1!: 


(A-10) 


The  square  matrix  in  Equation  (A-10)  is  the  required  [/f  which  completes  the  definition  of  the 
constrained  interface  modal  model. 

Therefore,  to  construct  the  constrained  interface  modal  model  following  information  must  be 


obtained  from  tests:  fA/fJ  [<I>] ,[/?],  [G]'1,  and  [7].  The  first  four  are  obtained  from  modal 
testing.  The  matrix  [G  must  be  measured  from  static  tests  using  a  nonredundant  subset  of  {x}B 
constrained.  The  rigid  body  transformation  matrix  [T]  can  be  obtained  from  geometry  alone. 

\2  Residual  Flexibility  Modal  Model 

The  dynamic  flexibility  G„  is  the  response  at  coordinate  q  due  to  excitation  at  coordinate  p.  It  can 
be  expressed  in  terms  of  the  modal  stiffness  k,  ,  mass  m,  ,  and  mode  shape  vector  {<j>}r  by  the 
following  summation  over  the  number  of  modes: 


E 


r=lkr-X2mr  T=mKtlkr -A2mr 


(A-ll) 


where  the  subscript  r  denotes  the  r*  mode  and  mK  denotes  the  number  of  retained  modes.  The 
second  term  on  the  right  hand  side  in  the  above  equation  is  called  the  residual  flexibility  due  to 
modes  deleted  in  a  simulation.  If  the  magnitude  of  the  forcing  frequency  w,  is  within  the  range  of 
lower  frequencies  of  the  component  containing  the  number  of  modes  retained  for  synthesis, 
o><<uv;r>n*  +  l,  a  first  order  approximation  to  the  flexibility  of  omitted  modes  is  obtained  as 

00  A.  A „ 

r-  vi  WpY*j 

1 

»  =  **  + 1 

Experimentally  the  term  GBpt  can  be  obtained  by  exciting  the  component  at  DOF  q,  measuring  the 
response  at  DOF  p,  and  subtracting  the  first  sum  in  Equation  (A-ll),  which  can  be  evaluated  after 
the  modal  analysis  has  determined  the  properties  of  the  first  mk  modes. 

The  residual  mass  term  [H]BB  in  Equation  192  must  however  be  obtained  from  the  knowledge  of  the 
residual  flexibility  and  the  component  mass  matrix  as  no  method  has  yet  been  demonstrated  for  its 
determination  from  direct  measurements.  Note  that  the  computation  of  the  residual  mass  requires 
the  full  mass  matrix  of  the  component  in  physical  coordinates  and  also  the  full  nynB  partition  [G]* 
of  the  residual  flexibility  matrix.  This  includes  residuals  for  both  boundary  [G  ]*»  and  other  [G  )OB 
DOFs. 

A.3  Inertia-Loaded  Interface  Modal  Mode! 

Let  the  experimentally  measured  modal  parameters  of  the  loaded  interface  component  be  stiffness  fc,, 
mass  mr,  modes  [<t>],  frequency  04,  and  [AAf]  be  the  matrix  of  added  mass.  The  procedure  to 
analytically  unload  the  component  and  obtain  the  resulting  modal  parameters  is  as  follows. 

Consider  the  equations  of  motion  of  loaded  component  in  modal  coordinates 


^{<7}+  %{<]  }  =  {(?} 


(A-12) 


where  {q}  is  the  generalized  coordinate  such  that 


A-3 


% 


W  =  [*]{<?> 


(A-13) 


{x}  being  the  physical  displacement  vector  of  the  component.  The  equation  of  motion  of  the  added 
mass  can  be  written  as 


[A M  ]{*}„={/}„ 


(A-14) 


where  {x}„  is  a  subset  of  the  component  boundary  DOF  {*}B  where  extra  masses  are  attached 
during  test,  and  {F}0  is  the  force  on  the  added  mass  from  the  component  at  the  attachment  points. 
Analytical  unloading  of  the  component  is  accomplished  by  combining  the  dynamical  equation  of  the 
loaded  component  and  the  negative  of  Equation  (A-14).  Thus, 


-  [AM  1  fO]lp«l  .  [[0]  [0]lpa]  ['fa 

[°]  siicM  KIW  ic/ 


Displacement  compatibility  at  the  attachment  interface  requires 


{*}a  =l$]u{fl} 


(A-15) 


(A- 16) 


where  is  the  partition  of  [$]  corresponding  to  the  attachment  degrees  of  freedom.  Thus, 


(A-17) 


The  use  of  Equation  (A-17)  in  Equation  (A-15)  leads  to  the  equation  of  motion  of  unloaded 
component  as 


[  fmn  -  A2  [*![„  [AM ][*] ,„]  {q }  ♦  ^  {</}=[-  [*]fc {F}0  +  { Q }]  =  0 


Solve  the  above  eigenvalue  problem 


[  \  -  A?  ( ^  - 1*]^  [A M  ]  [*]„)]  U>}r  =  {0} 


where  \  and  {#}r  are  respectively  the  eigenvalue  and  mode  shapes  of  the  unloaded  component.  The 
modal  properties  of  the  unloaded  component  are  then 

['i'l7  U:  ['I']  =  Generalized  stiffness 


[♦I7  (AM  ]  [$]a|  [i]  =  Generalized  mass 


and  since 


{<7}  =  [*]{»} 

where  {»>}  is  the  modal  participation  factor  with  respect  to  [4»]. 
and 


{x}  =  [*]{<?} 

the  mode  shapes  of  the  unloaded  component  are  given  as  the  product 

[*][»] 

Thus  in  effect,  the  kinetic  energy  of  the  added  mass  has  been  subtracted  from  the  mass  loaded 
component  modal  data  If  the  added  mass  possesses  flexibility  as  well,  then  the  unloaded  generalized 
stiffness  is  given  as 

[*r  [V[*]I»[a a#  ][*i][*] 

where  [A/C  ]  is  the  stiffness  of  added  mass  systems  referred  to  attachment  DOFs. 


APPENDIX  B:  SOFTWARE  DEVELOPMENT 


The  primary  computer  software  developed  in  the  component  dynamic  synthesis  (Section  5)  task  of 
this  contract  is  a  computer  program  COMSYN  for  the  synthesis  and  analysis  of  dynamics  of  built-up 
systems  from  component  subsystem  test  or  analysis  data  The  objective  of  this  appendix  is  to  describe 
the  design,  operation,  and  usage  of  the  computer  program.  The  mathematical  formulation  on  which 
COMSYN  is  based  is  contained  in  Section  5.2  of  this  report.  The  overall  organization  of  the  program 
is  discussed,  together  with  procedures  for  operation  and  selected  programming  information  required 
to  operate  and  modify  the  program.  A  description  of  the  organization  of  the  program  is  first  given, 
followed  by  details  concerning  the  usage  of  the  program.  Input  instructions,  sample  data  and 
procedure  to  execute  the  program  is  given  in  Appendices  C  and  D  respectively. 

All  COMSYN  coding  is  done  in  ANSI-77  FORTRAN,  and  the  present  version  is  configured  to  run 
on  VAX/VMS  machines.  The  program  interfaces  to  component  analysis  and/or  modal  testing 
software  through  data  files  described  in  Appendix  C. 

B.l  Program  Organization 

A  block  diagram  of  the  COMSYN  system  is  shown  in  Figure  B-l.  The  program  operation  is  divided 
into  three  phases;  component  dynamics  data  processing;  superelement  reduction;  and  system 
synthesis  solution  and  back  substitution. 

In  the  component  data  processing  phase,  the  following  operations  are  performed:  component  and 
system  sizing  parameters,  topology,  and  other  control  information  is  read  from  user  input  files;  and 
component  coordinate  rotations  and  modal  truncations  are  performed  and  the  processed  data  is 
written  to  sequential  access  internal  files. 

In  the  reduction  phase,  the  component  dynamics  models  are  transformed  to  superelement 
coordinates  and  the  reduced  dynamic  matrices  are  written  to  internal  files  for  assembly. 

In  the  synthesis  phase,  the  following  operations  are  performed:  assembly  of  superelement  matrices; 
system  eigensolution  and  modal  transformations;  and  back  substitution  to  recover  component 
displacements,  stored  energy,  etc. 

The  flow  of  information  between  various  modules  is  through  unformatted  low  speed  sequential 
system  files  indicated  in  Figure  B-l. 

COMSYN  Call  Tree: 

The  hierarchy  of  COMSYN  subroutines  is  given  in  the  following  tree.  Here  the  order  in  which  the 
subroutines  are  called  is  given  by  the  row  number.  The  level  at  which  the  subroutines  are  called  is 
given  by  the  column  number.  All  utility  subroutines  are  grouped  under  the  fourth  level. 
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COMPiO 


RDPRM 

GETFL  RDUVF 

RDiO  RDMAT 

WRMAT 
ROTMAT 
SINTLZM 
SCALEM 
RMVSR 
RMVSM 
COLX 
ROWX 

3  SEiO  RDPRM  TRIPLE 

MTRAN 

MULT 

ADDSM 

TRANS 

MULT 

COPY 

SELECT 

ZNEW 

INVERT 

ERRMSG 

PRINTP 

PRINTSM 

INVERT 

SYMINV 

INITCM 

PRIMAT 

4  ASM  ASMBL  PRIVEC 

SKM 

5  SOLVE  EIGEN  REDUC 

TQL2 

TRED2 

REBAK 

ENRG 

6  RECOVR  DISPiO 


The  overall  logical  flow  of  the  program  and  the  function  of  the  major  program  segments  is  given 
below.  Third  and  fourth  level  subroutines  perform  rudimentary  operations  such  as  read,  write,  matrix 
manipulations,  etc.  and  are  documented  in  the  program  listings.  The  program  uses  VAX  INCLUDE 
utility  to  maintain  consistency  between  labeled  common  blocks  used  to  pass  data  between  different 
modules  of  the  program. 

CMS 

CMS  is  the  main  driver  program.  It  reads  problem  definition  and  system  parameters  and  directs  calls 
to  component  processing,  superelement  reduction,  synthesis,  and  backsubstitution  segments  of  the 


program. 


COMPiO 

Processes  component  dynamics  matrices  of  component  type  i  (i=l:  finite  element  component;  i=2: 
free  interface  modal  component  with  residual  flexibility  attachment  modes;  i=3:  free  interface 
component  with/without  interface  loading;  i=4:  constrained  interface  component).  The  subroutine 
RDPRM  reads  component  sizing  and  topological  information  from  the  system  input  file.  Subroutine 
GETFL  assigns  logical  units  from  which  component  dynamic  data  is  read.  Subroutine  RDiO 
performs  the  component  dynamic  data  manipulations  such  as  modal  truncation,  coordinates  rotation, 
matrix  partitioning  etc. 

SEiO 

Reduces  component  type  i  dynamic  matrices  to  superelement  coordinate.  The  reduction  procedures 
use  the  fourth  level  subroutines. 

ASSIGN 

Assigns  system  equation  numbers  to  superelement  degrees  of  freedom  for  assembly. 

ASM 

Forms  system  dynamic  matrices  from  superelement  matrices. 

SOLVE 

Solves  system  eigenproblem.  Subroutine  EIGEN  calls  the  EISPAC  routines  REDUC,  TRED2, 
TQL2,  and  REBAK  for  eigencalculation.  Subroutine  MTRAN  computes  system  generalized  mass, 
stiffness,  and  damping  matrices. 

RECOVR 

Back  transforms  superelement  partition  of  system  displacement  vector  to  component  generalized  as 
well  as  physical  coordinates.  The  component  stored  energy  is  computed  in  subroutine  ENRG. 

The  subroutine  RDUVF  provides  an  interface  of  the  COMSYN  with  the  SDRC  modal  software 
generated  component  universal  data  file.  Component  dynamic  data  generated  by  any  other  means 
such  as  NASTRAN  or  ANSYS  may  be  supplied  in  the  Universal  File  format.  The  program  also 
accepts  the  component  dynamic  data  in  an  alternate  simpler  format  called  Neutral  file  format 
described  in  Appendix  C.2. 

Solution  Method: 

COMSYN  uses  the  similarity  transformation  based  QL  method  of  eigensolution.  The  system  mass 
and  stiffness  matrices  are  symmetric  and  fully  populated  with  possibly  singular  stiffness  matrix  and  a 
positive  definite  mass  matrix.  The  eigenproblem  is  solved  as  follows.  Consider  the  system  equations 

(B-l) 

Since  M  is  symmetric  and  positive  definite,  its  Choleski  factorization 

\M\-\L]\Lf 


B-4 


exits  enabling  the  Equation  (B-l)  to  be  transformed  to  an  equivalent  standard  symmetric 
eigenproblem 

[k]{Y}  =  \  {y}  (B-i) 

with 

[K]  =  [L  ]'T  and  {Y}  =  [Lf{^} 

The  eigenvalues  of  the  Equation  (B-2)  are  obtained  by  first  reducing  the  symmetric  matrix  [K]  into  a 
tridiadonal  form  using  the  orthogonal  similarity  transformation.  Starting  with  j=N  the  elements  in 
the  j*  row  to  the  left  of  the  diagonal  are  first  scaled.  The  sum  of  squares  of  these  scaled  elements  is 
next  formed.  Then  a  vector  {U}  and  a  scalar  H,  where 

H={U}t  {U}/ 2, 

define  an  operator, 

[P}  =  [I]-{U}{U}t/H, 

which  is  orthogonal  and  symmetric  and  for  which  the  similarity  transformation  [P]T[K][7>]  eliminates 
the  elements  in  the  row  of  [K]  to  the  left  of  the  subdiagonal  and  the  symmetric  elements  in  the 
column.  The  eigenvalues  and  eigenvectors  of  the  symmetric  tridiagonal  matrix  [X]*  are  found  using 
QL  transformation,  where  [Q]  is  an  orthogonal  matrix  and  [L]  is  a  lower  triangular  matrix  such  that 

[k]‘=[Q][L] 

A  sequence  of  symmetric  tridiagonal  matrices  is  formed  which  converges  to  a  diagonal  matrix.  The 
rate  of  convergence  of  the  sequence  is  improved  by  shifting  the  origin  at  each  iteration.  The 
similarity  transformations  are  the  orthogonal  eigenvectors  of  the  matrix  [K],  Backward  substitution 
in 

[Lf  {Z}={Y} 

then  leads  to  the  eigenvalues  of  the  matrix  system  of  Equation  (B-l).  The  subroutines  required  for 
the  above  operations  are  extracted  from  EISPACK  ^75l 

Program  Compatibility: 

COMSYN  is  written  for  use  with  the  VAX  standard  operating  system  to  run  on  the  Wright- 
Patterson  Air  Force  Base  VAX-11/780.  All  programming  is  done  using  VAX  FORTRAN  77, 
machine  language  is  not  used.  All  I/O  is  accomplished  with  VAX  FORTRAN  77  statements  and 
does  not  use  as  an  integral  part  of  the  operation  program  any  on-line  card  reader,  card  punch,  or 
printer.  The  standard  system  logical  unit  5  (SYSS1NPUT)  for  card  reading  and  unit  6 
(SYSSOUTPUT)  for  printing  is  used  to  facilitate  compatibility  with  other  computers.  Fortran  logical 
unit  5  is  equated  to  system  input  file  NSIN  and  Fortran  logical  unit  6  is  equated  to  output  file 
NOUT.  The  disk  file  organization  type  is  sequential  access  available  with  VAX  Record  Manager. 
The  program  uses  standard  VAX  issued  basic  function  routines  such  as  SQRT,  SINE,  etc. 


B-5 


Tape  Capability 

The  program  accepts  unlabeled  or  ANSI  standard  labels  in  ASCII  format  recorded  at  either  800, 
1600,  or  6250  bpi  density.  However,  it  is  recommended  that  the  tapes  be  9  track  (1600  characters  per 
inch),  unlabelled,  with  block  size  less  than  or  equal  to  5720  characters. 

Dynamic  Storage  Allocation 

Due  to  the  wide  potential  range  of  component  and  resulting  system  sizes  and  the  desirability  of 
maximizing  efficient  usage  of  the  central  memory  by  minimizing  page  faults,  dynamic  storage 
allocation  is  important  to  the  practical  utility  of  COMSYN.  Implementation  of  dynamic  storage 
allocation  is  accomplished  by  using  the  main  program  segment  to  set  up  pointers  to  the  start  of  every 
dynamic  array.  The  basic  advantage  of  the  scheme  is  that  except  for  this  program  segment  which 
establishes  pointer  set,  the  array  dimensions  are  dummy  argument  as  opposed  to  fixed.  The  means 
for  accomplishing  the  dynamic  dimensioning  is  by  passing  the  pointer  set  as  dummy  arguments  from 
the  main  program  segment  to  lower  subroutines.  These  lower  level  subroutines  then  refer  to  each 
pointer  by  the  actual  array  name  which  that  pointer  references.  The  complete  process  involves  the 
following  steps:  read  in  data  describing  problem  size;  use  this  data  to  establish  pointers  to  all 
dynamic  arrays,  compute  total  size  of  dynamic  memory;  increment  pointer  set  to  point  to  locations 
within  new  memory  block;  pass  pointer  set  to  lower  level  subroutines.  Sufficiently  large  memory 
block  must  be  included  in  the  blank  common  to  accommodate  the  total  dynamic  memory  required. 
Oversizing  this  dynamic  block  does  not  adversely  impact  the  VMS  system  in  any  way,  since  pages 
which  are  never  accessed  do  not  impact  the  program  in  any  way. 

B.2  Data  Input  Instructions 

The  following  describes  the  preparation  of  input  for  the  COMSYN  program.  The  program  requires 
two  sets  of  data: 

Set  1:  Sizing  and  Topological  Data 
Set  2:  Component  Dynamics  Data 

The  first  set  of  input  contains  the  definition  of  the  problem,  sizing,  and  other  control  information  for 
the  system  as  well  as  the  components  and  component  connectivity  definition.  The  second  set  of  input 
data  contains  the  component  dynamics  data  as  obtained  in  independent  component  tests  and/or 
analyses.  Because  each  component  may  be  analyzed  and/or  tested  independent  of  other  components, 
a  separate  data  file  for  each  component  is  required.  The  specification  formats  for  these  data  sets  are 
described  in  Appendix  C.  Specific  order  is  required  to  data  input. 

B.3  Program  Usage 

The  usage  of  the  program  involves  preparation  of  the  program  input  data  and  execution  of  the 
program.  Component  dynamics  data  in  most  cases  is  obtained  from  separate  test  and/or  analysis 
software.  The  program  accepts  this  data  in  the  industry  standard  universal  file  format  and  also  in  the 
neutral  file  format  as  described  in  Appendix  C.  The  program  is  executed  by  entering  @CMS  PI  P2 
P3  which  accesses  the  executable  binary  and  all  the  needed  component  data  files;  the  parameters 
P1,P2  and  P3  are  respectively  the  name  of  the  executable  binary,  system  input  and  system  output  files. 
Appendix  D  describes  the  procedure  to  use  the  program  through  an  example  problem. 


APPENDIX  C:  COMSYN  INPUT  INSTRUCTIONS 


COMSYN  data  consists  of  two  parts  :  (a)  sizing  and  topological  data,  and  (b)  component  dynamics 
data.  Contents  of  the  two  data  types  are  described  below. 

C.1  SIZING  AND  TOPOLOGICAL  DATA 


Record 

Column 

Format 

Data  Item 

Description  Notes 

1-3 

1-80 

20A4 

TITLE1 

System  Header  Data 

4 

1-5 

15 

NG 

No.  of  System  Nodes 

(1) 

6-10 

15 

NDFN 

No.  of  Degrees  of  Freedom 
at  each  Node 

(2) 

5 

1-5 

15 

NC 

No.  of  Components  in  the 
System 

6 

1-5 

15 

IRES 

Solution  Type  Flag 

(3) 

7 

1-5 

15 

METHOD 

Method  of  Eigensolution 

(4) 

6-10 

15 

NVEC 

No.  of  System  Modes  to 
output 

8 

1-10 

A10 

KEYWORD 

"COMPONENT  " 

(5) 

11-15 

15 

KID 

Component  Identification 
Number 

9 

1-5 

A5 

TYPE 

"TYPE  " 

1-5 

15 

KTYPE 

Component  Dynamics  Data 
Type 

(6) 

10 

1-5 

15 

KN 

No.  of  Component  Grid 
Points 

(7) 

11 

1-5 

15 

KB 

No.  of  Component 

Boundary  Nodes 

(8) 

12 

1-80 

1615 

KBLI5T ( I ) , 
1=1, KB 

List  of  Boundary  Nodes 

(9) 

13 

1-80 

1615 

KSLIST ( I ) , 

Component  to  System 

1=1, KB 

Correspondence  Table 

(10) 

! 
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Record 

Column 

Format 

Data  Item 

Description 

Notes 

14 

1-5 

15 

KM 

No.  of  Modes  Supplied 
in  input  data 

(11) 

15 

1-5 

15 

KK 

No.  of  Modes  to  be 

Used  in  Synthesis 
(KK.LE.KM) 


16 

1-80 

1615 

KKLIST ( I) , 

List  of  Modes  to  be 

(12) 

1=1, KK 

Used 

17 

1-60 

13 

KFLAG(I) , 

Flags  for  Specifying  Type 

1=1,20 

and  Data  Storage  Mode 

(13) 

61-70 

E10.3 

EPS 

Residual  Mass  Parameter  if 

KFLAG ( 8 ) =0 

18 

1-30 

3F10.0 

ANGL(I) , 

Three  Sequential  Angles 

1=1,3 

(Roll,  Pitch,  Yaw) 
Defining  Orientation  of 
of  Component  Axes  with 
Respect  to  System  Axes 
(see  Figure  5) 

19 

1-13 

A13 

KEYWORD 

"DYNAMICS  DATA" 

14-15 

12 

IOPT 

Dynamic  Data  File  Option 

(14) 

16-25 

A10 

FNAME 

File  Name  for  Dynamic 

Data  File 

(15) 

26-30 

15 

IDEN 

Tape  Density 

(16) 

31-35 

15 

IRDWR 

Read/Write  Flag 

(16) 

3-45 

A10 

MGTPFN 

Magnetic  Tape  File  Name 

(16) 

46-50 

15 

IBLF 

Blocking  Factor 

Record  numbers  12, 13,  and  16  may  run  into  several  cards  each. 


Repeat  record  numbers  8  through  19  for  each  component. 

Data  for  each  component  must  be  in  sequence;  that  is,  component  i  must  precede  component  i+1  and 
for  each  component  the  data  must  be  ordered  as  per  instructions  above. 

NOTES: 

1.  Number  of  system  nodes  includes  all  the  connection  nodes  plus  the  number  of  nodes  which  are 
to  be  restrained  in  system  analysis. 

2.  Valid  input  for  the  number  of  degrees  of  freedom  is  one  of: 

1 :  One  translational  degree  of  freedom  per  node 

2  :  Two  translational  degrees  of  freedom  per  node 

3  :  Three  translational  degrees  of  freedom  per  node 

6  :  Three  translational  and  three  rotational  degrees  of  freedom  per  node, 
referred  to  three  mutually  orthogonal  axes 


3.  IRES=1  for  eigensolution. 

4.  Method  of  eigensolution  used  is  QL  method. 

5.  Repeat  record  numbers  8  through  19  for  each  component.  Data  for  each  component  must  be 
in  sequence;  that  is,  component  i  must  precede  component  i+1  and  for  each  component  the 
data  must  be  ordered  as  per  instructions  above. 

6.  Component  data  types  allowed  are: 

10  :  Finite  element  stiffness  mass  and/or  damping  data 

20  :  Free  interface  normal  modes  plus  residual  flexibility  data 

30  :  Free  interface  normal  mode  data 

40  :  Constrained  interface  normal  mode  data 

7.  Component  grid  points  must  include  the  component  partition  of  system  nodes;  that  is,  all  of  its 
connection  points  and  any  other  node  which  is  to  be  carried  into  system  synthesis.  Each  node 
has  NDFN  degrees  of  freedom  (see  Note  2). 

8.  Boundary  nodes  include  nodes  connected  to  adjacent  components  plus  nodes  that  will  be 
restrained  in  system  analysis.  The  node  labels  are  those  used  in  component  grid  point 
numbering. 

9.  Record  type  12  is  to  be  repeated  in-place  as  necessary,  if  the  number  of  component  boundary 
nodes  (KB)  is  greater  than  16. 

10.  Connectivity  of  the  component  is  specified  by  the  component  to  system  correspondence  table 
that  lists  the  system  node  number  for  each  of  the  nodes  specified  in  KBLIST. 

Record  type  13  is  to  be  repeated  in-place  as  necessary,  if  the  number  of  component  boundary 
nodes  (KB)  is  greater  than  16. 

11.  Number  of  nodes  supplied  on  the  component  data  file.  For  type  10,  component  KM=0.  For 
Type  20  and  Type  30,  the  number  of  nodes  include  rigid  body  nodes  (if  any)  plus  elastic  normal 
nodes. 

12.  Record  16  is  to  be  repeated  in-place  as  necessary,  if  the  number  of  modes  to  be  used  (KK)  is 
greater  than  16. 

13.  Flags  refer  to  the  type  and  storage  mode  for  component  dynamic  data,  in  the  following  manner: 


ENTITY  KFLAG  VALUE  MEANING 

Stiffness  Matrix  Kx  1  0  Absent 

1  Vector  of  diagonals 

2  Full  matrix,  row-wise 

3  Full  matrix,  symmetric 
upper  triangular  (row) 

4  Full  matrix,  symmetric 
lower  triangular  (row) 

Mass  Matrix  Mx  20 

1 

2  (same  as  for  KFLAG (1)) 

3 

4 

Damping  Matrix  Cx  3  0 

1 


*>  v  wyffvwyvTro?  vr^  v  v  v  ’aw  w 


Modal  Stiffness  K, 


Modal  Mass  Mg  5 


Modal  Damping  Cq 


Residual  Flexibility  7 
#gB 


Residual  Mass  HB 


Residual  Damping  CB  9 


Constraint  Mode 
Stiffness  K^c 


10 


2 

3 

4 


0 

1 

2 

3 

4 


0 

1 

2 

3 

4 


0 

1 

2 

3 

4 


0 

1 


0 

1 


2 

3 


0 

1 


2 

3 


0 

1 

2 


C-5 


(same  as  for  KFLAG(l)) 


(same  as  for  KFLAG ( 1 ) ) 


(same  as  for  KFLAG (1)) 


(same  as  for  KFIAG(l)) 


Absent 
GrB  partition  only. 

n  T _  *  _ rt  *1  *  _ i.  . 


Size  KBFxKBF,  listed  row¬ 
wise. 


’9: 


KNFxKBF  size 


listed  row-wise. 


Absent 

HpB  partition  only  KBFxKBF 
size,  listed  row-wise. 

Hp  partition  KNFxKBF  size, 
listed  row-wise 
Compute  Gb.Mx.Gb  (full  mass 
matrix  in  physical 
coordinates  must  be 
supplied  as  per  KFLAG (2)) 


Listed 


Absent 

CBB  partition  only, 
row-wise 
Full  CB  partition.  Listed 
row-wise . 

Compute  GB.CX.GB  (Full 
damping  must  be  supplied 
as  per  KFLAG (3) ) 


Absent 

Full  matrix,  row-wise 
Full  matrix,  symmetric, 
upper  triangular  (row) 
Full  matrix,  symmetric. 


_  >  ,V 


.V.  -V  • 


<r.  s.  v. 


1 


SI 


* 


1 


A 


1 


3? 


3 


3 


■■1 


'j-, 


1 


?! 


lower  triangular  (row) 


Constraint  Mode 

11 

0 

Absent 

Mass  Mq£ 

1 

Full  matrix,  row-wise 

2 

Full  matrix,  symmetric, 
upper  triangular  (row) 

3 

Full  matrix,  symmetric, 
lower  triangular  (row) 

Normal  Mode-Constraint 

12 

0 

Absent 

Mode  Coupling  Mass  Mnc 

1 

Full  matrix  row-wise 

Hysteretic  Damping  Hx 

13 

0 

i 

X 

2 

3 

(same  as  for  KFLAG(l)) 

Modal  Inverse  Method 

14 

0 

Generalized  inverse  using 
Gaussian  elimination 

1 

Generalized  inverse  using 
singular  value 
decomposition 

KFLAG(15)  through  KFLAG(20)  are  used  for  controlling  level  of  printed  output  in  various 
modules.  Thus  a  unit  value  of  the  respective  flag  outputs  data  in  subroutines  COMPiO,  SEiO, 
ASSIGN,  ASMBL,  SOLVE  and  RECOVR. 

14.  Dynamic  data  may  be  read  in  two  ways:  IOFT=I  Read  dynamic  data  from  user  supplied  file  in 
Neutral  file  format  defined  in  Appendix  C.2.  IOPT=2  Read  dynamic  data  from  universal  file 
format  defined  in  Appendix  C.3. 

15.  On  VAX  the  program  automatically  accesses  the  file  named  FNAME.DAT. 

16.  The  tape  density,  read/write  flags,  magnetic  tape  file  name,  and  the  blocking  factor  refer  to  the 
SDRC  software  created  universal  file. 


C.2  COMPONENT  DYNAMICS  DATA:  NEUTRAL  FILE  FORMAT 


The  specification  format  for  component  dynamic  data  differs  with  the  type  of  component  dynamic 
model  as  described  below. 


Finite  Element  Component:  Type  10 

Record  Column  Format  Data  Item  Description 


Notes 


1 

1-80 

A80 

TITLE2 

Component  Header  Data 

2 

1-80 

8E10.3 

*x 

Stiffness  Matrix  in 
Physical  Coordinates 

(1) 

3 

1-80 

8E10.3 

Mass  Matrix  in  Physical 
Coordinates 

(1) 

4 

1-80 

8E10.3 

cx 

Viscous  Damping  Matrix 
in  Physical  Coordinates 

(1) 

Free 

Type 

Interface 

20 

Modal  Component 

with  Residual  Flexibility: 

Record  Column 

Format 

Data  Item  Description 

Notes 

1 

1-80 

20A4 

TITLE2 

Component  Header  Data 

2 

1-5 

15 

MOD 

Mode  Number 

(15) 

3 

1-10 

11-20 

21-30 

E10.3 

E10.3 

E10.3 

gK 

°M 

GC 

Modal  Stiffness 

Modal  Mass 

Modal  Damping 

(2,15) 

4 

1-4 

6-65 

14 

6E10.3 

NUM 
U(I)  , 
1=1,6 

Node  Number 

Node  Point  Modal 
Coefficients 

(3) 

(14,15) 

5 

1-80 

8E10.3 

gB  or 

gbb 

Residual  Flexibility 

(4) 

6 

1-80 

8E10.3 

Hb  or 
hbb 

Residual  Inertia 

(5) 

7 

1-80 

8E10. 3 

MX 

Mass  Matrix 

(6) 

8 

1-80 

8E10.3 

CB  or 
CBB 

Residual  Damping 

(7) 

9 

1-80 

8E10 . 3 

CX 

Damping  Matrix 

(8) 

Free  Interface  Modal  Component  With  or  Without  Interface  Loading: 
Type  30 

Record  Column  Format  Data  Item  Description  Notes 


1 

1-80 

20A4 

TITLE2 

Component  Header  Data 

2 

1-5 

15 

MOD 

Mode  Number 

(15) 

3 

1-10 

11-20 

21-30 

E10.3 

E10.3 

E10.3 

gK 

gM 

GC 

Modal  Stiffness 

Modal  Mass 

Modal  Damping 

(9) 

(15) 

4 

1-4 

6-65 

14 

6E10.3 

NUM 

U(I)  , 
1=1,6 

Mode  Number 

Node  point  Modal 
coefficients 

(10) 

(14) 

(15) 

5 

1-80 

8E10. 3 

Kq 

Generalized  Stiffness 
Matrix  (KM*KM  size) 

(9) 

6 

1-80 

8E10 . 3 

Mq 

Generalized  Mass  Matrix 
(KM  *  KM  size) 

(9) 

7 

1-80 

8E10.3 

CQ 

Generalized  Damping 

Matrix  (KM  *  KM  size) 

(9) 

Constrained  Interface  Modal  Component:  Type  40 


Record  Column  Format  Data  Item  Description  Notes 


1 

1-80 

20A4 

TITEL2 

Component  Header  Data 

2 

1-5 

15 

MOD 

Constrained  Normal  Mode 
Number 

(11) 

(15) 

3 

1-10 

E10.3 

gK 

Constrained  Normal  Mode 
Stiffness 

(2) 

(15) 

4 

11-20 

E10.3 

°M 

Constrained  Normal  Mode 
Mass 

(2) 

(14,15) 

5 

21-30 

E10. 3 

GC 

Constrained  Normal  Mode 
Damping 

(2) 

(16) 

6 

1-4 

6-65 

14 

6E10.3 

NUM 

U(I)  , 
1=1,6 

Node  Number 

Node  Point  Moda 
Coefficients 

(10) 

(16) 

7 

1-5 

15 

MOD 

Constraint  Mode  Number 

(12) 

8 

1-4 

6-65 

14 

6E10.3 

NUM 

U(I), 

1=1,6 

Node  Number 

Node  point  Modal 
Coefficient 

(10) 

9 

1-80 

6E10.3 

KCC 

Constraint  Mode  Stiffness 
Matrix  (KBF  *  KBF  size) 

(13) 

10 

1-80 

6E10.3 

Mcc 

Constraint  Mode  Mass 
Matrix  (KBF  *  KBF  size) 

(13) 

11 

1-80 

8E10.3 

%c 

Normal  Mode-Constraint  (13) 

Mode  Inertia  Coupling  Matrix 
(KM  *  KBF  size) 

NOTES: 

1.  The  specification  of  the  matrix  coefficients  A(IJ)  should  be  consistent  with  the  corresponding 
flags  specified  in  record  17  of  sizing  data  described  in  Appendix  Cl.  Each  record  may  run  into 
several  data  cards.  The  matrices  should  conform  with  the  free  interface  boundary  conditions  of 
the  component. 

2.  Physical  units  of  the  generalized  mass,  damping,  and  stiffness  are  at  the  user’s  discretion;  any 
consistent  units  may  be  selected. 

3.  The  node  point  modal  coefficients  should  be  ordered  as  u,v,w,rx,ry,rz,  where  u,  v,  w,  and  rx,  ry, 
rz  are  respectively  the  three  translational  displacements  and  three  rotations  with  respect  to 
Cartesian  axes  x,  y,  z. 

4.  Residual  flexibility  matrix  coefficients  are  input  row-wise  as  ((GBB(IJ)J=1,KBF),I=1,KBF)  if 
KFLAG(7)=1  and  as  ((GB(U)J-1,KBF),I=1,KNF)  if  KFLAG(7)=2.  KBF  and  KNF  are 
respectively  the  number  of  boundary  DOFs  and  total  nodal  DOFs  of  the  component. 


5.  Residual  inertia  matrix  coefficients  are  input  row-wise  as  (HBB(U),J=1,KBF),I=1,KBF)  if 
KFLAG(8)=1  or  ((HB(IJ)J=1,KBF),I=1,KNF)  if  KFLAG(8)=2.  Omit  this  record  set  if 
KFLAG(8)=0  or  3. 

6.  If  KFLAG(8)=3,  physical  mass  matrix  and  the  full  G  partition  of  the  residual  flexibility  matrix  is 
required  ((KFLAG(7)  must  be  2)  to  compute  residual  mass  matrix.  Input  mass  matrix  as  per 
KFLAG(2).  Omit  this  record  if  KFLAG(8).NE.3. 

7.  Residual  damping  matrix  coefficients  are  input  row-wise  as  ((CBB(IJ)J=1,KBF),1=1,KBF)  if 
KFLAG(9)=1  or  ((CBB(U)J=1,KBF),I=1,KNF)  if  KFLAG(9)=2.  Omit  this  record  set  if 
KFLAG(9)=0  or  3. 

8.  If  KFLAG(9)=3,  physical  damping  matrix  and  the  full  G  partition  of  the  residual  flexibility 
matrix  is  required  ((KFLAG(7)  must  be  2)  to  compute  residual  damping  matrix.  Input 
damping  matrix  as  per  KFLAG(3).  Omit  this  record  if  KFLAG(9).NE.3. 

9.  Physical  units  of  the  generalized  mass,  damping,  and  stiffness  are  at  the  user’s  discretion;  any 
consistent  units  may  be  selected. 

When  nonorthogonal  modes  are  used  one  or  more  of  the  generalized  matrices  may  be  fully 
populated.  Full  matrices  are  specified  in  records  4  through  6  as  per  KFLAG(4),  KFLAG(5), 
and  KFLAG(6). 

10.  The  node  point  modal  coefficients  should  be  ordered  as  u,  v,  w,  rx,  ry,  rz,  where  u,  v,  w,  and  rx, 
ry,  rz  are  respectively  the  three  translational  displacements  and  three  rotations  with  respect  to 
cartesian  axes  x,  y,  z. 

Record  type  4  must  be  repeated  as  many  times  as  there  are  nodes  in  the  component;  that  is, 
KN  times,  where  variable  KN  has  been  set  in  record  type  10  in  the  Sizing  and  Topological  Input 
Data  File. 

11.  Constrained  normal  modes  are  normal  modes  of  the  component  obtained  with  connection 
interface  degrees  held  fixed.  The  input  formats  for  the  constrained  normal  mode  data  is  similar 
to  that  of  the  normal  mode  data. 

12.  The  constraint  modes  are  static  deflection  shapes  as  defined  in  Equations  183  through  186.  The 
component  and  nodal  degrees  of  freedom  arrangement  and  the  mode  shape  data  input  is  same 
as  that  of  the  normal  mode  shape  data.  Record  types  5  through  6  must  be  repeated  for  each 
constraint  mode.  The  number  of  constraint  modes  equals  the  number  of  boundary  degrees  of 
freedom,  NBDG.  This  can  be  calculated  by  multiplying  the  number  of  boundary  nodes,  KB, 
times  the  number  of  degrees  of  freedom  per  node  NDFN.  That  is,  NBDG  =  (KB*NDFN), 
where  KB  and  NDFN  are  the  values  specified  in  the  Sizing  and  Topological  Data  Input  File  on 
records  1 1  and  4,  respectively. 

13.  Stiffness,  mass  and  damping  matrices  associated  with  the  constraint  modes  are  defined  in 
Equation  200.  Note  the  matrix  partitions  required. 

14.  Record  Number  4  should  be  repeated  as  many  times  as  there  are  number  of  nodes  in  the 
component. 

15.  Record  2  through  4  should  be  repeated  for  each  mode. 

16.  Records  5  and  6  should  be  repeated  for  each  constraint  mode. 

C.3  COMPONENT  DYNAMICS  DATA:  UNIVERSAL  FILE  FORMAT 

The  universal  file  format  is  a  widely  accepted  form  for  interfacing  data  between  different  software. 

Details  of  the  universal  file  are  given  in  Reference  76.  Subroutine  RDUVF  translates  the  universal 
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file  contents  and  creates  a  neutral  file  for  the  component.  The  data  requirement  for  a  certain  type  of 
component  is  the  same  as  those  defined  in  Appendix  C.2;  only  the  specification  format  differs.  The 
following  describes  the  universal  file  formats  for  all  the  dynamic  data  types  needed  by  COMSYN. 

The  general  characteristics  of  a  universal  file  data  set  are  as  follows: 

•  The  first  record  of  the  data  contains  "-T*  right  justified  in  columns  1  through  6.  Columns  7 
through  80  of  this  physical  record  is  blank. 

•  The  second  record  of  the  data  set  contains  the  data  TYPE  number,  numeric  range  1  through 
32767,  right  justified  in  columns  1  through  6.  The  remaining  columns,  7  through  80,  are  blank. 

•  The  last  record  of  the  data  set  contains  "-1"  right  justified  in  columns  1  through  6.  Columns  7 
through  80  of  this  record  are  blanks. 

The  following  lists  the  universal  data  sets  for  the  four  component  types  and  the  specifics  of  each  data 
set. 

Universal  Data  Set  for  Physical  Coordinate  Component  (TYPE  10) 

Dataset  Type  Description 

241  Component  Header  Data 

250  Stiffness  Matrix  (IMAT=9) 

250  Mass  Matrix  (IMAT=6) 

250  viscous  Damping  Matrix  (IMAT=7) 


Universal  Dataset  for  Residual  Flexibility  Component  (TYPE  20) 

Dataset  Type  Description 

241  Component  Header  Data 

55  Eigenvector  and  Modal  Parameters 

250  Residual  Flexibility  Matrix 

250  Residual  Mass  Matrix  (if  KFLAG(8)=l  or  2) 

250  Total  Mass  Matrix  (if  KFLAG(8)=3) 

250  Residual  Damping  Matrix  (if  KFLAG(9)=1  or  2) 

250  Total  Damping  Matrix  (if  KFLAG(9)=3) 


Universal  Datas  it  for  General  Coordinate  Component  (TYPE  30) 

Dataset  Type  Description 

241  Component  Header  Data 

55  Mode  Shape  Vector  and  Modal  Parameters 

250  Generalized  Stiffness  Matrix  (if  KFLAG(4)>0) 

250  Generalized  Mass  Matrix  (if  KFLAG(5)>0) 

250  Generalized  Damping  Matrix  (if  KFLAG(6)>0) 


Universal  Dataset  for  Constrained  Interface  Component  (TYPE  40) 

Dataset  Type  Description 

241  Component  Header  Data 

55  Eigenvector  and  Modal  Parameters 

55  Constraint  Mode  Vector 

250  Constraint  Mode  Stiffener  Matrix  (K  partition) 

250  Constraint  Mode  Mass  Matrix  (M  partition) 

250  Normal  Mode-Constraint  Mode  Mass  Coupling 

(M  partition) 


Dataset  Type  55 

Record  Column  Format  Data  Item  Description 


1-5 

1-80 

80A1 

TITLE 

6 

1-10 

110 

1 

11-20 

110 

1 

2 

21-30 

110 

3 

6 

31-40 

110 

8 

41-50 

110 

2 

51-60 

110 

NDV 

7 

1-10 

110 

1 

2 

11-20 

110 

1 

4 

21-30 

110 

MOD 

31-40 

110 

MOD 

8 

1-13 

E13. 

5  w 

14-26 

El 3  . 

5  °M 

27-39 

E13. 

5  Gc 

9 

1-10 

110 

NOD 

10 

1-78 

6E13 

.5  U(I)  , 

1=1, NDV 


Identifying  Title  (5  lines) 

Structural  Model  Flag 
Constraint  Mode 
Normal  Mode  Shape 
3  DOF  (Translation)  per  Node 
6  DOF  (Translations  + 
Rotations)  per  node 
Flag  for  Displacement  Vector 
Real  Coefficients 
No.  of  Data  Values  per  Node 
(=data  in  field  2  above) 

No.  of  integer  data  values  for 
constant  mode 

No.  of  integer  data  values  for 
normal  mode 

No.  of  real  data  values  for 
constant  mode 

No.  of  real  data  values  for 
normal  mode 

Constraint  Mode  Number 
Normal  Mode  Number 

Normal  Mode  Frequency  in  Hz 
(blank  for  constraint  mode) 
Modal  Mass  (normal  mode) 
(blank  for  constraint  mode) 
Modal  Viscous  Damping 
(blank  for  constraint  mode) 

Node  Number 

Node  Point  Modal  Coefficient 


NOTES: 


1.  Records  9  and  10  are  repeated  for  each  node. 

2.  Title  lines  may  not  be  blank.  If  no  information  is  required  the  word  "NONE"  must  appear  in 
columns  1-4. 

3.  The  order  of  data  values  U(I),I=1,NDV  is 

3  DOF  Vector  X,Y,Z 
6  DOF  Vector  X,Y,Z jBX.eY.8Z 

4.  Modal  stiffness  is  obtained  computed  in  the  program  as  GK  =  4.0  ir3  w8  Gm,  where  w  is  natural 
frequency  in  Hertz. 


Dataset  Type  250 


Record 

Column 

Format 

Data  Item 

Description 

1 

1-10 

110 

IMAT 

Matrix  Identification 

2 

1-10 

110 

MDTYPE 

Matrix  Data  Type  (=2  for  real) 

11-20 

110 

MFORM 

Matrix  Form  (=3  General 
Rectangular 

21-30 

110 

NROWS 

Number  of  Rows 

31-40 

110 

NCOLS 

Number  of  Columns 

41-50 

110 

MKEY 

Storage  Key  (=1  Row,  =2  Column) 

3 

1-10 

110 

ISR 

Starting  Row  for  Submatrix 

11-20 

110 

ISC 

Starting  Column  for  Submatrix 

21-30 

110 

NR 

No.  of  Rows  in  Submatrix 

31-40 

110 

NC 

No.  of  Columns  in  Submatrix 

41-50 

110 

MFORMS 

Submatrix  Form  (=3  General 
Rectangular,  =5  Diagonal) 

51-60 

110 

MKEYS 

Submatrix  Storage  Key  (=1  Row, 
=2  Column) 

4 

1-80 

4E10 . 

12  A ( I , J) 

Real  Matrix  Data 

•  Repeat  record  4  as  necessary  to  fulfill  requirements  of  record  3. 

•  Records  3  and  4.  as  a  group,  are  repeated  as  necessary  to  define  all  non-zero  submatrices. 


APPENDIX  D:  EXAMPLE  JCL  AND  INPUT  DATA 


The  following  describes  a  sample  VAX/VMS  command  procedure  to  execute  program  COMSYN. 
The  example  shown  in  Figure  D-l  uses  input  data  files  which  are  all  resident  on  disk  as  permanent 
files.  The  example  problem  uses  one  system  input  file  and  four  component  dynamic  data  files.  These 
are  declared  as  LV.DAT,  COMPl.DAT,  COMP2.DAT,  COMP3.DAT  and  COMP4.DAT, 
respectively.  The  procedure  file  and  the  input  data  file  are  listed  in  Figures  D-2  through  D-7. 

The  system  considered  is  that  of  the  longitudinal  motion  of  a  launch  vehicle-like  structure.  It  is 
partitioned  into  four  components  as  shown  in  Figure  D-l.  In  order  to  cover  the  full  modeling 
capabilities  of  the  synthesis  program,  a  different  type  of  dynamic  model  is  used  to  represent  the 
component  dynamics.  Thus  the  four  components  are  modeled  respectively  as  free  interface  normal 
modes  plus  residual  flexibility;  free  interface  modes  alone;  finite  element  matrices  in  physical 
coordinates;  and  constrained  interface  modes.  The  system  input  and  component  data  file  are  as 
follows.  The  component  dynamic  data  files  are  input  to  the  program  in  universal  file  format.  In 
generating  the  component  data,  the  connection  interfaces  are  left  unrestrained.  Only  those 
component  nodes  are  restrained  in  component  model  generation  which  correspond  to  the  restrained 
system  nodes.  Thus  the  node  1  of  the  component  1  is  restrained  because  the  built-up  system  is  also 
restrained  at  the  same  location.  In  the  following  the  background  information  needed  to  create  the 
data  is  also  given  for  the  sake  of  completeness. 

Component  1  Residual  Flexibility  Component 

The  normal  modal  and  static  residual  flexibility  properties  of  the  component  are  derived  from  the 
following  stiffness  and  mass  matrices.  Node  2  is  the  boundary  node  and  is  kept  unrestrained. 

1  -1  0  01  [1  0  0  & 

Stiffness  Matrix:  [*L  =  q  \  ,  Wlo  =  0  0  2  0 

0  0  -1  lj  [o  0  0  1 

Normal  modes  and  modal  properties  are: 

0  0  0  1 
r  xi  -.288675  .57735  .288675 
l$lAr  "  -.5  0  -.5 

-.57735  -.57735  .57735. 

r.  13397  0  0 

[K]„=  0  1.0  0 

0  0  1.86603 


\M), 


r  1.0  o  oi 
0  1.0  0 
0  0  1.0 


Residual  Flexibility  Properties: 


Modes  to  be  kept 

0  0  1 
_  -.28865  .57735 
l$1* =  -.5  0 

-.57735  -.57735 


1.0  01 
0  1.0J 


Total  Static  Flexibi'ity 


0  0  0  C 
0  111 


[G]  =  [*tf=  o  1  2  2 


0  12  3 


Static  Flexibility  of  Truncated  Mode  (Residual  Flexibility) 


0  0  0  1 
.044658  -.07735  -.089316 
-.07735  .133975  -.154701 
.089316  -.154701  .178633 


Residual  Flexibility  Attachment  Mode: 


0 

.089316 

-.154701 

.178633 


The  universal  data  set  for  the  above  component  COMP1.DAT  is  listed  in  Figure  D-3. 


Component  2:  Free  Interface  Normal  Mode  Component 


The  modal  properties  are  derived  from  the  stiffness  and  mass  matrices  in  physical  coordinates  given 
for  the  Component  1.  the  mode  shapes,  modal  stiffness,  and  modal  mass  are  given  as 


r-408  .577  .577  -.408 

.408  -.2887  -.2887  .408 
.408  -.2887  -.2887  -.408 
.408  .577  .577  .408 


0  0  0 


[KT  =  0  0.5  0  ,  [M] 


0  0  1.5 


ri.  o  oi 
],=  0  1.  0 
0  0  1. 


The  universal  data  set  for  the  above  component  COMP2.DAT  is  listed  in  Figure  D-4. 


Component  3:  Finite  Element  Component 


The  stiffness  and  mass  matrices  for  this  component  are  as  shown  above  for  the  Component  1.  The 
universal  data  set  for  this  component  is  listed  in  Figure  D-5. 


Component  4:  Constrained  Interface  Modal  Component 


The  component  is  constrained  at  the  nodes  1  and  4.  The  Node  4  is  treated  as  a  boundary  node  since 
it  is  to  be  represented  in  the  system  equation  in  physical  coordinates.  The  normal  mode  properties 
are  obtained  from  the  solution  of  the  following  eigenproblem. 


m  MS0] ]{;:}■«« 


(WO 


leading  to 


0  01 

[H-[o  u]-  1°1 

.0  OJ 

The  constraint  modes  are  derived  by  solving  the  two  statics  problems  obtained  by  (a)  constraining 
node  1  and  imposing  a  unit  displacement  to  node  2,  and  (b)  constraining  node  2  and  imposing  a  unit 
displacement  to  node  1.  The  two  constraint  nodes  are 

1.  0 
667  .333 
333  .667 
0  1. 

and  the  corresponding  stiffness  and  mass  properties  are 


[K]cc-[_'333  .333*]’  Wbc  -  [  2.1in]’  3X1(1  \-m\hc  =  [.333  ..333] 


l< 

I 


i 

* 

\ 

i 

r 

f 

s' 


Ll'kl'Ll  ktlll  A 


with  one  DOF  at  each  node.  The  dynamics  data  type,  sizing,  connectivity,  and  other  control 
information  for  each  component  is  also  shown  in  Figure  D-7. 
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System 
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COMPONENT  3 


COMPONENT  4 
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O 

4  - 


COUPLED  SYSTEM 


System  DOFs 


Figure  D-l.  Example  Problem 


$! Command  Procedure  File  CMS.COM  : 

$ ! 

$ !  PROCEDURE  TO  RUN  COMSYS  INTERACTIVELY 

$  1 

$ ! - SYSTEM  INPUT  FILE 

$! 

$  ASSIGN/USER_MODE  ' P2 '  FOR005 
$ ! 

$ ! - OUTPUT  PRINT  FILE 

$ ! 

$  ASSIGN  USERMODE  ' P3 '  FOR006 
$ ! 

$  ! - EXECUTION  STEP 

$ ! 

$  RUN  'PI' 

$  EXIT 

TO  Start  Execution  Enter: 

@CMS.COM  PI  P2  P3 
Where : 

PI  =  File  Name  Containing  COMSYN  Executable  Binary 
P2  =  File  Name  Containing  System  Input  Data 
P3  =  File  Name  for  Output  File 


Figure  D-2.  Procedure  File  for  COMSYS  Execution 


-1 

241  (Component  Header  Data) 

6 

KM01 

DYNAMIC  DATA  FOR  COMPONENT  1  -  UNIVERSAL  FILE  INPUT 
15-OCT-86 

1  6 
-1 
-1 

55  (Mode  Shape  Data  at  Nodes:  Mode  1) 

COMPONENT  1:  RESIDUAL  FLEXIBILITY  COMPONENT 

NORMAL  MODE  1 

NONE 

NONE 

NONE 

12  18  2 

2  4  11 

.0582538  1.0  0.  0. 

1 

.0  0.  0.  0. 

2 

-.288675  0.  0.  0. 

3 

- .  5  0 .  0 .  0 . 

4 

-.577350  0.  0.  0. 

-1 

-1 

55  (Mode  Shape  Data  at  Nodes  :  Mode  2) 

COMPONENT  1:  RESIDUAL  FLEXIBILITY  COMPONENT 

NORMAL  MODE  2 

NONE 

NONE 

NONE 

12  18  2 

2  4  12 

.15915  1.0  0.  0. 

1 

.0  0.  0.  0. 

2 

. 57735  0 .  0 .  0 . 

3 

.0  0.  0.  0. 

4 

-. 577350  0.  0.  0. 

-1 


Figure  D-3.  Universal  File  for  Component  1 


-1 

55  (Mode  Shape  Data  at  Nodes  :  Mode  3) 
COMPONENT  1:  RESIDUAL  FLEXIBILITY  COMPONENT 
NORMAL  MODE  3 
NONE 
NONE 
NONE 

1  2 

2  4 

.21745  1.0 

1 


.0 

.288675 

-.5 


2 

3 

4 


0. 

0. 

0. 


0. 

0. 

0. 


-.57735  0.  0. 

-1 

-1 

251  (Residual  Flexibility  DOF  Data) 
1 

1  1 


-1 


-1 

250  (Residual  Flexibility  Matrix  Data) 
149 


2  3  4  1 

11  4  1 

0.  .089316 

.  178633 
-1 


2 


0. 

0. 

0. 

0. 

0. 


1 

3 

-.154701 


Figure  D-3  (continued) 


241  (Component  Header  Data) 

6 

KM02 

DYNAMIC  DATA  FOR  COMPONENT  2  -  UNIVERSAL  FILE  INPUT 
15-OCT-86 

1  6 
-1 
-1 

55  (Mode  Shape  Data) 

COMPONENT  2:  FREE  INTERFACE  MODAL  COMPONENT 


NORMAL  MODE  1 

NONE 

NONE 

NONE 

1 

2 

1 

8 

2 

2 

4 

1 

1 

0. 

1.0 

0. 

0. 

1 

.408248 

0. 

0. 

0. 

2 

.408248 

0. 

0. 

0. 

3 

.408248 

0. 

0. 

0. 

4 

.408248 

0. 

0. 

0. 

-1 

-1 

55  (Mode  Shape 
COMPONENT  2:  FREE 
NORMAL  MODE  2 

NONE 

NONE 

NONE 

1 

Data) 

INTERFACE 

2 

COMPONENT 

1 

8 

2 

2 

4 

1 

2 

.  5 

1.0 

0. 

0. 

1 

.  57735 

0. 

0. 

0. 

2 

-  .288675 

0. 

0. 

0. 

3 

-.288675 

0. 

0. 

0. 

4 

. 57735 

0. 

0. 

0. 

-1 


Figure  D-4.  Universal  File  for  Component  2 


55  (Mode  Shape  Data) 

COMPONENT  2  FREE  INTERFACE  MODAL  DATA  (TYPE  30) 

NORMAL  MODE  3 

NONE 

NONE 

NONE 

12  18 

2  4  13 

1.5  1.0  0.  C 

1 

.57735  0.  0.  C 

2 

-.288675  0.  0.  C 

3 

-.288675  0.  0.  C 

4 

.57735  0.  0.  C 

-1 
-1 

55  (Mode  Shape  Data) 

COMPONENT  2  FREE  INTERFACE  MODAL  DATA  (TYPE  30) 

NORMAL  MODE  4 

NONE 

NONE 

NONE 

12  18 

2  4  14 

2.0  1.0  0.  ( 

1 

-.408248  0.  0.  ( 

2 

.408248  0.  0.  ( 

3 

-.408248  0.  0.  ( 

4 

.408248  0.  0.  ( 

-1 


Figure  D-4  (continued) 
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-1 

241 

\ 

6 

KM03 

DYNAMIC  DATA  FOR  COMPONENT  3  -  UNIVERSAL  FILE  INPUT 
15-OCT-86 


(Mass  Matrix) 
6 

2  3 

1  1 

1.0 


(Stiffness  Matrix) 
9 


2 

1 

1.0 

-1.0 

0.0 

0.0 


4 

4 

-1.0 

2.0 

-1.0 

0.0 


1 

3 

0.0 

-1.0 

2.0 

-1.0 


(Viscous  Damping  Matrix) 
7 


2 

1 

0.05 

-0.05 

0.0 

0.0 


4 

4 

-0.05 

0.1 

-0.05 


1 

5 

0.0 

-0.05 

0.1 

-0.05 


0.0 

0.0 

-0.05 

0.05 


Figure  D-5.  Universal  File  for  Component  3 


vxv--:-/-.-- 


-1 

241 

6 

KM04 

DYNAMIC  DATA  FOR  COMPONENT  4  -  UNIVERSAL  FILE  INPUT 
15-OCT-86 

1  6 
-1 
-1 
55 

COMPONENT  4  CONSTRAINED  INTERFACE  MODAL  COMPONENT 

NORMAL  MODE  1 

NONE 

NONE 

NONE 


COMPONENT  4  CONSTRAINED  INTERFACE  MODAL  MODEL 

NORMAL  MODE  2 

NONE 

NONE 

NONE 


Figure  D-6.  Universal  File  for  Component  4 


S  ^  sv 


COMPONENT  4  CONSTRAINED  INTERFACE  MODAL  MODEL 
CONSTRAINT  MODE  1 


NONE 

NONE 

NONE 


1 

1 

0. 

1 

1. 

2 

666666 

3 

333333 

4 

0. 


COMPONENT  4  CONSTRAINED  INTERFACE  MODAL  MODEL 
CONSTRAINT  MODE  2 


NONE 

NONE 

NONE 


1 

1 

0. 

1 

0. 

2 

333333 

3 

666666 

4 

1. 


142  (Kcc  partition  of  constraint  mode  stiffness  matrix) 
2  3  2  2  1 

1  3  2  2  3  1 

.333333  -.333333  -.333333  .333333 

-1 


Figure  D-6  (continued) 
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-1 

250 

134  (Mcc  partition  of  constraint  mode  mass  matrix) 


2 

3 

2 

2  1 

1 

3 

2 

2  3 

1 

2. 

-1 

-1 

250 

mu 

.888889 

.888889 

2.11111 

132 

(Knc 

2 

partition 

3 

of  constraint  mode  mass 
2  2  1 

matrix) 

1 

1.0 

3 

1.0 

2 

2  3 

.33333 

1 

-.33333 

-1 


Figure  D-6  (continued) 


D-14 


LAUNCH  VEHICLE  SYSTEM  DYNAMICS  SYNTHESIS  USING  COMPONENT  MODELS 
COMP.  1:  RESIDUAL  FLEXIBILITY  COMP.  2 : FREE  INTERFACE  NORMAL  MODES 
COMP.  3:  CONSTRAINED  INTERFACE  MODES  COMP.  4  DISCRETE  SPRING-MASS 
4 

1  4 

COMPONENT  1 

TYPE  20 

4 
1 
4 
1 

3 

1  2 

00011020000000000000  .001 

0.  0.  0. 

0. 

DYNAMICS  DATA  2  C0MP1.DAT 
COMPONENT  2 

TYPE  30 

4 

2 

1  4 

1  2 

4 

4 

12  3  4 

0  0  0110000000 

0.  0.  0. 

0. 

DYNAMIC  DATA  2  C0MP2.DAT 
COMPONENT  3 
TYPE  10 

4 

2 

1  4 

2  3 
0 

0 

0 

2  1000000000 

0.  0.  0. 

0. 

DYNAMICS  DATA  2  COMP3.DAT 
COMPONENT  4 

TYPE  40 

4 
2 


Figure  D-7.  System  Input  File 
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0  2 

0. 


DYNAMIC  DATA  2  C0MP.DAT 


Figure  D-7  (continued) 


